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Molecular Emission Lines 
From the 

Orion Molecular Cloudy 

Paul J. Viscuso, Ph.D. 

Cornell University 1986 

The Orion Nebula has been the subject of intensive study for over one hun- 
dred years. Far-infrared (FIR) molecular line observations of CO in the shock region 
surrounding the infrared source IRc2 have suggested that the molecular hydrogen 
density in the shocked and post-shock gas is roughly 3 x 10® cm -3 . The temperature 
of this gas is on the order of 750-2000 K. IRc2, like other nearby infrared sources 
within the Nebula, is thought to be a site of recent star formation. This object 
is apparently at the center of a massive bipolar molecular outflow of gas, which is 
producing a shock front where it meets the ambient molecular cloud surrounding 
IRc2. Study of such regions is important for the understanding of the chemical and 
physical processes that are involved in the formation of stars from molecular clouds. 

Recently, several far-infrared transitions among the low-lying levels of OH 
have been observed toward IRc2. OH is thought to be abundant, and it plays an 
important role in the chemical evolution of the shock and post-shock regions. The 
OH emission serves as a sensitive probe of the temperature and density for the 
shock-processed gas. 

A rigorous treatment of the radiative transfer of these measured transitions is 
performed using the escape probability formalism. From this analysis, we determine 
the temperature of the OH-emitting region to be on the order of 40 K. This suggests 
that the gas is part of the post-shock gas that has cooled sufficiently, most likely by 
way of radiative cooling by CO. Such cooling from shock temperatures of several 
thousand degrees can be accomplished in 100 years. A molecular hydrogen density 
of 3 x 10 7 cm -3 and an OH column density of 1.0 x 10 17 cm -2 is found. The beam 
filling factor is determined to be 36%. 



To better understand the role of OH in the Nebula, more observations will 
have to be made at higher spectral and spatial resolution. A continuously tunable 
grating Fabry-Perot spectrometer has been built to achieve spectral resolution of 
15 km s -1 . Instruments such as these will help in determining the dynamics of the 
OH gas from line profiles. Higher spatial resolution must await the next generation 
of infrared telescopes..-^, __ 
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Nel mezzo del cammin di nostra vita 
mi ritrovai per una selva oscura 
che la diritta via era smarrita. 

Ahi quanto a dir qual era e cosa dura 
esta selva selvaggia e aspra e forte 
che nel pensier rinova la paura! 

Dante Alighieri, Inferno 


Midway in the journey of our life 
I found myself in a dark wood, 
for the straight way was lost. 

Ah, how hard it is to tell what that 
wood was, wild, rugged, harsh; 
the very thought of it renews the fear! 

Dante Alighieri, Inferno 
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I. Introduction 


The infrared part of the electromagnetic spectrum spans three orders of 
magnitude in wavelength, from 1 to 1000 microns (/im)-.- We may roughly divide 
this range into three regimes: near-infrared (1-10 /xm), mid-infrared (10-100 nm), 
and far-infrared ( 100-1000 nm). The portion of the far-infrared (FIR) from 300 
to 1000 nm is also known as the sub-millimeter regime. Heterodyne techniques 
are usually employed in sub-millimeter work, whereas incoherent instruments such 
as grating spectrometers, interferometers, and Fabry-Perot spectrometers are the 
mainstay for observing the rest of the infrared spectrum. In chapter II, we discuss 
the instrumentation used to collect most of the data presented in this thesis. The 
optical hardware includes our grating spectrometer and lamellar-grating interferom- 
eter. In chapter III, we introduce our novel, continuously tunable, hinged grating 
Fabry-Perot spectrometer. 

Continuum emission in the infrared is observed from a large variety of astro- 
nomical sources, including Hn regions, molecular clouds, planets, comets, circum- 
stellar shells, extended stellar atmospheres, galaxies, and quasars. This emission 
tells us much about these objects, including such physical parameters as tempera- 
ture and density. Dust particle sizes can also be inferred from such observations. 
Infrared continuum emission studies have also been used to determine the hydrogen 
to helium ratios for in the atmospheres of Jupiter and Saturn. 

Atomic and molecular line emission in the infrared has been observed from 
most of these objects as well. Ionized and neutral atomic transitions, for exam- 
ple, can tell us much about the Hn regions surrounding hot O and B stars and 
white dwarfs associated with planetary nebulae (Shure 1985). Hyperfine atomic 
transitions of ionized carbon and ionized and neutral oxygen are among the prin- 
cipal coolants in the interstellar medium (Stacey 1985). Infrared observations of 
molecular transitions from the gaseous planets such as Jupiter and Saturn have led 
to a better understanding of the gas dynamics of these systems (see, for example, 
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Haas et al. 1985; Tokunaga et al. 1980). 

The visual extinction caused by intervening clouds of dust along the line of 
sight is not a problem in most of the infrared spectrum. We are able to see star 
formation regions embedded in dust clouds using infrared measurements while the 
nascent stars remain optically hidden from view behind their cloak of dust and gas. 

The earth’s atmosphere is opaque over most of the far-infrared spectrum, 
owing to absorption by water vapor and other trace constituents such as CO2 and 
O3. This makes it necessary to make all of the FIR observations at high-altitude. 
The National Aeronautics and Space Administration (NASA) provides two airborne 
observatories for the purposes of infrared and sub-millimeter astronomical research. 
The Lear Jet contains a 12 inch, f/6.5 telescope while the Kuiper Airborne Observa- 
tory (KAO) houses a 91 cm, f/17 telescope. The observations that are discussed in 
this thesis were obtained with the KAO, using our grating spectrometer in conjunc- 
tion with the interferometer to obtain an instrumental resolution, A/AA, of up to 
4000. Stressed gallium-doped germanium detectors cooled to 2 K give us sensitivity 
in the spectral range 100-195 /im. 

One of the most intensely studied star formation regions is the Orion Nebula 
and its associated molecular cloud, OMC-1. Being only 500 parsecs distant, it is also 
one of the closest of such regions. The nebula is thought to be quite young, perhaps 
several tens of thousands of years old, and there are indications that star formation is 
still happening deep within. One of the young stellar objects that has been found is 
an infrared source known as IRc2. This object seems to be at the center of a massive 
bipolar outflow which is pushing into the surrounding ambient molecular cloud, 
giving rise to shock structures on the order of 10 17 cm in size with temperatures in 
the range 750 to 2000 K (Erickson et al. 1982; Watson 1982). The temperatures and 
densities were derived from molecular hydrogen and CO observations (Beckwith et 
al. 1978; this work, chapter V; and references therein). The molecular hydrogen 
density in the shocked and post-shock gas is roughly 3-10 x 10 6 cm~ 3 . 

At these temperatures and densities, the CO is collisionally excited by the 
H2, especially the higher J rotational levels. The optically thin transitions in the 
FIR part of the spectrum are not very sensitive to densities above 10 5 cm -3 and 
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t em peratures above their excitation temperatures. These transitions are thus good 
yardsticks by which we may measure the amount of CO present in the shocked 
gas. We report in chapter V the first astronomical detection of the J — 15 — ► 14 
rotational transition of CO at 173.631 j*m toward OMC-1. The 173 fxux line flux is 
consistent with other CO rotational line fluxes measured and reported previously 
by Stacey et al. (1982), Stacey et al. (1983), and Watson (1982). The data support 
the view that the CO-emitting region is contained within the post-shock gas of 
molecular hydrogen density — 3 x 10® cm -3 and kinetic temperature of — 750 K. 
The mass of hydrogen that is associated with this region is between 1.0 and 1.5 M®. 

Another important molecular constituent of the gas, according to chemical 
equilibrium models, is OH. It is possible to determine the physical parameters of 
the OH-emitting region by using the same line of analysis that was used for the 
CO. We need to collect a data base of as many rotational transitions as is possible, 
and then try to determine the temperature and density of the emitting region using 
conventional radiative transfer theory. To this end, we have made observations of 
the 2 n 3/2 J — \ \ rotational transitions of OH at 84.42 and 84.60 pm and of 

the 2 n 1 y 2 J = § — ► j rotational transitions at 163.12 and 163.40 pm. Coupling 
these data with the 2 II 3 / 2 J = § — < ► \ (119.23 and 119.44 /xm) data of Watson et 
al. (1985), we are able to construct a model of the OH-emitting region in which 
the above data are consistent. This model, discussed in chapter V, makes use of 
the new experimental H 2 -OH collisional cross-sections of Andresen, Hausler, and 
Lulf (1984) and Schinke and Andresen (1984), which are significantly different from 
earlier values calculated by Dewangan and Flower (1981, 1983). 

Since we have measured rotational transitions that involve the 5 lowest rota- 
tional energy levels, the data are not very temperature sensitive for temperatures 
above about 100 K. We cannot, therefore, determine the temperature of the OH- 
emitting region. If the temperature is 100 K, then we axe talking about the post- 
shock gas which will form a column roughly 10 17 cm deep (Genzel 1985). However, 
if the temperature is near 1000 K, then we axe looking at the gas near to or within 
the shock-front. Since we expect an enhancement of about a factor of 40 of the 
OH within a thin sheet about 10 14 cm deep (Draine and Roberge 1982), then for 
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high-temperatures, this thin region of enhanced OH should dominate in the line 
emission. 

The molecular hydrogen density that we derived from the OH data is quite 
insensitive to temperature. Whether the OH is coming from the shock front or from 
the post-shock gas, the density of H 2 in the emitting region is »1.0 x 10 7 cm -3 . This 
compares favorably with the somewhat lower estimate from the CO data (Watson 
1982). New and more accurate H 2 -CO collision rates have been calculated by 
Schinke et al. (1985), and suggest that the older rates overestimated the large AJ 
cross-sections. We expect that the newer rates will require a slightly higher density. 

The model used in chapter V makes use of somewhat simplified and idealized 
assumptions concerning the treatment of radiative transfer. In order to make better 
use of the data, we have constructed a numerical radiative transfer code based upon 
the escape probability formalism first introduced by Sobolev (1957). In chapter VI, 
we discuss the results of this code. We find that the three sets of infrared data (84, 
119, and 163 /*m) can be fitted by a one-component model. For high temperatures 
(T>50K), the observed 163 ^m flux is always too bright, usually by an order of 
magnitude. This rules out the shock front or hot post-shock gas as the source of 
the OH radiation. We find that the temperature of the gas must be about 40 K to 
produce the observed line ratios. The H 2 density is about 4 x 10 7 cm -3 and the OH 
column density is on the order of 10 17 cm -2 . This fits in with the present thinking 
of Melnick (1985) concerning OH in OMC-1. 

The failure of the simple model used in chapter V to predict these low tem- 
peratures is a result of several factors. One reason for the incongruent conclusions of 
the two methods of analysis is that not all of the radiative processes were taken into 
account, and these processes are significant. For example, not every OH molecule 
that is excited into the 2 n 1/ / 2 ladder will decay and emit a 163 /xm photon. There 
are other routes of decay back down to the ground state. Another source of error 
was the inability to reliably estimate the level populations in the simplified model. 
A third effect is that the simple model assumes that photons may be trapped and 
reradiated throughout the cloud, whereas we know that the radiative processes must 
be local phenomena owing to the velocity gradient that exists in OMC-1. 



Whereas the simple method of analysis presented in Chapter V is useful in 
the interpretation and understanding of the physics, it is good only as a first-order 
analysis of the data. For OMC-1, it is necessary to perform more exact numerical 
radiative transfer calculations. These calculations allow us to better determine the 
physical parameters of the OH-emitting region. 

Better spectral and spatial resolution would greatly enhance our ability to 
understand the OH data. We have developed a unique Fabry-Perot spectrometer 
which will allow us to have spectral resolutions as small as 15 km s -1 . Whereas this 
resolution is not unique to FP systems in and of themselves, the novel aspect of 
this instrument is that a grating is placed between the two highly reflecting plates. 
With this arrangement and proper placement of the two plates and the grating, 
it is possible to continuously tune the instrument over the range of 138-170 ^m. 
This tuning is achieved by a simple rotation of the grating. The plates never move, 
thus making the system a very sturdy and stable one that is well suited to space 
missions. This hinged grating Fabry-Perot is discussed in chapter HI. 

Finally, in chapter VII, we list what needs to be done for future work in order 
to resolve the OH issue. Also discussed are plans for redoing the CO calculations 
of McKee et al. (1982) and Watson (1982) using more accurate collisional cross- 
sections with H 2 computed by Schinke et al. (1985). It is expected that these new 
models may result in an upward revision of the temperature of the CO-emitting 
region, owing to the smaller collision rates for large A J. Future plans for using the 
new Fabry-Perot spectrometer are also discussed. 



II. Instrumentation 


2.1 Introduction 

The detection of FIR radiation is essential to the study of a wide variety of 
celestial objects. The central engines powering quasars could release much energy in 
the FIR. Abundances of neutral and ionized atomic species and physical conditions 
of the Galactic center can be deduced from FIR observations. This can also be done 
for active galactic nuclei and starburst galaxies such as M82 and NGC 1068 (Craw- 
ford et al. 1985). For example, recent infrared observations of Nell emission from 
the central 2 parsecs of the Galaxy revealed evidence for a massive black hole in the 
center (Serabyn and Lacy 1985). Fine-structure lines from highly ionized species 
cam determine temperatures amd densities in the ionized regions of plametary neb- 
ulae amd Hll regions (Shure 1985; Herter, Heifer, amd Pipher 1983). Temperature 
profiles of planetary atmospheres cam be deduced from FIR continuum measure- 
ments. Emission data of molecular amd atomic tramsitions in the FIR cam tell us 
about the physical conditions in the shock regions associated with mass outflow 
from nascent stars. Circumstellar dust shells cam best be studied in the infrared. 
Galaxies that were formed in the early universe will also have much of their energy 
output redshifted into the infrared pant of the spectrum. Virtually the whole gamut 
of objects that are known to exist cam be seriously studied in the far-infrared. 

Although broad-band photometry, such as that done by the infrared astro- 
nomical satellite (IRAS), is important in these studies, it is necessauy to have high 
spectral resolution in order to gain a deeper understanding of the physical condi- 
tions within these regions. For example, the O I 3 Pj — ► 3 P 2 fine structure transition 
at 63.17#im is one of the major coolants in the interstellar medium (ISM). With 
broadband photometry, it is impossible to deduce the contribution of this line to the 
astrophysics of the ISM. However, with high spectral resolution instruments, this 
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line can be separated from the continuum. The instruments which can provide us 
with this kind of resolution are grating spectrometers; Michelson, lamellar-grating, 
and Fabry-Perot interferometers; and heterodyne receivers. 

The energy of a FIR photon is —0.01 eV, two orders of magnitude below the 
average energy involved in chemical reactions. For this reason, no photographic 
emulsions are available that can detect FIR photons. The most efficient method 
of detecting these photons is the use of extrinsic photoconductive material such as 
gallium-doped germanium (Ge:Ga). More will be said about detectors later. 

2.2 Why Airborne Astronomy? 

2.2.1 Atmospheric Transmission 

The detection of FIR radiation is hampered by the water vapor in the at- 
mosphere. Carbon dioxide, CO2, and ozone, O3, play a somewhat lesser role in 
the absorption of FIR photons, as well. In figure 2.1, we see the transmission of 
the atmosphere for wavelengths between 120 and 170 nm at 4 different altitudes 
(Traub and Stier 1976). At an altitude equivalent to that of the top of Mauna Kea 
in Hawaii (4 km), we see that the transmission is ~0% for nearly the entire range 
of wavelengths. This makes ground-based work in the FIR virtually impossible. 

As we move up into the atmosphere, the transmission begins to greatly im- 
prove. For unmanned balloon packages at 28 km, there are rather narrow atmo- 
spheric features, but most of the spectrum is relatively free from absorption. For 
manned observations, the highest altitude available for reasonably sized instruments 
is that reachable by aircraft. It is important to also be above the troposphere, the 
section of the atmosphere closest to the earth’s surface. It is in the troposphere 
where convection occurs and, hence, weather. Due to the convection, the water 
vapor content at the higher altitudes in the troposphere can still be quite high. 
The tropopause, the top of the troposphere, is usually at an altitude between 33000 
and 41000 feet. An altitude of 14 km provides us with many clear windows in the 
FIR. Pressure-broadened water vapor absorption features, however, do still provide 
a significant amount of blockage of FIR photons from astronomical sources. 
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Wavenumber 


FIGURE 2.1 Atmospheric transmission profiles for four altitudes: 41 km, 28 km, 
14 km, and 4 km. From Traub and Stier (1976). 
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The advantages of airborne packages over balloon-borne telescopes are better 
pointing and target-tracking accuracy and the control that the observer actually has 
over the observations. One of the disadvantages is that at the lower altitudes of 
airplanes such as the Kuiper Airborne Observatory (KAO), one is still bothered by 
telluric features and, at times, also by clouds. 

2.2.2 Airplane Environment 

Since we must make use of an observing platform that can penetrate the 
stratosphere, we use the KAO. The KAO is a Lockheed L300 configured as a pro- 
totype version C-141A Starlifter jet transport that is operated by NASA/ Ames 
Research Center from Moffett Field Naval Air Station in California. Sitting in a 
cavity in the port side of the forward section of the plane is a 0.91 m Cassegrain 
telescope. The optics constitute an f/ 17 system. 

The KAO is large enough to comfortably accomodate a retinue of observers 
and staff onboard. Six or seven observers can be on a flight (as opposed to only 
two on NASA’s Leu Jet; Stacey 1985). Five flight personnel are also on board to 
operate the telescope and associated support systems. They include the following. 
The mission director coordinates interaction of the KAO staff, the investigator team, 
and the flight crew. He also assures optimum support of research activities. The 
telescope operator activates and controls telescope stabilization while the systems 
monitor is responsible for the proper operation of the oscillating secondary mirror 
and the telescope cavity. The Airborne Data Acquisition and Management System 
(ADAMS) is run by the ADAMS operator. The tracking operator acquires the 
astronomical objects of interest in the closed-circuit video system which allows us 
to guide on the object or suitable nearby guide stars. Since the telescope has a 
limited range of motion in azimuth (±2°), it is necessary for the plane to fly in a 
specified direction in order to observe a certain object at a specific time. For this 
reason, a flight plan which lays out the routes that must be taken in order to observe 
the objects of interest must be worked out and filed at least one day before flight. 
The flights generally begin and end at Moffett Field and last up to 7.5 hours, of 
which 6.5 hours can be dedicated to observing. Optimizing the flight plans involves 
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rearranging the observing sequence of sources, choosing the latitude and longitude 
where the observations can start (within reasonable distance of Moffett, that is), 
and choosing the time of takeoff which gives best results. Other considerations to 
take into account in the optimization of the flight plan are several restricted military 
warning areas over which flying could be hazardous that need to be avoided as well 
as the added bonus that the winds that the aircraft will encounter are generally 
not very well known one or more days before takeoff. These could alter the flight 
plan significantly if they are drastically different from the assumed winds. All of 
these factors add up to making the optimization of a flight plan a highly evolved 
art. NASA provides the users with several navigators for just this purpose. 

2.2.3 The Telescope 

The telescope that is onboard the KAO consists of a primary paraboloid of 
diameter 36 inches (91.5 cm). The effective focal length is 610 inches yielding a 
focal ratio of f/ 17. The diffraction limit at 150 /im is roughly 6.2 mm. At a scale of 
14 /; /mm at the focal plane, a 4 x 4 mm aperture corresponds to a » l' x l' beam 
size. 


The telescope is supported in the airplane by an air bearing. This spheri- 
cal bearing allows the telescope to float with a minimum amount of friction and 
influence from any airplane rotations. The system is gyro-stabilized to allow ac- 
curate tracking of celestial objects. Most investigator packages use the Nasmyth 
instrument cavity. Once the investigator’s instruments are attached to the Nas- 
myth cavity, the telescope must be balanced to insure the proper operation of the 
gyros. 

Since the telescope is in contact with the ambient air at altitude, this will 
cool down the telescope. It has been found that an image shift of up to 2 1 will result 
from this cooling. In order to minimize the shift during flight, it is possible to pre- 
cool the telescope before flight to a temperature of -30° C. In reality, the pre-cool 
may get the temperature of the primary mirror down to 0°C. The efficiency of the 
pre-cool is dependent on ambient temperature and relative humidity. 
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2.2.4 The Chopper 

At an altitude of 14 km and a wavelength of 150 ^m, the emissivity of the 
atmosphere is 0.13 (Traub and Stier 1976). The ambient temperature is -50° C, 
or 223 K. With an instrumental resolution of 200 (A/AA), this corresponds to an 
incident flux of 2.8 x 10 -8 W cm -2 str -1 . This is at least a factor of 5 stronger than 
any flux we may hope to receive from the most luminous astronomical sources. Not 
only do we see the atmosphere, but the telescope is also a source of FIR photons. 
All of these background sources will dominate the signal. 

It is desirable to subtract out this background, and this can best be done 
by a method called ‘chopping’. Chopping consists of switching the beam from on- 
source to off-source and then to on-source again. This is usually done at a frequency 
of 25-35 Hz since the fluctuations in the background can occur on relatively short 
time-scales. Chopping is accomplished by having the secondary mirror of the tele- 
scope wobble. In chopping, the on-source beam will have the ordinary background 
contribution plus that of the source. The off-source beam will only have the back- 
ground. By subtracting the two, we find the source contribution. For the KAO, the 
chopper throw can be as large as 15 ; , although for amplitudes of this size, the chop- 
per efficiency (ratio of cycle time in the beams to the cycle time between beams) 
decreases. For most sources a throw of 4-6 ; is sufficient. 

2.2.5 Tracking and Guiding 

Once the object to be studied is acquired in the field of view of the telescope, 
it must be tracked. This is done using the KAO tracker camera and tracking 
computer. The pointing accuracy of the camera is about ±5". Since most FIR 
sources are not in themselves optical sources, it is necessary to use a guiding system 
which enables one to guide on a nearby star while the telescope boresight is looking 
at the FIR source. We have used Professor James R. Houck’s offset-guider coupled 
to the KAO tracking system in order to track objects that we have studied. The 
mean guiding error is about 5". The relative positions of the appropriate guide 
stars to the sources are calculated before flight. 
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2.2.6 Data Acquisition 

The data acquisition software that is run on our IMSAI computer was de- 
veloped by Charles Fuller at Cornell University. The software package will run 
a motor drive box which controls the spectrometer grating and the interferome- 
ter movable mirror position. The program also will record and average the signal 
strengths. The drive box was constructed by Scott Smyers in Ithaca. Scott also 
wrote the fast-fourier transform package that is used to generate spectra from the 
interferometer data. 

Figure 2.2 shows a schematic layout of our system as it is interfaced with the 
KAO. The instrumentation is attached to the telescope and controlled by the data 
acquisition software. The telescope position is controlled by the tracking computer. 
The IMSAI acknowledges when the tracking computer has lost track of the guide 
star and, if this occurs, will stop taking data until the tracking computer has reac- 
quired the guide star. The IMSAI is also interfaced with the ADAMS computer. 
From ADAMS, it acquires such housekeeping data as the object name, rotation 
angle, boresight water vapor readings (in microns of precipitable H 2 O), universal 
time, etc. 

The chopped signal from the detectors is demodulated by Ithaco lock-in am- 
plifiers. These amplifiers use the telescope secondary driver as a frequency reference. 
The output of the Ithacos is an analog voltage signal between —10 and +10 volts. 
This signal is then routed to an analog-to-digital converter board in the IMSAI 
computer. The data acquisition program will sample this input to determine the 
signal strength. The data are sent to ADAMS to be stored on tape. It is also 
written to floppy disks for post-flight reduction on the IMSAI. 

2.3 Detectors 

2.3.1 Overview 

There are several different types of detectors that may be used for infrared 
astronomy. Bratt (1977) gives a detailed account of the physics of infrared detec- 
tors. A short summary follows here. These may be separated into two classes: 
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thermal and quantum detectors. Quantum detectors may be further divided into 
intrinsic and extrinsic detectors. Another division of quantum detectors is into 
photoconductive and photovoltaic detectors. 

Thermal detectors, or bolometers, have been widely used in IR astronomy. 
Bolometers absorb the incident radiation which heats up the crystal. This increase 
in temperature will change a characteristic of the detector, such as its resistance. 
The most popular bolometer is the low-temperature germanium bolometer (Low 
1961). The germanium bolometer is sensitive to radiation from 1 pm to 1 mm. For 
a more detailed discussion of thermal detectors, see Gillett, Dereniak, and Joyce 
(1977) and references therein. 

In quantum detectors, also known as photon detectors, the incident photons 
interact to produce free charges in the semiconductor crystal. There are two modes 
of operation: the photoconductive and the photovoltaic mode. In the photovoltaic 
mode, the freed carriers are separated by the intrinsic electric held that is generated 
by a diode junction. One can then measure a change in one of the diode character- 
istics, which could be a current change at zero bias voltage or a voltage change at 
zero current. An example of photovoltaic detectors is the indium-antimonide (InSb) 
detector used for A < 5.5 pm. 

For photoconductors, two types of electronic excitation are utilized: i) in- 
trinsic excitation and ii) extrinsic excitation. For intrinsic photoconductors, the 
transition occurs in the pure semi-conductor material between the valence band 
and the conduction band. This transition can only occur if the photon has an en- 
ergy greater than the band gap energy, Eg. The temperature of the detector must 
also be below the band gap energy. Examples of intrinsic photoconductors and 
their band gap energies (in microns) are Si (1.2pm), Ge (1.8pm), InSb (5.6pm), 
and Hgi_ x Cd x Te (8-15 pm) (Lena 1978). 

Most of the observations that are reported in this work have been made using 
extrinsic Ge:G a. Extrinsic photoconductors consist of semiconductor material in 
which impurities are implanted during the crystal growth. Germanium or silicon 
usually serves as the matrix into which either group Ilia elements (B, Al, Ga, In, 
or Tl) or group Va elements (P, As, Sb, or Bi) we implanted. In the crystal lattice, 
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one of the dopant atoms will replace one of the matrix atoms. 

2.3.2 Ge:Ga Detectors 

Let us consider our Ge:Ga detectors. Since Ge has 4 valence electrons, it will 
form covalent bonds with 4 neighboring Ge atoms by sharing an electron. At an 
impurity site, the Ga has only 3 valence electrons. (The Ga is called an acceptor 
impurity since it can accept an electron. Group Va elements have 5 valence electrons 
and are called donors.) With only 3 valence electrons, there is one electron missing 
that would complete the bonding arrangement. In effect, a hole is created which is 
loosely bound to the impurity site. 

The existence of impurities actually introduces discrete levels just above the 
valence band. When a photon is absorbed by the detector, an electron is excited 
from the valence band into one of these impurity acceptor levels. Alternatively, a 
hole is excited into the valence band. When an electric field is applied across the 
detector, the hole can now move until it recombines at another impurity site. The 
number of passes that the hole makes through the detector before it recombines is 
known as the photoconductive gain, Gp. 

The fraction of incident photons which actually excite the charge carriers in 
the detector is known as the quantum efficiency, rj. It is desirable to have a high 
quantum efficiency; viz., rj « 1. This requires having a high dopant level; however, 
if the dopant level is too high, then the energy levels due to each impurity atom will 
begin to overlap. This phenomenon is known as impurity banding. With too many 
impurity sites, the photoconductive gain will decrease since the mean free path for 
the charge carriers will be small and the resulting responsivity will also be small. 
Therefore, the doping must be kept at a low level. Since the quantum efficiency, 
and hence the cross-section for absorption of a FIR photon, will be relatively small 
with the low doping levels, it is advantageous to have a larger detector crystal. 
With typical doping levels, the mean free path of an impurity site ionizing photon 
is roughly 5 mm . This means that detectors should be on the order of 5 mm in size. 
Unfortunately, detectors are rather susceptible to cosmic rays (Houck and Briotta 
1982). The bigger the detector, the larger the cross-section for absorption of cosmic 
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rays. For this reason, detector dimensions are kept small. 

The dimensions of our Ge:Ga detectors are 1 x 1 x 3 mm. Since these sizes 
do not provide for adequate opacity to FIR photons, it is necessary to place the 
detectors in integrating cavities. These cavities allow the incoming photons to 
be trapped and pass through the detector until most have been absorbed by the 
detector. 

The acceptor impurity energy levels for Ge:Ga are about 10~ 2 eV above 
the valence band with doping levels of ~ 2 x 10 14 cm -3 . This corresponds to a 
cutoff wavelength of about 120 /xm. Price (1961) and Hall (1962) have reported 
that uniaxially stressing the detector along the [100] crystallographic axis splits the 
four-fold degeneracy into two doubly degenerate bands, one with slightly higher 
energy and one with slightly lower energy. The splitting is proportional to the 
applied stress, and a stress of 6.6xl0 7 kgf m -2 will result in a gap between the 
valence band and the acceptor levels of 6meV (Haller, Hueschen, and Richards 1979; 
Kazanskii, Richards, and Haller 1977). This results in a photoconductive response 
out to wavelengths of »205 /xm. We have achieved sensitivity out to roughly 195 /xm 
with our detectors. It should be noted that for low signal levels, Ge:Ga used as a 
photoconductor is 2 or 3 orders of magnitude more sensitive at these wavelengths 
than bolometers. 

Our detector material was ordered from Eagle-Picher in the form of 2.75" 
diameter wafers. These wafers were cut 3mm thick along the [100] axis. The ger- 
manium was doped with 2.1 x 10 -14 atoms cm -3 of gallium. The residual impurity 
level was less than 10 13 cm -3 . Haller, Hueschen, and Richards (1979) describe the 
process of making ohmic contacts. Boron ions are implanted into the two faces of 
the wafer, first at a surface density of 10 14 ions cm -2 at an energy of 25 keV, and 
then at a dosage of 2 x 10 14 ions cm' 2 at 50keV. Next, the wafers are annealed 
in an argon atmosphere at 200 °C in order to remove implantation damage. This 
process also activates the boron. The faces are then wetted with 500 A of palladium 
upon which is deposited a 1 /xm layer of gold for electrical contact. 

Since a stress needs to be applied to the detector, the detector must be 
mounted in a housing. We also need an integrating cavity, so the housing will also 
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serve this purpose. The detector housings are made of hard brass. This material 
has the hardness that is needed in order to minimize the cold flow that would be 
associated with aluminum housings (Watson 1982) and also has a high coefficient 
of thermal expansion. Most of the stress that a detector experiences is due to the 
differential thermal contraction between the housing and the detector upon being 
cooled down to liquid helium temperatures (4.2 K). Hence, the housing must be 
made of a material that will undergo a substantial contraction during cooldown. 
Stacey (1985) discusses in detail the manufacture and design of housings and detec- 
tors. Other descriptions of stressed detector construction can be found in Watson 
(1982) and in Haller, Hueschen, and Richards (1979). 

The first astronomical detection made with stressed gallium-doped germa- 
nium detectors was in 1980 (Russell et al. 1980). This observation was of the 157 pi n 
C n line in Orion. This line is one of the dominant coolants in the ISM. 


2.3.3 Detector Performance 


The total number of photons absorbed per unit time interval can be written 


as 


XT _ I v C\AT v rj v 6u _ x 
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( 2 . 1 ) 


where I u is the specific intensity in erg s -1 cm -2 str -1 Hz -1 , Cl is the solid angle of 
radiation incident upon the detector, A is the collecting area upon which the radi- 
ation is incident, r v is the transmission of the optics, rj u is the quantum efficiency, 
and 6u is the bandpass. In a diffraction limited system, the product Afl is equal 
to 4.7A 2 . If we use the system resolution, R, in place of the bandpass, then we can 
rewrite (2.1) to be 


N u 


4.7 1%/ A 2 Ty TJy 
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( 2 . 2 ) 


For the present, let us consider only the background radiation incident upon 
the detectors. In this case, at an altitude of 14 km, the detector will be looking at 
radiation coming from the earth’s atmosphere, which behaves like a grey body with 
temperature T v and emissivity e„. With an emissivity of 0.13 at 150 /zm (Traub 
and Stier 1976) and a temperature of about 220 K, the number of photons absorbed 
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by the detector is given as 

N„ ~ 2 x 10 11 R -1 photons s _1 , (2.3) 

where we have assumed an instrument transmission of roughly 30% and a quantum 
efficiency of 0.15 (Stacey 1985). The fluctuation in the number of detected photons 
from the background radiation is given as 

AN V = [2Ni/B(l , (2*4) 

where B is the electrical bandwidth and is the photon occupancy which is simply 
( e WkT — i)~ l (Boyd 1982). The first term on the right side of (2.4) is known as 
‘shot noise’ and is due to the uncorrelated arrival rate of photons. The second term 
known as ‘excess noise’, and is due to the correlated bunching of photons. This 
excess noise is important when hu < kT. Using (2.2) in (2.4), we find that 

A N y « 3 x 10® [^^(1 + 1 . 8 r 7 p)]I . (2.5) 

It is clear that in order to decrease the photon fluctuations which contributes 
to the noise of the system, we can either decrease the throughput, decrease the quern- 
turn efficiency, or increase the resolution. Decreasing the throughput will have the 
effect of also decreasing the signal strength. Since the signal strength is proportional 
to r„ and the noise is proportional to the square root of r„ , the signal-to- noise ratio 
would be proportional to v . Decreasing the throughput will therefore degrade 
the overall performance of the system. The same argument may be used for the 
quantum efficiency. It is only by increasing R that we can improve the system. This 
assumes that photon noise dominates all other sources of noise. 

Other sources of noise include phonon-generated noise in the detector, back- 
ground radiation from warm filters and other optical parts in the light path, and 
electronic noise from the rest of the circuit. Phonon noise can be minimized by 
having high-purity detector material and by also cooling the detector down to 2 K 
or less. With cold optics and filters, the background emission will be limited to that 
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caused by the atmosphere. With careful wiring and design, the electronic noise can 
also be minimized. In theory, one should be able to achieve background limited 
photon noise. 

During actual observations, we have found that the background fluctuations 
are on the order of 5 or 10 times larger than those due to photon noise. These fluc- 
tuations are probably caused by a combination of changes in the air mass, changes 
in the sky temperature and emissivity, and turbulence inside or outside the tele- 
scope cavity. It is still the case that higher resolution will lower the background 
noise. For more discussions on background-limited detections in the FIR, see Boyd 
(1982); Gillett, Dereniak, and Joyce (1977); and Keyes and Quist (1970). 

2.4 Dewars 

Since it is advantageous to have the optics and the detectors at temperatures 
close to or below that of liquid helium (4 K), we place the spectrometer, filters, and 
detectors within a helium-cooled dewar. The dewar consists of a 3.5 liter stainless 
steel cryogen can which has a thick copper block, or work surface, at one end to 
which are affixed the detector housings and the spectrometer. The copper acts as 
a good conductor and provides for adequate cooling of the detectors through their 
brass housings. The thickness of the copper provides for mechanical as well as 
thermal stability. 

Surrounding the cryogen can are about 16 layers of aluminized mylar (super- 
insulation). The insulation reflects most of the radiation coming from warmer parts 
of the dewar and helps to keep the helium from boiling at an unacceptable rate. 
The average hold time for 3.5 liters of liquid helium is roughly 22 hours. 

Surrounding the cryogen can and the spectrometer is a cylindrical shell of 
stainless steel that is wrapped with another 30 layers of super-insulation. This is 
then all placed inside an outer can which is evacuated. The cylindrical shell provides 
a radiation shield from the room temperature radiation emitted by the outer can. 
Evacuating the dewar minimizes the heat path from the cryogen to the outside 
world. 

Inside the dewar is placed activated charcoal wrapped in aluminum foil with 
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holes punctured in it. This serves as the getter, which absorbs the heavier molecules 
left in the vacuum jacket of the dewar and which are not cryogenically frozen out 
at liquid nitrogen temperatures (77 K). The outer can must be connected to the 
helium can in some fashion. At the spectrometer end, the connection is via a 
double-triangular fiberglass system to the outer can. The fiberglass offers low heat 
conductivity as well as provides for several points of stabilization. On the cryogen 
can end, the connection to the outside world is made via a thin-walled neck. This 
neck is made of 0.010 11 stainless steel and welded onto the top of the can. In years 
past, it was attached with a silver-solder joint. However, these joints began to fail 
regularly several weeks before a flight series and we have since gone to using heli-arc 
welds. Most likely, with the silver-solder joints, not all of the flux was cleaned out 
of the joint and it slowly corroded. The weakened joints would then not be able to 
handle the thermal shock sustained upon cooldown. The welds are stronger and do 
not require the use of flux. 

The cooldown process cam take anywhere from 12 to 24 hours. Once the 
dewar is assembled, the vacuum jacket needs to be evacuated. Gases will adsorb 
onto the superinsulation, making it necessary to pump on the jacket for up to 10 
hours. New cement or paint may also outgas for several hours as well. Once the 
pressure in the jacket is down to 30 millitorr or less, the cryogen can is filled with 
liquid nitrogen. The nitrogen will cryo-pump out most of the rest of the gas in the 
vacuum jacket. The nitrogen also provides a cheap and effective way to cool down 
the spectrometer and cryogen can as far as possible before using liquid helium. 
Helium prices are about thirty times higher than those for nitrogen. Once thermal 
equilibrium is reached in about 2 hours, we then pump on the nitrogen in the dewar. 
By lowering the ambient pressure over the liquid nitrogen to about 100 torr, we can 
cool the nitrogen down to roughly 55 K. At this point, the nitrogen is frozen solid. 
To effectively decrease the temperature of the spectrometer, it is necessary to have 
the dewar holding nitrogen for at least 6 hours. 

After this time, the solid nitrogen can be partially melted by the addition 
of liquid nitrogen. This is then dumped and liquid helium is introduced into the 
cryogen can after one makes sure that all of the nitrogen is out of the can. Otherwise, 
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the nitrogen can form a sheet of ice between the helium and the work surface and 
lower the thermal conductivity enough to adversely affect the performance of the 
detector. It takes the fill and evaporation of about 15 to 20 liters of helium before 
thermal equilibrium is reached. Since Ge:Ga requires temperatures on the order 
of 2 K or less, we need to lower the pressure over the helium in order to lower the 
helium temperature, as was done with the nitrogen. The hold time of the cryogen 
can when we pump on the helium is 14 to 16 hours. 

2.5 Filters 

Stressed gallium-doped germanium detector material is sensitive to radiation 
of A < 200 pm. This includes photons in the mid- and near-infrared regimes, as well 
as in the visible part of the spectrum. The grating spectrometer is capable of passing 
a relatively narrow bandwidth (A/ A A ss 210). However, if the grating is efficiently 
diffracting radiation of wavelength A, then it will also efficiently diffract radiation of 
wavelength A/2, A/3, and so on. It is desirable to filter out the unwanted radiation 
that can be diffracted. This is done by using several different types of filters. 

The first optical material that the incident radiation encounters is a win- 
dow on the side of the dewar. The window material is composed of 0.25mm thick 
polyethylene, which is transparent to radiation from about 16 fim well into the 
millimetric regime. Polyethylene serves as a rugged window material that is non- 
permeable (except slightly to helium gas) and non-soluble in water, hence preserving 
the vacuum inside the dewax. The polyethylene window is at room temperature. 

The next filter that the incident beam encounters is a diamond dust scattering 
filter which is placed in the aperture of the spectrometer. This filter is composed 
of a monolayer of diamond dust that is pressed into melted polyethylene. The 
diamond dust is 8-16 /xm in size and is purchased from the Dunnington Comp., 
Chester Springs, PA. The resulting filter is 0.15 mm thick and 0.5" in diameter. 
The dust scatters visible and near infrared radiation out of the beam, and works 
best when it is far away from the detector and in a slow (high focal ratio) beam. 
See Watson (1982) for more details on diamond dust scattering filters and their 
construction. 
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The radiation that goes on unimpeded by the diamond dust scattering filter 
will pass through a 1mm thick BaF 2 crystal. The BaF 2 has a very high reflectivity 
for wavelengths between 18 and 58 nm (Armstrong and Low 1973). This band 
corresponds to the ionic vibration band within the crystal. Alkali crystal filters 
that display this phenomenon are known as reststrahlen (residual rays) filters. 

The index of refraction for the alkali crystals is ~ 2.4. With such a high 
index, the reflection losses at the interfaces are ~ 30%. In order to minimize these 
losses, the filters are coated with a thin layer of polyethylene (n; ss 1.5). This allows 
roughly 81% of the incident radiation to pass through. The BaF 2 filters transmit 
~ 90% of the radiation at 157 /im. 

Another problem of these filters is the Fabry-Perot action that results if the 
faces are parallel. This is not a particularly welcome attribute to the system, since 
this action will allow only narrow bands through. This can be avoided by hand 
polishing the faces of the crystals. If the faces are wedged or non-planar, which is 
the usual result from hand polishing, then the Fabry-Perot action will be minimized. 

The last filter through which the radiation passes before it goes to the detec- 
tor is another reststrahlen filter of AgCl that is about 0.5 mm thick. The band stop 
for this filter is from 42-100 /im and it transmits ~80% at 157 fim. All of the filters 
are cooled down to liquid helium temperature in order to minimize the amount 
of background radiation introduced into the beam. It has also been observed by 
Hadni et al. (1965) that the cut-on and cut-off wavelengths for a reststrahlen filter 
axe much sharper at 4 K than at room temperature. This effect is due to the change 
in phonon population for kT <C hi/. 

Other reststrahlen filters and their stop bands are listed in Table 2.1. For 
the Fabry-Perot spectrometer, a film of polyethylene impregnated with carbon is 
used in place of the diamond dust scattering filter. This ‘black poly’ transmits « 
70% at 153 /im. The cut-on wavelength for the black poly depends on the density 
and size of the carbon particles. For the material that we have used, this cut-on is 
around 50 n m. An additional reststrahlen filter of CsBr is used in order to cut out 
radiation shortward of 140 ^m. CsBr has a stop band between 87 and 125 /im at 
4 K. 
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Table 2.1 Infrared Filter Properties at 4 K 

Sample 

Thickness (in mm) 

Tmax 

Stopband(^m) 

Sapphire 

1 

0.59 

7-18, 20-23 

Crystal quartz 

0.86 

0.78 

20-27 

LiF 

2 

0.65 

13-34 

CaF2 

1 

0.68 

13-38 

SrF2 

1 

0.71 

20-45 

BaF2 

1.5 

0.65 

18-58 

KC1 

2.4 

0.8 

37-72 

NaBr 

1.2 

0.88 

42-78 

KBr 

1 

0.53 

52-85 

AgCl 

1 

0.60 

42-100 

KI 

2.2 

0.88 

60-100 

CsBr 

1 

0.68 

87-125 

Csl 

1 

0.65 

115-155 

T1C1 

1.5 

0.70 

52-190 

KRS-5 

1.3 

0.82 

83-300 


Properties of infrared filters at 4K. Included is the thickness of the sample, the 
maximum transmission, and the stop band, that range of wavelengths for which the 
reflectivity is close to unity. Data from Armstrong and Low (1973). 
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Fixed Fabry-Perot etalons can be used as filtering devices. In the FIR, it is 
possible to also use metal meshes. These meshes will transmit or reflect radiation 
with a wavelength on the order of and longer than the spacing on the mesh, depend- 
ing on whether the grid pattern is inductive or capacitive. For more information on 
filters, their construction and properties, see, for example, Stacey (1985), Watson 
(1982), Armstrong and Low (1973), and Kimmitt (1970). 

2.6 Electronics 

The charged carriers that are excited by the absorbed photons in the detector 
will add to the bias current that is already flowing through the detector. This change 
in current is measured by using a transimpedance amplifier (TIA). Figure 2.3 shows 
a schematic of the circuit. The circuit is composed of a source follower which feeds 
into an operational amplifier (op amp). The feedback loop goes to a load resistor 
which is connected to the detector. The source follower is a 2N 6484 junction field- 
effect transistor (JFET). One gate is connected to ground. The source belonging 
to that gate is then connected to the positive feedback pin of an LM308 op amp. 
The other source is connected to the negative feedback pin. The op amp will act 
to keep the voltage at the negative feedback pin the same as the voltage at the 
positive feedback pin. It does this by controlling the voltage output, V ou t. Since 
the JFET is matched, the other gate that is connected to the detector and load 
resistor is at a virtual ground. What this means is that the current flowing through 
the load resistor is equal in magnitude but opposite in direction to the current 
flowing through the detector. This can be written as 

Vb i as _ Vout , 0 

= ^ ' (2 ' 6) 

We see that if is 0, then Vo U t should also be 0. However, the two 

channels on the JFET may not be perfectly matched and phonon and Johnson noise 
will create a dark current flowing through the detector. In order to accomodate for 
this, a variable resistance to ground is placed between the two feedback pins of the 
op amp. By adjusting the variable resistor, we may zero out the dark current. 


Detector 


tOOK 


300 K 
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When FIR radiation is absorbed by the detector, an additional current due 
to the photo-excited charged carriers drifting in the bias potential is produced. In 
order to maintain the virtual ground, an equal but opposite change in current must 
be induced in the load resistor. This is done by varying V ou t. This change in 
V oa t can then be measured. Since the amount of radiation falling onto the detector 
varies periodically with the chopper, the measured change in V ou t is accomplished 
with a de-modulation circuit. We use an Ithaco lock-in analyzer in order to find 
the amplitude of the sinusoidal change in 

Sources of electronic noise in the system include the detector, the load resis- 
tor, the source follower, and the rest of the electronic circuit. Before 1983, G118AL 
MOSFETs were used as source followers. The MOSFETs worked well at liquid 
helium temperatures and so could be placed as close as possible to the detector and 
load resistor. This is important since the resistances of both are on the order of 10 9 
to 1 0 10 Cl. With impedances this high, shunt capacitive effects become important 
and pickup and microphonics will become a problem if wires are too long. Hence, 
the G118AL transistors were used since they could be placed on the work surface 
next to the detector. Unfortunately, the electronic noise from the MOSFETs dom- 
inated the other sources of noise. JFETs, on the other hand, are electrically much 
quieter. The only problem with silicon JFETs is that they do not operate below 
temperatures of roughly 60 K. Since we need the FET to be close to the detector/ 
load resistor, it must be somehow isolated from the work surface and heated. The 
optimum operating temperature is between 100 and 150 K (Lengeler 1974), where 
the transductance is a maximum. 

The method used to isolate and heat the JFET is a relatively straightforward 
and simple one inspired by George Gull. In the spectrometer, there are two detectors 
and hence, the need for two source followers. A FET is mounted with rubber cement 
at the end of an ‘L’ section made of thin-walled 1/4" fiberglass rods. The fiberglass 
is mounted to the work surface as close to the detectors as possible while keeping 
the FET outside of the spectrometer. This is necessary due to the requirement that 
the FETs be heated and the copious infrared radiation produced by these relatively 
warm objects could dominate the background radiation field and hence the system 
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noise. The heating is accomplished by using a 1.1 kH Allen-Bradley carbon resistor. 
A copper wire is wound around the body of the resistor and then soldered to the 
two drains on the FET. The wire and resistor are enshrouded in a blanket of rubber 
cement to insure good thermal contact. A current of 4 to 5 milliamperes is sufficient 
to keep the FETs at operating temperature. 

The fiberglass posts are notched in order to keep the heat conductivity low 
between the work surface and the JFET. The wires from the JFETs are wound 
tightly around the fiberglass posts and tacked down with rubber cement in order to 
minimize the microphonics and to provide a long thermal path between the JFETs 
and the cold work surface. 

The rest of the circuitry which is outside the dewar contributes negligible 
amounts of noise to the system. I have incorporated an electronic voltage source for 
the bias voltage. Before 1983, the bias voltage was supplied by batteries which were 
capable of delivering a very stable and electrically quiet potential. The electronic 
bias voltage source made from off-the-shelf components has been found to also be 
quite stable. It is buffered by an RC circuit with a time constant of 10 seconds. This 
insures against high frequency oscillations that may be induced by vibrations in the 
aircraft structure during operation. The measured noise of the bias source is about 
15x10-® V Hz -1 / 2 at a frequency of 35 Hz. This is comparable to that measured 
for the batteries. The 60 Hz pickup is minimal in the circuit. One significant 
advantage of the electronic bias source is that it is maintenance free. One need not 
worry about forgetting to check the batteries. The noise contributed by the rest of 
the circuit is on the same order as that of the bias source. 

The load resistor will contribute to the measured noise of the system. The 
noise due to the resistor is known as Johnson noise, and its source is the random 
thermal motions of the electrons in the conducting material (Nyquist 1928). It is 
given by the formula 

v = (4JfcTRA/)* volt Hz - * , (2.7) 

where k is Boltzmann’s constant, T is the absolute temperature, R is the resistance, 
and A/ is the electrical bandwidth. For a load resistor with R = 3 x 10 9 ft at 
T = 2 K, the Johnson noise will be roughly 6x 10 -7 volt Hz -1 / 2 . This is a factor of 
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40 times the noise from the TIA circuitry. 

The detector will contribute several types of electrical noise. The resistance 
of the detectors at 2 K is on the order of 1000 MO. Johnson noise from the detector 
will therefore be on the order of 6 x 10~ 7 volt Hz -1 / 2 . Ohmic contact noise due to 
the interface between the semiconductor material and the gold contacts on the faces 
of the crystal can be significant. Advances in the manufacturing of the detectors has 
led to significant reduction in the contact noise. Impurities in the detector material 
will also contribute to the noise. Impurities within the material we have used are 
less than 10“ 13 cm -3 . Haller has achieved impurity densities one to two orders of 
magnitude less than this in recent years (Stacey 1986, private communication) . At 
a bias potential of zero volts, we have been able to measure Johnson noise. 

The wires that carry small currents to high resistance components are sus- 
ceptible to microphonic pickup. This is a shunt capacitance problem for which the 
most simple solution is to tie or glue down the wires to the work surface so that they 
are not going to move. We use a combination of Kimwipe tissue and rubber cement 
to glue the wires down in an orderly fashion on the work surface. The wires that 
connect the work surface to the outside world are run alongside the cryogen can 
and then wrapped around the neck on tissue before being soldered to the connector 
pins. These wires are also glued down to the surfaces they touch. Microphonics 
has not been much of a problem. Most microphonics are high frequency, which are 
filtered out with a low-pass filter that has a 3dB point at 88 Hz. Coaxial cables are 
used to shield the wires that run from the dewar to the pre-amp box. This distance 
should be as short as possible to avoid any unnecessary pickup. 

It was mentioned earlier that cosmic rays can be absorbed by the detectors. 
The absorption of that much energy will momentarily saturate the detector and 
cause a spike in the output voltage. This spike will have a very small rise time due 
to the nature of the event. Houck and Briotta (1982) designed a circuit that would 
look for a large differential rise in voltage and activate a sample- and- hold circuit for 
a period of time that one may adjust. This circuit has been built and employed in 
our pre-amp box. We have used an LF198 sample-hold device in place of the LM-2 
found in the original design due to its superior characteristics. We have chosen a 
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gain of 75 for the output voltage. The in-flight noise has been reduced about 10% 
due to the use of this circuit. 

2.7 The Spectrometer 

2.7.1 Grating Theory 

For spectroscopic work in the far-infrared, we have used a combination of a 
grating spectrometer and a lamellar-grating interferometer. The spectrometer has 
a bandpass of roughly 0.7 /im, yielding an instrumental resolution of « 210. 

Monochromators use gratings in order to disperse the incoming radiation by 
diffraction. If one considers a planar reflection grating, one essentially has a multi- 
mirrored surface. For FIR work, the separation of the rulings on the grating is 
usually on the order of the size of a wavelength, so each ruling acts as the source 
of one ray of light. Consider collimated radiation incident upon the grating. If we 
follow the light to each individual facet, and then follow it as it is diffracted, then 
we see that there will be a phase shift between adjacent rulings. In order for the 
diffracted light to still be in phase when it is collected, this phase shift must be an 
integral number of wavelengths. This constraint leads to the grating equation: 

mA = jjfsin a + sin 0] . (2.8) 

Here, m is the diffraction order, A is the wavelength of interest, g is the spacing 
between rulings, and a and 0 are the incident and diffraction angle relative to the 
grating normal, respectively. One sees that in 0 th order, the angle of diffraction is 
just the reflection angle (a = -/?). This is true regardless of wavelength. Hence, 
in 0 th order, all of the incident light will be diffracted in this direction. That is 
to say, there is no spectral dispersion in 0 1 * 1 order. The amount of radiation that 
is allowed to go into the higher orders is small, making the grating most efficient 
where it does not disperse. 

It is possible, however, to increase the efficiency of one of the orders by 
simply blazing the grating. A blazed grating, instead of consisting of alternate 
reflecting and non-reflecting lines, is composed of a series of steps. The angle of the 
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step relative to the plane of the grating is the blaze angle. A blazed grating will 
still diffract radiation in the same manner as an ordinary grating, but for certain 
wavelengths, most of the power will be directed to an order other than 0. This 
means that the grating will be highly efficient in a non-zero order where the desired 
dispersion is obtained. 

The grating is only efficient in passing the radiation with electric field vectors 
perpendicular (s-plane) to the grating grooves (Loewen 1977). This causes a 50% 
loss of signal that can be measured by the detector. The grating efficiency for 
the s-polarization depends on the blaze angle. Calculations and experiment have 
shown that the efficiency is ~95% for 0.9< A fg < 1.8 if the blaze angle is above 22° 
(Loewen 1977). For smaller blaze angles, the efficiency is not constant in this range, 
but peaks for A/$r ~ 0.6 and then decreases for longer wavelengths. Our grating has 
been blazed at 26° 45 / for 150 nm radiation. The grating was ruled and blazed by 
Bausch & Lomb Co., Rochester NY. 

2.7.2 Spectrometer Design 

The spectrometer is based on a Littrow-mode design that was developed by J. 
R. Houck and Dennis Ward (Houck and Ward 1979). A schematic is shown in figure 
2.4. The original spectrometer operated at f/6.5, matching the Lear Jet telescope 
beam. Melnick (1981) improved the original design and converted it to f/ 17 optics. 
An f/ 17 beam of light enters the collimator from the telescope. The beam passes 
through the entrance aperture at A and encounters optical and near infrared filters 
at B. The FIR radiation continues onward, hitting a 45° mirror at C, and then 
impinges upon the secondary mirror at D. From here, the convex secondary will 
diverge the beam to a size about 3 inches in diameter at the primary mirror, E. 
This is important since the spectral resolution will depend on the number of rulings 
illuminated, and therefore on the beam size when it hits the grating. 

At the primary, the beam is collimated and sent onto the grating to be 
dispersed. In Littrow mode, the angle of incidence and the angle of diffraction are 
the same (or nearly so). The diffracted light will then retrace the incident light’s 
path. When the light reaches the secondary'for the last part of its trip through the 
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system, it will fall upon two small diagonal mirrors at F which send the dispersed 
radiation into two detector integrating cavities. The detectors axe separated by 
about 4 mm, which corresponds to a spectral separation of — 0.8/xm. 

The grating is controlled by a gear system that is connected to a stepping 
motor via a fiberglass rod. The rod is made of fiberglass in order to insure a 
minimum heat load on the cryogenically cooled spectrometer and work surface. 
The amount of torsional backlash that the rod experiences is equivalent to almost 
a tenth of a revolution. We compensate for this when the grating is being driven to 
a certain position. The vacuum inside the dewar is maintained by placing an o-ring 
around the fiberglass shaft where it exits the outer can of the dewar. 

The spectrometer is constructed entirely of stress-relieved aluminum. This 
allows for uniform contraction while cooling to liquid helium temperatures. The 
spectrometer is bolted to the work surface, which is the copper bottom of the 
cryogen can. A thin film of N grease (Apiezon Products Ltd., London) is sandwiched 
between the work surface and the spectrometer flanges in order to provide better 
thermal contact. 


2.7.3 Optics and Limitations 

The primary and secondary glass mirrors that were used in the f/ 17 system 
were replaced in the fall of 1983 with aluminum mirrors in order to acquire better 
stability and ease of operation of the system. Since the spectrometer is bolted to the 
work surface and cooled down, the whole instrument undergoes thermal contraction. 
If the mirrors are made of glass, then the difference between the thermal contraction 
of glass and aluminum will shift the focus of the instrument several centimeters. 
Because the entrance windows to the detectors’ light cones are 4x4 mm, or the size 
of the diffraction lobe, it is important that the focus be at these apertures. This is 
difficult to assure if the focus shifts during cooldown procedures. 

If the whole system is constructed of the same material, then the spectrome- 
ter will undergo uniform thermal contraction during cooldown and the focus of the 
optical system will remain unchanged. That way, one may align the instrument at 
room temperature and accurately place the focus where it should be. 
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The convex secondary mirror is 3.43 cm in diameter with a radius of curvature 
of 5.258 cm. The concave primary mirror is 8.33 cm in diameter with a radius of 
curvature of 27.06 cm. These figures were arrived at after careful consideration of 
baffling and light loss due to the Cassegrain system. Because of the Cassegrain 
design, the central portion of the beam is blocked so that undispersed light from 
the secondary will not fall onto the detectors. For this system, the amount of light 
loss we have due to blockage of the central portion of the beam is 33%. About 90% 
of the stray radiation which comes into the optical path should be blocked. We 
have normally measured stray radiation through the system to be approximately 
3% of the maximum signal. 

The aluminum mirrors are overcoated with a thin layer of SiO to protect them 
from oxidation. The interior walls of the spectrometer are painted with Nextel- 
Black Velvet paint manufactured by 3M. The FIR absorption of this paint is about 
80% (Melnick 1981). Care is taken to cover up any openings with silvered tape so 
that stray light may be kept at a minimum. 

Operation of the spectrometer in Littrow mode has the advantage of a com- 
pact design which can easily fit into a 6 inch diameter housing and also accomodate 
two detectors. If the spectrometer were operated in an off-axis mode whereby the 
dispersed light did not retrace the path of the incident beam, then it would not 
be necessary to cut out the central 33% of the beam. An off-axis system, however, 
does not easily accomodate multiple stressed detectors and fit into a small package. 

2.8 Spectrometer Operation 
2.8.1 Calibration 

Once the dewar is assembled and cooled down, the system must be calibrated. 
The dewar is placed upon a stand that is known as the Mock Lear. The test stand 
can deliver a range of beams, from f/ 6.5 (simulating the Lear Jet telescope) to f/ 17 
(for the KAO optical system) and slower. 

In order to detect FIR radiation, we need to chop. This is accomplished in 
the lab by using an aluminum chopping wheel which allows us to look alternately at 
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reflected 300 K radiation and then at glass spheres that are cooled to liquid nitrogen 
temperature (77 K). The glass is opaque at FIR wavelengths and therefore acts as a 
blackbody source. The amount of radiation falling into the dewar from a blackbody 
of temperature T is given by 


F *~ \4 e hc/XkT^ AR 1 ’ ( 2 * 9 ) 

where A is the aperture area, D is the solid angle of the beam as seen by the 
aperture, and R is the instrumental resolution. For the Mock Lear setup, we have 

F x = 5.2 x 10 -8 ( — . * — 77 ^ ) — Watt , (2.10), 

A V e O «/A 10 o_i e 1.87/Ajoo — \ ) R Ajoq V ' 

where Axoo is the wavelength in units of 100 (im. At 153 ^m, the flux reaching the 
dewar is 2.3 xlO -8 R ~ 1 W. The dewar assembly sits 17 inches from the chopper, in 
a cradle in which it may be twisted and tilted to align it with the beam. Once the 
dewar has been peaked up in this fashion, the electronics may be tested. 

If the bias voltage across the detector is zero, then there is no electric field 
to move the excited charge carriers. As the bias voltage is slowly increased, the 
charges will begin to move further, but will still move slow enough for most of them 
to recombine. As the voltage is increased still further, the response of the detector 
will level off and stay constant until near the breakdown bias voltage. We find that 
the point of leveling off for our detectors is at a bias voltage near 160 mV. 

The gain of the system, G, is defined as the differential change in V ou t with 
respect to V bi u . From (2.6), we see that this is simply the ratio of the resistivity of 
the load resistor to that of the detector. We find that as the bias voltage increases, 
this ratio also increases, which suggests that either there is ohmic heating of the 
detector or that the resistance of the detector is decreasing due to the better mobility 
of the charged carriers (and impurities) at higher voltages. Knowing the specific 
heat of Ge, the temperature change due to ohmic heating is expected to be less them 
10 -4 K. It is therefore most likely that the increased gain is due to the increased 
conductivity of the detector material. 
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It should be remembered that the measured signal is proportional to the 
response of the detector and to the gain. Since the gain increases while the response 
is constant for increasing bias, the signal will also increase with increasing bias. The 
most im portant measure is the signal-to-rioise ratio, S/N. For bias voltages below 
a certain value, the signal increases faster than the noise. Above this value, the 
S/N decreases. We have found that for the Ge:Ga detectors, the optimum bias 
value is around 0.22-0.24 volts. Noises are typically on the order of 5/zV for a 1 Hz 
bandwidth. 

Now that everything has been optimized, the grating is rotated via a stepper 
motor. Use of a gear on a 10-turn potentiometer in the external drive train pro- 
vides us with a solid reference point from which the grating may be accurately and 
repeatedly driven to any desired position. Figure 2.5 shows a typical calibration 
spectrum. The absorption features are mainly due to water vapor in the optical 
path. The central wavelength of these features are well known, and hence knowing 
the position of the stepping motor gives us a calibration of the grating in wave- 
length vs. step position. It should be noted that the absorption feature at 152 pm 
is due to the AgCl filter that is used, and not to any telluric features (Stacey 1985). 
Close, but distinguishable features in the spectrum assure us that the instrumental 
resolution of the spectrometer is around 210. 

Many of the water vapor features are pressure broadened and saturated. Be- 
cause of this, it is sometimes difficult to exactly determine the position of the center 
of the absorption line due to asy mm etries in the wings and probable overlapping of 
lines. It is therefore desirable to have an evacuated chamber in which to calibrate 
the system. To this end, Charles Fuller designed and built an evacuated box which 
provides us with a chopped source, the choice of f/6.5 or f/17 optics, and several 
apertures at the source. The source is also movable in x and y coordinates, which 
allows us to do an aperture map of the spectrometer. The movement of the source 
is controlled by stepper motors and can be done under the direction of the IMSAI 
in one tenth the time that it used to take to do manually. The chamber can be evac- 
uated to ~500 mTorr, which corresponds to roughly 1/1500 atmosphere. Figure 2.6 
shows the spectrum obtained using this chamber. For comparison, Figure 2.7 is a 
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FIGURE 2.5 Laboratory calibration spectrum on the Mock Lear using the FIR spectrometer. The radiation passes 
through roughly 20 inches of air at atmospheric pressure and room temperature. 
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FIGURE 2.6 Same as Figure 2.5, except the air pressure in the optical path is roughly 10 3 atmosphere. 
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spectrum that was taken using a collimated source constructed by Dr. E. Erickson 
of NASA/Ames Research Center and mounted on the KAO secondary. The path 
traveled by the radiation through air is about 3 or 4 meters. One sees how fast the 
increasing water vapor content suppresses the FIR transmission. 

It is necessary to calibrate the spectrometer after every cooldown for several 
reasons. Mechanical couplings are not placed in the same position, and it is therefore 
necessary to account for this. The first diagonal mirror in the spectrometer which 
the radiation encounters has also proved to be of a temperamental lot, shifting 
from cooldown to cooldown. If it does move, then the path of the radiation will 
change in the spectrometer and, hence, the angles of incidence and diffraction. It is 
therefore possible that, at a certain grating position, the measured wavelength can 
vary by almost a micron for different cooldowns. For these reasons, it is extremely 
important to calibrate. We have never noticed any shift of the calibration while 
the dewar was cold for extended periods of time, and so the spectrometer is highly 
stable while kept cold. The signal levels that we measure on different cooldowns 
may vary, but usually by not more than 10%. Diurnal and seasonal changes in 
humidity may contribute to these variations. 

2.8.2 Spectrometer NEP 

A measure of how well the system is performing overall is the noise equivalent 
power, or NEP. The NEP is defined as the amount of radiation which must be 
incident upon the system in order to provide a S/N ratio of 1 in one second for a 
bandwidth of 1 Hz. The NEP can be calculated quite simply by noting the amount 
of radiation incident upon the system at a particular wavelength and for a particular 
beam, and the S/N ratio that is measured at that wavelength. If A is 153 and 
the resolution is 210, then we have 1.1 x 10 -10 W falling onto an f/17 system. 
With a S/N ratio of 2300 for a 1Hz bandwidth, this corresponds to an NEP of 
5 x 10 14 W Hz -1 / 2 . 
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2.9 The Interferometer 

2.9.1 Goals and Inspirations 

Moderate resolution spectroscopy possible with the grating spectrometer has 
made significant contributions to our understanding of the Galaxy and star forma- 
tion regions within. There are, however, some disadvantages to moderate resolution 
[R ~ 200 — 500) instruments. For the moment, consider detecting an emission line 
at 140 pm. Let us look at the Orion Nebula in the direction of the BN-KL region. 
Furthermore, let the resolution of the instrument be 210, the line width be 30 km 
s -1 , and the line flux be 1.0 x 10 -17 W cm -2 . 

The continuum radiation from the dust grains in the nebula acts as a 70 K 
grey body with an emissivity of 20/ A, A being in microns (Werner et al. 1976). With 
an instrumental resolution of 210, the flux from the continuum is 5.3 x 10' 16 W 
cm" 2 . In order to see the line, we need to be able to measure the signal levels to 
within one part in 250 for a S/N ratio of 5 in the line. This is not always easy to 
do, especially when one has shifting amounts of water vapor in the boresight. 

If we could suppress the continuum, the line strength won’t change, but the 
line to continuum ratio (LCR) will improve. At a resolution of 1000, the LCR is 
now 0.1. At a resolution of 4000, this value goes up to about 0.4. The line width, 
being one part in 10 4 , will still be situated well within the instrument profile and, 
therefore, will not be suppressed. With a LCR of 0.4, we need to measure a 40% 
change in the continuum level, as opposed to looking for a 2% change for a resolution 
of 200. Clearly, the chance of seeing the line is much better at higher resolution. 

In low continuum sources, the line-to-continuum ratio is not much of a prob- 
lem. However, another problem that exists for all sources is the fact that several 
lines may be close together. In particular, if we consider the CO J = 16 — ► 15 
rotational transition, and the OH 2 n 1 / 2 J = 3/2 + — ► 1/2“ and J = 3/2“ — ► l/2 + 
transitions, then we find that these lines occur at 162.812, 163.124, and 163.397 ^tm, 
respectively (Brown et al. 1982; Todd et al. 1976). Moderate resolution instruments 
cannot resolve these lines, and so if they are detected, confusion precludes one from 
determining which molecular line has been seen, as well as the flux levels of each 
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line. Only when the instrumental resolution is greater than 1200 can we resolve 
these lines. 

Another reason why high-resolution instruments are desirable is that if the 
line profile can be resolved, then much information can be obtained concerning 
the kinematics of the gas. Profile information can, for example, provide data on 
turbulent line widths or thermal line widths of sources, and can even allow us to 
distinguish several sources that may be in the line of sight. 

2.9.2 Design 

It was with these considerations in mind that the FIR interferometer was 
designed and built by the Harwit group at Cornell (Harwit et al. 1981). In figure 
2.8, a schematic of the interferometer is presented. Light from the telescope enters 
an aperture in the interferometer wall. The converging beam will hit a dichroic- 
coated glass which reflects the infrared and transmits the optical radiation. The 
optical radiation then goes to an offset guider which enables us to guide on a star 
as far away as 15' from the source. The guider was designed and constructed by J. 
R. Houck and his group at Cornell, and is kindly lent to us for our flights. 

From the dichroic, the FIR radiation beam converges to a focus at the 4 mm 
size aperture at mirror A. The beam then diverges upon reflection and heads toward 
mirror B. A and B are flat front-silvered mirrors. The radiation then goes to mirror 
C, where it is collimated into a parallel beam which is reflected up to D. By this 
time, the image is about an inch in size. The radiation is then directed toward 
mirrors E and F, where it is split between the two. Mirror F can move relative to 
mirror E. After this, the beam is recombined, being reflected back to mirror D and 
then to mirror G, which takes the collimated beam and converts it back to an f/17 
system while sending the light to the spectrometer. The system acts as a lamellar 
grating instrument, with mirrors E and F serving as the lamellae. 

Stacey (1985) describes the original interferometer. There were some design 
changes that were implemented in 1984 January. Mirrors E and F were replaced 
with corner reflectors. Originally, E and F were side by side, and the beam directed 
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onto these mirrors was incident at an angle. This caused walkoff problems and 
limited the travel of the movable mirror to 2 inches. With the addition of the corner 
reflectors, this problem is now eliminated. The reflectors are optically polished solid 
aluminum, overcoated with SiO to prevent tarnishing. 

Another improvement that was made to the interferometer was the increase 
in the length of travel for corner mirror F. As we will see, the resolution of the 
instrument depends on the optical phase delay between the two mirrors. The greater 
the phase delay, the higher the instrumental resolution will be. Originally, the 
distance travelled was 5 cm, but with the addition of a longer drive, that distance 
is now 16 cm. The drive stage consists of a stepper motor which turns a ball screw. 
The ball screw is coupled to the mirror mount through a mechanism which decouples 
any motions due to the ball screw that are not in a forward or backward direction. 
This insures that no wobbles due to a warped ball screw are translated to the 
movable mirror. The mirror mount, in turn, moves on a straight bar of stainless 
steel. Most of the design work was done by Noel Kurtz. 

The stepper motor is geared to the ball screw so that one step of the mo- 
tor corresponds to a 1.6 fxm movement of mirror F. The ball screw is guaranteed 
accurate to 1 part in 10000 along the length of travel, and we have found this to 
be so. The actual step size seems to be more like 1.6001 n m from measurements of 
several CO rotational lines in M42, well within the tolerances and uncertainties of 
the measurements. No detectable wobbling of the image from the movable mirror 
is observed throughout the entire length of travel. 

In order for this system to work, half of the beam must be imaged onto 
mirror E and the other half on mirror F. If the gap between the two mirrors is wide, 
then one will begin to lose signal and interference. We have found that most of 
the interference is due to the central part of the beam and it is therefore necessary 
that the mirrors be placed as close together as possible without touching. We have 
managed to place the mirrors ~ 0.004" apart, or the thickness of a piece of paper. 

The interferometer is constructed out of aluminum, while the ball screw is 
stainless steel. The stepper motor is a source of heat, and it can raise the tempera- 
ture of the ball screw by 5 or 10 degrees Celsius. Since the ball screw is supported 
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by ball bearings at both ends, it is not in good thermal contact with the rest of the 
interferometer, only with the stepper motor through the gear train. The large ther- 
mal mass and effective cooling area of the interferometer prevents it from heating 
up very much, but the ball screw has little thermal inertia and poor thermal contact 
with the rest of the instrument. Hence, the ball screw will heat and expand. In 
so doing, it effectively has a temperature only about 5 degrees Celcius higher than 
the interferometer. The ball screw is not allowed to physically move back and forth 
for obvious reasons and is thus prevented from doing so at the stepper motor end. 
Hence, when the screw expands, it results in moving corner reflector F toward the 
dewar by about 16 /xm. Thermal equilibrium is achieved after the stepper motor 
has been on for an hour or so. Until then, however, the position of the movable 
mirror is changing. The resultant change in the step size due to the expansion of 
the screw is only 1 part in 10000. 

2.9.3 Theory 

The radiation from the telescope first encounters the interferometer and then 
is passed to the grating spectrometer before ultimately hitting the detector. Let us, 
for the moment, assume that 5(A) is the spectral profile of the radiation coming 
from the telescope. The grating spectrometer instrument profile is given by 'I'(A). 

For radiation of wavenumber u, incident on the corner reflectors, we have 
half of the beam falling on one reflector and half on the other. For an optical path 
difference between the corner reflectors, x , the value of the electric field vector upon 
recombination of the two beams is 


£(x) = Y + Y (2*11) 

where Eq is the amplitude of the electric field for the entire beam. The intensity of 
the radiation is simply EE* , or 

I(x) = ^[1 + cos(2 nux)\ , (2-12) 

St 


where Iq is Eq. What the detector will see is this intensity convolved with the in- 
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strument profile and the spectral profile. If <f>(x) is the intensity hitting the detector, 
then it follows that 

4>(x) = 77 / dt/ 't(iz) S{u) [1 + cos(27ri/z)] . (2.13) 

2 Jo 

The first term on the right hand side is a constant d.c. value we call Aq/ 2. 
Ao is the average of 4>(x) over all x. The second term is the one which contains all 
of the information about the incident radiation. Since the argument in the integral 
is symmetric, the integral may be rewritten as 

= £ + T r iv *(») S M • (214) 

2 4 y_oo 

If we let $(x) = <j>(x) - Ao/2 and T(i/) = Jo 5(i/) \&(i/), then (2.14) may be rewritten 
as 

$(x) = - [°° dvl{v) e- 12 ™* . (2.15) 

4 J-oo 

We recognize the integral to be the Fourier integral of T (u), or in other words, the 
signal at the detector is 1/4 the Fourier integral of T(i/). Conversely, we can write 

r(i/) =4 ?(j/) , (2.16) 

where $(t/) is the Fourier transform of <£(x). 

If we know the instrumental profile, '^(t'), then by calculating T(t/) from 
$(x), we can determine the spectral profile, S(u). This profile is actually a convo- 
lution of the atmospheric transmission profile and the source profile. Usually, the 
transmission profile is relatively flat over the spectral range of interest, since that 
range is small. Therefore, we can assume that S(u) is proportional to the source 
profile. 

In practice, it is impossible to perform Fourier integrals since that implies 
talcing data over an infinite and continuous range of optical path differences. We 
can sample the spectrum by moving the mirror in bursts of A/2 (for an optical path 
difference of A), and do this for 2n points centered Mound zero path difference. The 
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resulting interferogram will consist of a discrete series of data points, $(x r ). We 

can represent this function as a Fourier series, 

/ 

n— 1 

$(x) = Ao + 2 y: {A m cos(27Tn wox) + B m sin(27nm/Qx)} +A n cos^Tmt/px) . (2.17) 

m=l 

Here, x = r A, n < r < n — 1, and vo is the fundamental wavenumber equal to the 
inverse of the total optical path difference between the first and the last mirror 
positions, i.e., i/q = l/2nA. It should be emphasized that 4(x) is only an estimate 
of $ (x). 

By using the orthogonality properties of the cosine and sine functions, we 
can easily determine the coefficients: 

Am = ^ Y1 *(**) cos (^r) » ( 2 - 18 ) 

r=— n 

Bm = ^ H $ ( Xr ) sin (^^) » ( 2 - 19 ) 

r=— n 

for m = 0, l,...,n. One consequence of a finite Fourier series is that the spec- 
trum covers a finite range. We see that the range of wavenumbers covered is 
v = 0,t/o,2t'o> • • • >ni/(). This range, which is equal to 1/2A, is the free spectral 
range, or FSR. Owing to the periodic nature of the sine and cosine functions, a 
mirror image of the spectrum will exist for 1/2 A < u < 1/A. This is also the case 
for the order on the other side. The spectrum will then alternate mirror images of 
itself in wavenumber space. 

It is necessary to alias the orders so that we are only observing radiation 
coming from one order. That is to say, the radiation falling onto the detector 
cannot have a bandwidth Su > 1/2A. If 6u is larger than the FSR, then the 
radiation which falls outside of the FSR will be folded back onto the spectrum. 
Aliasing is accomplished by using the grating spectrometer. For more information 
on Fourier series, integrals, and spectroscopy, see Gaskill (1978), Bell (1972), and 
Jenkins and Watts (1968). 



One thought that is usually neglected is the fact that the radiation incident 
upon the two lamallae is not necessarily in phase with itself. If this is the case, 
how does this affect the interference patterns? Consider a beam of radiation that 
contains a random phase distribution. If half of the radiation is incident upon each 
of the lamellae, then each corner reflector should be illuminated by a subset of the 
photons with random phase distributions. Let us further assume that there are 2N 
wave packets in the beam, each packet having a phase <t>}. Let the optical path 
phase difference between the mirrors for radiation of wavenumber v be e. Then we 
find that the total electric field upon recombination is 

N N 

Em « £e-*J + ^e-^i +< ) . (2.20) 

j=l j=l 

To find the intensity, we find EE*. We find that 
N N 

/ a J2 J2 {2e- i( *J-*k) + + e _i Wr* +<) } . (2.21) 

j=l k=l 

Owing to the randomness of the individual phases, they are not correlated with 
each other and so the only non-zero terms that we find are for j = k. This reduces 
(2.21) to 

/« 2(1 + cost). (2.22) 

Throwing in the appropriate constants of proportionality, we find that this is equiv- 
alent to (2.12), which was the result for a single wave-packet. 

2.9.4 Determining Burst Sizes 

The bandpass of the grating spectrometer is « 0.7 /xm FWHM. If we wish 
to effectively alias the orders we obtain from the interferometric action, we need 
to consider two things. First, the grating bandpass should be centered within the 
order. If A is the wavelength corresponding to the peak of the grating bandpass, 
then it follows that 


1 


1 


1 
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where m is the spectral order. Secondly, the grating bandpass should fit within the 
FSR. 

If we choose A and m, then by (2.23), we know the optical burst size, A. The 
mirror will move only half that amount, or 

_ M m + ?) 

^mirror “ ^ 

Knowing the mirror burst and the distance travelled per step of the stepping motor, 
we can determine the number of steps to drive the stepper motor. 

It is apparent that the FSR is determined solely by A. For a fixed A, the 
resolution, then, is determined by the number of data points that are taken. If a 
32 point interferogram is taken, then the instrumental resolution will be the FSR 
divided by 16. With a 64 point interferogram, the resolution will be finer by a 
factor of 2 and so on. Hence, it is the total distance travelled that determines the 
resolution. By proper choice of A and m for a particular wavelength A, as well as 
the number of data points, one can properly alias the spectral order and realize 
high resolution. 

The actual operation of the interferometer is accomplished by taking one- 
sided interferograms and then reflecting them about zero-path. In order to do 
this effectively, zero-path must be known and maintained. Once the instrument is 
allowed to come to thermal equilibrium, the drift of zero-path throughout a flight 
may be no more than 6 nm, which corresponds to errors of ~3% in the determination 
of the spectrum. Reflection of the spectra can be done since we may assume that 
mirror images would be obtained if we did sample on both side of zero. Reflecting 
saves half of the integration time. 

2.9.5 The Fellgett Advantage 

It would seem to be much more straightforward, in the quest for high spec- 
tral resolution, to simply build a bigger grating instrument. This has been done 
by Erickson et al. (1985). Aside from the difficulties of making a large grating (at 
150 fxm, the grating must be at least 50 cm long for a resolution of 5000) and an effi- 
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cient cryogen can to keep it and the detectors cold, there is a distinct disadvantage 
if the system noise is primarily electronic and not background limited. 

Consider two systems, one a grating and the other an interferometer. Allow 
both to look at a spectral region of width A u in increments of Su. Let the total 
integration time be T. Then, for the grating, M steps must be taken, where M 
= AufSu. The total integration time per band Su is simply T/M. If the noise is 
independent of signal strength, then the noise of each spectral band is given as 

n„oc(T/M)*. (2.25) 


For the interferometer, 2M data points must be taken in order to obtain 
comparable resolution. However, since each band is encoded in each point of the 
interferogram, the effective integration time per Su is T. The noise per band is 
then proportional to T 1 / 2 . The S/N ratios for gratings and interferometers may be 
written as 


(S/N /Ms* 

(S/N) G[ „ V 8 J ’ 


(2.26) 


This is known as the Fellgett advantage. The factor of 8 comes into play because of 
the fact that the signal is modulated by a cosine with an amplitude equal to 1/2. 
The root mean square value for this term is simply l/%/8. 

If, on the other hand, the noise is proportional to the signal strength, the 
Fellgett advantage is lost. That can be seen by the following argument. For a 
background limited system, let I{Su) be the intensity of radiation in the band Su. 
Then, other things being the same as above, for the grating, the noise should be 
proportional to the square root of the integration time per point and the band 
intensity, or 


(S/N) g ,„«[(T/M)/(H]5 • (2-27) 


For the interferometer, the signal in a band is proportional to TJ(6i/), assuming 
that I{Su) is constant over the spectral range. At any one mirror position, all of the 
bands will contribute to the noise. In performing a Fourier transform, the noise will 
stay about the same, and be proportional to TM Hence, the signal to noise 
ratio for the interferometer is also proportional to [(T/M)/(<5z/)]5. That is to say 
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that the Fellgett advantage is lost once the photon statistics dominate the source 
of noise. For more information on the Fellgett advantage, see Harwit and Sloane 
(1979), Bell (1972), and Carli (1972). 

2.10 Operation and Performance 

2.10.1 The Setup 

In section 2.9.4, we discussed the choosing of the proper burst size and res- 
olution for an interferogram. One of the constraints is that we properly alias the 
spectral orders. This is done using our grating spectrometer, which has a spectral 
bandpass of roughly 0.7 ^m FWHM. We have found that a FSR of 2 (im or more 
allows us to properly alias with little or no folding over of the grating profile. 

For flight operations, the interferometer is mounted onto the backplate of the 
telescope at the Nets myth port. The mounting ring is fitted with o-rings on both 
sides, as are the top and bottom lids of the interferometer, to allow a good seal. 
This is important since the telescope cavity is at ambient pressure during flight at 
41000 feet altitude while the plane’s interior is kept at a pressure equivalent to an 
altitude of 8000 feet. If the seal was not tight, then moisture laden air from the 
plane interior would be travelling through the optical path, creating an additional 
source of noise. 

The interferometer is also electrically isolated from the telescope by a thick 
ring of polyethylene and plastic sheaths and washers for the mounting bolts. The 
telescope focus is adjusted so that the focus of its f/ 17 beam is at the 4 mm aperture 
at mirror A in the interferometer. The dewar is then attached to the interferometer, 
as is the offset guider. 

A typical observing program dictates that several emission lines of interest be 
observed during any one flight. In order to properly do this, the grating is calibrated 
relative to its stepper motor so that it may be positioned for the wavelengths of 
interest. The burst sizes for the interferometer are calculated for each wavelength, 
usually for several resolutions and free spectral ranges. Zero-path must also be 
found for the interferometer. 
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A data acquisition program written by Charles Fuller for the IMS AI computer 
is used to take the data and drive the stepper motors. Once the FIR source is 
obtained and centered in one of the chopped beams, the data are taken by setting 
up the grating and then integrating for 5 seconds or more at each mirror position 
of the interferometer. Once one interferogram is acquired, the data are written to 
floppy disk and also sent over to the ADAMS computer onboard. The interferometer 
is stepped backwards for the next run. This saves time in driving the mirror all 
the way back to zero path (upwards of 90 seconds for the full travel) and beginning 
again. The IMSAI is dedicated during this time to the data acquisition. 

With data in hand on the ADAMS system, the Fourier transforms can be 
performed at an operator’s terminal linked to the ADAMS computer. The program 
performs a fast Fourier transform (FFT) and was written by Scott Smyers with 
some modifications later added by various members of our group. The FFT routine 
is an algorithm that can be applied to an interferogram with 2 n points, and takes 
only a fraction of the computational ability that a normal discrete Fourier transform 
would tak The algorithm is an adaptation of that published by Cooley and Tukey 
(1965). 

Rapid data reduction permits flexibility to change observing plans during 
flight. Thf FFT program uses the raw data, and so gives us a relatively good idea 
of what we are actually seeing. Post-flight data reduction includes searching the 
individual integrations for glitches due to cosmic rays or radio frequency interference 
and redoing the Fourier-transforms again. Various lag windows may be used at this 
point to effect an apodization. These are discussed in Stacey (1985). 

2.10.2 Fringe Visibility and the Spectrometer 

According to theory, the signal level for a monochromatic beam should fall 
to zero at the places of maximum destructive interference. In practice, the signal 
level doesn’t. There are several reasons for this. The first and most obvious one 
is that the beam is not monochromatic. The grating resolution of 210 limits the 
amount of total destructive interference that one can have for a certain optical path 
difference, A. However, this limit still provides a maximum-to-minimum ratio, or 
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fringe visibility Q t of ~ 1000, and can be shown to be equal to 



(2.28) 


where R is the grating resolution and m is the spectral order for the interferogram. 

Misalignment of the interferometer could allow one of the corner reflectors to 
have more than half of the radiation falling on it and the other, of course, to have 
less. This would provide an imbalance which would act to lower the modulation. If 
/? is the fraction of the radiation falling onto one of the mirrors, then Q is 


Q = 


l 

1 - 4/3 + 4(3 2 ’ 


(2.29) 


At 0 = 0.5, Q is infinite. At /? = 0.6, Q = 25. For 66% of the radiation falling onto 
one mirror, we have Q = 10. 

In the above instrument configuration, we have found Q ~ 4 at best. It is 
possible to measure the radiation coming from each mirror, and this suggests that 
this is not the limiting factor. We have seen above that the resolution is also not 
to blame for such a low Q. Another factor which we have not considered is the 
spectrometer itself. If we remember, the central 33% of the beam is blocked due to 
the Littrow-mode operation. It is most likely that up to one-half of the interference 
is contained within this part that is discarded. While using the hinge Fabry-Perot 
spectrometer (outlined in chapter III), we have found Q on the order of 12. This 
instrument is an off-axis system and, therefore, has no central obstruction of the 
beam. For this reason, we believe that the efficiency of the interferometer is limited 
mainly by the grating spectrometer. Considering (2.12), we can see that a low Q 
implies a smaller modulation. This, in turn, implies that the Fourier transform will 
give results that are smaller by some efficiency factor. It is easy to show that the 
efficiency of the interferogram is 1 — Q~ l . For a Q of 4, the efficiency is 75%. For 
a Q of 12, that number climbs to 92%. The efficiency is a slowly varying function 
of Q once Q is of order 4 or more. 



III. The Hinged Grating Fabry-Perot 
FIR Spectrometer 

3.1 Introduction 

In Fourier transform systems (FTS), the resolution is dependent upon the 
total distance of travel of a mirror. At 150 pm, a resolution of 5000 requires a 
total travel of 38 cm. In order to resolve fine detail in molecular line emission, it 
is often necessary to have resolutions of 25000 or more. For a FTS, this means 
that the mirror must be driven a total of at least 1.88 meters. Constructing such 
an instrument with such a long travel length is difficult at best. Not only must 
one obtain a driving mechanism which is accurate to several microns over the total 
distance, but one also needs a guide that is straight over this distance for the 
movable mirror. To guard against flexure, the instrument would have to be severely 
massive. Our present interferometer has a weight of roughly 30 pounds, which 
would imply that a FTS with a resolution of 25000 would have a weight of over 150 
pounds. An instrument of this size and bulk could not possibly fit on either the 
KAO or the Lear Jet. These factors alone would also lead to pointing and tracking 
problems for balloon-borne systems, and would preclude the use in space-borne 
infrared packages. 

Gratings offer different problems that also prohibit using them for high res- 
olution in the FIR. Again, if we desire resolutions of 25000 at 150 ^tm for a grating 
with 10 rulings per mm, we will need a grating that is at least 2.5 m in size. It is 
desirable to have the grating cooled along with the other optics and filters in order 
to minimize the background radiation. Building and cooling a spectrometer with a 
grating of this size is not an easy or practical proposition. 

Obviously, if the optical path can be folded many times, we can have a high 
resolution system which fits in a small space. Fabry-Perot spectrometers (FPS) are 
just such a system. It is possible to easily obtain resolutions of 25000 in the FIR 
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with a Fabry- Perot that is no more than 10 cm in size. 

Since a Fabry-Perot is operated much like a grating instrument, there is no 
multiplexing, or Fellgett, advantage. However, for the FIR, detector sensitivity is 
at present photon noise limited and we know that in this situation, the Fellgett ad- 
vantage is lost by a Fourier transform system, such as an interferometer. Therefore, 
a FTS cannot offer an advantage in the signal-to-noise ratio over a FPS. 

For these reasons, we have been inspired to design and construct a Fabry- 
Perot spectrometer system. This system, as will be shown in this chapter, is based 
upon placing a grating within the FP system itself. The etalon plates are fixed in 
position, and it is through a rotation of the grating that the FP is timed. Because 
of this rotating action, the instrument is known as ‘The Hinge’. 

We shall first examine ordinary FPS systems in order to understand their 
characteristics, and then shall present the analysis, both qualitative and quantita- 
tive, of the Hinge system. The last part of this chapter will deal with laboratory 
results that demonstrate the workability of this unique Fabry-Perot instrument. 

3.2 Fabry-Perot Spectrometers 

3.2.1 Theory 

FPS systems have been employed for many years, mainly in the optical and 
near infrared spectral regimes. For an overview of these FPS applications, see 
Roesler (1974) and Vaughan (1967). Basic FP theory is straightforward. We will for 
the most part adopt the nomenclature of Roesler (1974). Figure 3.1 demonstrates 
the Fabry-Perot action. Light from a collimated beam enters the Fabry-Perot etalon 
from the left. The etalon consists of two semi-transparent, optically flat windows 
that are parallel to each other and separated by a distance o. After the light enters 
the cavity, it begins to bounce back and forth between the two plates. With each 
bounce, some of the radiation is passed through one plate and the rest is reflected 
to the other plate. Collecting optics on the right side of the etalon will focus the 
rays at the detector. 

Let r and t be the reflection and transmission coefficients for the electric field 



55 



Yabry 


C oUVta» ted ^d^ ttU 


L VttW°' 


•^s^JSsSbrs*- 


l- *££*■ 

"'** nt «\e®^ »' 1 1 »« eUOi 'Sp 
\$X Ctfe*® 



56 


vector with wavenumber v. If Eo is the amplitude of the initial electric field vector, 
then tEo will pass to the inside of the cavity. Upon encounter with the right-hand 
plate, some of the light will be passed (t 2 Eq) and the rest (rtEo) will be reflected. 
Of this reflected light, a fraction rt of it will pass through to the collecting optics. 
This light will be phase shifted from the initial beam by <t>, where 


<f> = 4ir<ri/sec0 


(3.1) 


The variables refer to those marked in figure 3.1 and we are assuming that the index 
of refraction is unity. As we follow the individual reflections and transmissions, we 
see that each subsequent transmission through the etalon is multiplied by a factor 
of r 2 e -1 ^. Summing over all of the transmissions gives us a total electric field of 


oo 


E to t =£oXy r2n e _in * . 

n=0 


(3.2) 


Noticing that the sum is simply a geometric progression and solving for the intensity, 
we find that 

rp2 

(3.3) 


I{4>) h ^ j £ cos 4> + R 2 

where T is the transmission and R is the reflectance of one plate. Iq is equal to 
EqEq. Algebraic manipulation will turn (3.3) into the well known Airy equation, 


I{4>) = Iq 


j*2 


(1 - *)* 


1+ (TTfl)I s ' n <2> 


-1 


(3.4) 


If there is no absorption (J2 + T = 1), then (3.4) reduces to 


m = la 


4 R , 

1 + (T^Bp sm ( 2’ 


1-1 


(3.5) 


It is easy to calculate the reflected intensity from the etalon, 


f,.a(d) = /o 


4 R sin 2 (£) 


(1 - R) 2 + 4 R sin(^) 


(3.6) 



Equations (3.5) and (3.6) show that energy is conserved in the process if there are 
no absorption losses. 

When the sine term is zero, the transmitted intensity is 100% of the incident 
intensity. This occurs whenever 

2o v seed = m; m = 0,1,2, ... . (3.7) 

The transmitted radiation consists of a series of fairly narrow peaks that are uni- 
formly spaced in wavenumber (see figure 3.2). The spacing is called the free spectral 
range and is denoted by Q. From (3.1) and (3.7), we see that 

Q = cm 1 . (3.8) 

The instrumental resolution is not Q , but is the FWHM of one of the trans- 
mission peaks. One can easily derive that this St/ is given by 



For R ^ 0.8, (3.9) is approximated as 

e Ql-R /0 , ft . 

Su = — — t=-.. (3.10) 

* s/R ' 

The reflective finesse is the ratio of the free spectral range to the FWHM and is 
given by _ 

tfR = . (3.n) 

from which we find the theoretical resolution to be 

= ^ = mN * ’ (3 - 12) 

where m is the order number given by (3.7). The contrast, C, is the ratio of the 
maximum transmission to the minimum transmission, and is given by 
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FIGURE 3.2 A transmission profile for a Fabry-Perot etalon. 
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If continuum emission or closely spaced lines are present, it is important to have a 
high contrast in order to avoid confusion in line flux measurements. The contrast 
is a function only of the reflectivity. 

3.2.2 Limitations to the Finesse 


The above theoretical description of a FPS is for a perfect system. The 
usual measure of performance for a FP is its finesse. In a perfect system, the finesse 
is given by Wr. No system is perfect, however, and the limitations of various 
components will result in a lower resolving power and hence, in a lower value for 
the instrumental finesse. 

There are several factors that limit the finesse of a FPS. The first is the size of 
the entrance aperture. Since the aperture is extended, radiation incident upon the 
etalon is not totally collimated. Hence, if /? is the angle subtended by the aperture 
as seen by the etalon, then it follows from (3.7) that the range of wavenumbers that 
are passed iu a given order is 

= "Omt! ~ c °s(^)] , (3.14) 

where uo m b the wavenumber of the radiation that is passed through in order m for 
a perfect system. We assume here that the radiation is incident normal to the first 
etalon surface. Since the apertures are small, we can replace the cosine function 
with the first term in the Taylor expansion to obtain the familiar 

(3.15) 

The aperture finesse b then given by 


W _ P 2 . 

i'Om 8 


_ Q _ 4 _ 8Q 

Su A ~a/3 2 u 0m 0 2 vo m ' 


(3.16) 


Diffraction effects will also play a role in degrading the resolution. Consider 
an aperture of size a. Then we know from diffraction theory tha: the diffraction 
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pattern as a function of angle B is 


1(B) = 4/q sine 2 (tt i/o sin0) 


(3.17) 


At half-maximum intensity, the argument of the sine function is 2.24. This leads to 
an angular spread in the beam of 


^diff = 


1.43 


v a 


(3.18) 


Again, from (3.7) we find that the diffraction finesse is given by 

2a 2 i/ 

tfdiff = — - • 


(3.19) 


Another component of a FPS which affects the finesse is the flatness of the 
plates comprising the etalon. The finesse due to nonflatness can be calculated in 
the following manner. We divide the etalon into a collection of many microetalons. 
Each microetalon is characterized by a separation, 0 $. If the a x constitute a normal 
distribution about a mean, oq, with an error 6a, then it follows from (3.7) that 


6u 

u 


6a 

a 


(3.20) 


and that the nonflatness finesse is 


(2 - 21) 

The shape of the bandpass will be affected by the distribution of the nonflatness 
errors. Here, we simply assume that it is a normal distribution. 

If the plates in the etalon are not parallel to each other, then another degra- 
dation of the performance will occur. This tilt finesse can be obtained in the same 
manner as the nonflatness finesse above. Let one plate be tilted with respect to the 
second plate by an angle, Bt- Then it follows that the finesse due to the tilt is 


(3.22) 
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where D is the diameter of the beam as it enters the etalon. 

The FWHM contributions from each of these factors add in quadrature. This 
implies that the resulting instrumental finesse is given by 

J_ _ ^ i _L _j_ 1 _| L _|_ _L_ (3 23 ) 

N 2 Nl Nl Nis Nf ' 

We therefore find that the instrumental finesse can be no better than the worst 
finesse imposed by one of the above mentioned considerations. 

3.2.3 The Effect of Absorption 

One certainly will expect some losses from the plates due to absorption. For 
each plate, we let T be the transmittance and R be the reflectance. If we denote 
the absorption by A, then it follows that 

T + R + A=l. (3.24) 


Using (3.4) and substituting for A from (3.24), we find that 

m = * (i + t^r) [i + 


(3.25) 


It is important to note that the absorption will not change the characteristics of the 
system such as finesse or resolving power. It will, however, limit the transmission 
in a non-linear way. When A has the value of 0.3(1 — R ), or T = 0.7(1 — R), the 
transmission of a peak is down to 50%. Even if the absorption is one or two percent, 
this is an important factor if high resolution, and therefore, high reflectivity, plates 
are desired. 


3.2.4 Isolation of a Spectral Element 

As we saw in section 3.2.2, the FP will transmit a series of relatively narrow 
lines separated by the free spectral range. Changing the spacing of the plates will 
change the wavelengths of the transmitted lines as dictated by (3.7). In actual 
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usage, only one of these lines is to be isolated. The suppression of the other lines 
must be accomplished by some method, and several such methods shall be listed 
below. 

The simplest method to isolate a line is to use a combination of filters. It is 
difficult to obtain a narrow bandpass with a filter stack, and one should not expect 
to have a bandpass narrower than a few tens of microns. If one is using a high 
resolution FPS, then the finesse will have to be several thousand. This is not a 
practical finesse for the FIR, where finesse values on the order of a few tens to one 
hundred are more appropriate. The diffraction and aperture finesse values are not 
very large, and to make them so requires having a large entrance aperture located 
far from the etalon. 

Another problem with using filters to suppress the unwanted orders is that 
you are very limited in the range of wavelengths that you can observe. This is 
undesirable since not all FIR lines of astrophysical interest are contained within a 
narrow band of the spectrum. 

The second method of order sorting that is used is a combination of filters and 
a grating. The filters can cut out most of the radiation below a certain wavelength, 
and then the grating can be used to select a bandpass in the first order of the 
grating. Again, if one desires high resolution and the instrumental finesse is a few 
tens, then the grating resolution must be about one thousand or so. In order to 
achieve this resolution in the FIR, one will need a grating at least 10 cm in size. 

A third method that can be used is to combine filters and another FP that 
works in a different order in tandem with the primary FP. In this way, it is possible 
to effectively sort out one order over a rather broad wavelength range with the 
second etalon and then let the filters suppress any sidelobes further out. This can 
best be accomplished when the second FP is used in a very low order. Figure 
3.3 demonstrates this method. The primary FP in the system is set for a specific 
wavelength in a certain order and finesse. The second FP is set also for that 
wavelength, but in first order. The resulting instrumental profile is a convolution 
of the two profiles from the two FP etalons. This is shown in the bottom box of 
figure 3.3. 
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Wavelength (p m) 


F IGURE 3 . 3 Transmission profiles for a Fabry-Perot working in a high order (top) , 
for one working in first order [middle), and for a tandem FP (bottom), incorporating 
both of the above Fabry-Perot systems. 
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One notices right away that the suppression ability of this setup is very good. 
The central order is changed very little due to the wide envelope of the second FP. 
Since we are in first order, the second order wavelength is half of the first and for 
the FIR, the higher order radiation is well blocked by salt filters. 

For more information on Fabry-Perot systems, see Watson (1982), Born and 
Wolf (1980), Roesler (1974), and James and Sternberg (1969). 

3.3 The Hinge — A Qualitative Analysis 

3.3.1 Inspiration 

In a normal FPS, the spacing between the plates is adjusted using piezo- 
electric stacks. The spacing is then a function of the voltage applied across these 
stacks. Owing to the nature of such a design, the electronic circuitry can become 
very complicated. Maintaining parallelism of the plates ofttimes requires human 
intervention, a service that is missing if such a device is sent into space aboard a 
satellite. This is not to say that a FPS cannot be successfully launched and op- 
erated, only that such an instrument must be flown with a retinue of complicated 
circuitry and instrumentation in order to maintain alignment and parallelism, as 
well as determining the spacing between the plates. 

A tandem FP instrument is limited in the range of wavelengths it can cover. 
Usually, only a few narrow spectral regions can be seen, and the lack of continuous 
spectral coverage is a severe handicap. Because of this, the FP/grating combination 
is considered more desirable. However, the complications that are associated with 
a variable etalon are also found in the FP/grating instrument. 

If it were possible to vary the spacing between the plates while keeping the 
plates stationary, this would be quite helpful. For various reasons, simply tilting the 
etalon to vary 0 is not always am acceptable solution (Roesler 1974). If a reflecting 
surface could be placed within the etalon cavity, then movement of the surface could 
tune the cavity. However, there is still the requirement that the reflecting surface 
must move in some predetermined way. The easiest way to move a surface is to 
rotate it about an axis. If that surface were the grating, then it might be possible 
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to tune the FP and the grating at the same time. 

Not only can such an instrument be compact, but if there is only one moving 
part and that part rotates, then it is also extremely sturdy. Another advantage is 
that it should be tunable over a broad range of wavelengths. Such an instrument 
would be perfectly suited for high-resolution spectroscopy in FIR space missions. 

3.3.2 Principle of Operation 

In order to realize such a design, we should first try to understand a diffrac- 
tion grating. Consider the setup depicted in figure 3.4. A diffraction grating works 
on the principle of interference. Each ruling in the grating will act as a slit. If a 
collimated beam of radiation is incident upon the grating at an angle of incidence 
a and is diffracted at an angle /?, then by the grating equation (2.8), we know the 
wavelengths that will be passed through. The path difference between rulings turns 
out to be an integral number of wavelengths. In first order, the difference is one 
wavelength; two wavelengths for second order, and so forth. We can consider a 
fiducial ruling on the grating and construct from that ruling two planes that inter- 
sect the ruling and are perpendicular to either the incident beam or the diffracted 
beam (see figure 3.4). It follows that radiation perpendicular to both planes will 
add constructively only for the wavelengths that satisfy the grating equation. 

There was nothing special about the angles a and (3. Since they can change, 
we cam consider what will happen when we rotate the grating. The only change that 
occurs is that the wavelength of the radiation that is passed through the system will 
change. It is still true that radiation perpendicular to both plates will constructively 
interfere for the wavelengths that satisfy (2.8). 

If we place highly reflecting plates where these planes are, then we find that 
we have a tunable FP cavity. The radiation that enters the etalon will be reflected/ 
diffracted back and forth between the plates and the grating. 

The instrumental finesse of the system is a measure of the number of times 
that the light bounces inside the cavity. For a finesse of 10, that means that there 
will be roughly ten bounces. On each bounce, the light will hit the grating twice, 
and so in effect the beam will appear to have hit 20 times the number of rulings that 



FIGURE 3.4 Principle of operation of the Hinge. Light enters along the direction 
AC, and is diffracted by the grating along CE. For the two beams presented here 
(thin solid lines), the phase difference between them is represented by BC + CD, 
which is equal to one wavelength. The total optical path, AC + CE is equal to an 
integral number of wavelengths, since F is coincident with a grating step. If highly 
reflective meshes are placed along the planes represented by AF and EF, then one 
has a FP system convolved with a grating in the middle. 
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it illuminates. This means that if the grating has a normal resolution of about 200, 
then it will appear to have a resolution of roughly 4000. Hence, not only will the 
grating tune the FP cavity, but it will act as a very good spectral element isolator. 

3.4 The Hinge — Theory 

Placing a grating inside a FP etalon alters the characteristics of the etalon. 
Since each successive ruling of the grating produces a relative phase shift of one 
Wavelength, the cavity can be thought of as being composed of a series of mi- 
croetalons, each microetalon having a spacing one wavelength longer than the pre- 
vious one. In other words, each microetalon will operate in a different order, no 
two orders being the same. Because of this, the transmission profile of the hinged 
grating FPS will differ from that of a regular FP/grating system. 

To find out what the transmission of the Hinge should be, we consider the 
following construction. In figure 3.5, we have the grating and the two plates that 
make up the FP cavity. Since it is not possible to have the plates actually intersect 
at the grating axis, the beam illuminates the grating an integral number of rulings 
down from the fiducial ruling. 

We begin by considering each of four phase delays for the different rulings. 
Every part of the beam will undergo the two phase delays and ^4. Additional 
delays due to the grating are $2 and <f> 3. If the angle of incidence upon the grating 
is a, then for radiation of wavelength A, it follows that 

2tt 

<£l = — Ni psina , (3.26) 

where g is the ruling spacing on the grating and Ni is the number of the first ruling 
to be illuminated, starting with N = 0 for the ruling coincident with the rotation 
axis. Likewise, if the diffraction angle is /?, 

04 = yN 1 ysin/?. (3.27) 

The grating will also contribute to the phase delay. The radiation that falls 
on the m th ruling from the first illuminated one will undergo an additional phase 
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delay of 

. 2tt 

</»2 = — m g sina , 

due to the incident path difference and 


, 2tt . 

91 = y m^sin/? , 


(3.28) 


(3.29) 


due to the diffraction path difference. We see that for radiation to go from plate A 
to plate B, for each ruling the phase delay will be 


= 4>\ + <f>2 + 4>z + 


<t> 4 • 


(3.30) 


Using (3.26) to (3.29), we find that 

2 7T 

9 = -7- (Ni -l- m) g (sina + sin/3) . (3.31) 

In a FP etalon, the phase difference is due to the optical path from one plate to 
the other and back again. Obviously, if </> m is that phase difference, then <f> m = 2<f>. 
We use the subscript m to indicate that the phase difference will depend on which 
ruling is illuminated, given that the other parameters are fixed. 

Since each ruling defines a microetalon, it follows from (3.2) that the total 
electric field for this microetalon for the m 1 * 1 ruling is 


E m = E 0 ^ * 2 r2n e“ i(n+ 5^ n 


n=0 


= EqT 


1 - Re 1 *™ 


1 + .R 2 — 2R cos <(> 




m 


= £ ° T 1+ J ' ~2R cos* ' 3 in» 


(3.32) 


To find the total electric field from all of the microetalons, we sum over the M 
illuminated rulings to get 


Etot E o T 2_, + r2 _ 2R cog 

m=0 


(3.33) 
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The intensity is then given by 


I = IqT 2 


( r M_1 / 

v L m=0 v 


(1 - R) cos(^f) 

+ R 2 - 2 R cos <f> m 


)f 


+ 


rM-l / 

£ I 

L m=0 N 


(1 + iZ) sin(% ) 

+ i2 2 — 2R cos 4> m 


)]'} 


(3.34) 


The calculation of I must be done numerically, since it is most likely that (3.34) 
cannot be written in an analytic form. More will be said about the computations 
later. 


3.5 The Hinge — Construction 

3.5.1 Design Considerations 

One of the advantages of the Hinge is its compactness. The Hinge was 
designed to fit within the present dewars that we have been using for the past five 
years or so. These dewars have a diameter of 6 inches and can house a spectrometer 
up to 10 inches tall and 5.5 inches in diameter. This size constraint did not prove 
to be much of a problem. 

Figure 3.6 is a schematic of the Hinge assembly. The Hinge was designed 
as an off-axis spectrometer. Light enters the Hinge and is baffled by two entrance 
apertures. Black poly, CsBr, and BaF 2 filters are placed between the entrance 
apertures to block out unwanted optical and near infrared radiation. 

The diverging f/ 17 beam is directed by two planar diagonal mirrors along 
the perimeter of the spectrometer to a collimating mirror. From the collimating 
mirror, the radiation is directed to the FP cavity. For the FIR, the etalon plates 
are usually made of a grid or mesh. 

Since we are constrained to place the intersection of the mesh planes coinci- 
dent with the rotation axis of the grating, the meshes need to be accurately aligned. 
This can most easily be accomplished by constructing a mesh/grating holder which 
accurately places the grating pivot point in the mesh planes. The meshes are then 
allowed to rest upon the holder, which serves to define their planes. 
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Entrance Aperture/ ‘Detector 

and Diagonal Mirror and Cavity 


FIGURE 3 . 6 Schematic diagram of the FIR Fabry-Perot hinged grating spectrom- 
eter. An f/ 17 beam of radiation enters at the bottom, perpendicular to the page, 
and is then reflected towards the collimating mirror at top left. The collimated 
beam is then directed into the FP etalon, where it is reflected and diffracted a 
number of times before exiting to the post-optics, which focus the radiation down 
onto the detector. 
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The grating plane must intersect the two mesh planes at the pivot axis. Since 
the beam will not actually be hitting the grating near this position, but several 
centimeters away, it is not necessary to rule the entire grating. It is necessary, 
however, that the rulings are an integral number of spacings away from the pivot 
point. 

With the present size of the FPS, the grating can be rotated through 14°. 
A somewhat larger spectrometer could easily accomodate a larger rotation. As 
it is, the angle of incidence on the grating ranges from 26-40°, with the angle of 
diffraction being 38° larger than the angle of incidence. The corresponding range 
of wavelengths is from 140 to 170 /xm for a ruling spacing of 104 /xm. 

The grating is made of stress-relieved aluminum, as are the rest of the spec- 
trometer’s structural and optical members. This allows the entire structure to 
contract uniformly upon cool down, alleviating problems with unnecessary stresses 
and shifting focal planes. The bearing housing and the axle for the grating pivot are 
stainless steel, but owing to the circular symmetry, the differential contraction upon 
cooldown will not shift anything. The grating has a thin layer of gold evaporated 
upon its surface to minimize losses. 

Finally, after the FP cavity, there are post optics which take the collimated 
light beam and converge it upon the detector cavity. A AgCl filter is positioned at 
the focus just outside of the detector light cone. 

The spectrometer itself is only 10 inches high and 5.5 inches in diameter. It 
is actually somewhat smaller than the grating spectrometer described in chapter II. 
All of the mirrors and the grating are heat sunk in order to facilitate the cooling 
down time and to also make sure that all parts of the instrument are cold in order 
to insure a minimum of background radiation. 

As one cm tell, the package is quite compact. The entire dewar is less than 30 
inches in length with a diameter of 6 inches. The weight is less than 30 pounds. This 
configuration allows a liquid helium supply which lasts for roughly 20 hours. With 
the cryogen technology that has been used for IRAS, it is most certainly possible to 
equip a small telescope with a high-resolution Fabry-Perot of this type and launch 
it in a space vehicle. Short or long term missions comprised of individually targeted 
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spectr&l line observations or even whole sky surveys are possible with such a rugged 
instrument. 


3.5.2 The Meshes 

Many groups have constructed infrared FP systems using wire grids for their 
highly reflective plates. There are several ways in which the grids can be fabricated. 
Grids that are in the form of free standing screens are the most common. These 
sheets of mesh may be stretched flat across a ring and be made into an etalon plate. 

Other methods of constructing meshes is to deposit a mesh structure onto an 
optically flat transparent medium such as crystal quartz or silicon. The deposition 
can be done as either a positive or a negative. In either case, one will have a highly 
reflecting plate for a certain spectral regime. 

For grid deposition, we will only consider the positive deposition. Periodic 
structures in the form of crossed dipoles or tripoles have also been used. For more 
information concerning these, see Rhoads, Damon, and Munk (1982); Durschlag 
and Temple (1981); and Tomaselli et al. (1981). 

For the Hinge, free standing metal meshes were first considered for use. Since 
the mesh plane of the Hinge is defined by the mesh/grating holder (MGH), it is 
necessary to first mount the mesh on a mesh ring of some sort in order to maintain 
flatness, and then press that against the MGH. The mesh rings used are like those 
shown in figure 3.7. 

The mesh is made of nickel and is manufactured by Buckbee-Mears Co. in St. 
Paul, Minnesota. Nickel was chosen over copper due to its better cycling capability. 
Copper tends to not cycle well cryogenically and in time will become stretched. 

The mesh rings must define the flatness of the mesh. In order to accomplish 
this, the ring (figure 3.7) is constructed with an optically flat surface over which 
the mesh is stretched. A beveled outer surface gives one room to glue the mesh to 
the ring. The bonding agent is cyanoacrylate, also known as ‘Superglue’. Normal 
epoxies do not cycle well, but there are epoxies for low temperature work. The 
advantage to cyanoacrylate is that is can be readily dissolved with acetone to remove 
a mesh which is torn or improperly made. This is extremely helpful to the neophyte 
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FIGURE 3.7 Drawing of the mesh holders. These holders are constructed of 
stainless steel, the mesh is nickel. The narrow, flat annulus on the top of the holder 
is ground and polished to 1/ 100th of a wave, and then the mesh is stretched over 
the top and glued down on the wide beveled annulus. These free-standing meshes 
were later abandoned in favor of aluminum mesh photo-deposited on optically flat 
crystal quartz. 
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mesh-maker. 

Once the mesh is stretched and glued onto the ring, the ring may then be 
inverted so that the mesh is sandwiched between the MGH and the ring. The ring 
must be made of a material that has a lower coefficient of thermal expansion than 
Ni, else the ring will contract more than the mesh upon cooldown and the surface 
will no longer be flat. Stainless steel was chosen as the ring material. Since the 
mesh is sandwiched between dissimilar metals (stainless steel and aluminum) , it is 
possible to tear the mesh upon cooldown. A torn mesh can sag and the FP action 
can be lost. Care must therefore be taken in order to insure that the MGH surface 
is smooth and fiat, to minimize the chance of snagging. 

At first, #400 and #500 meshes were used. (#400 mesh has 400 lines per 
inch.) The results seemed to suggest that the Hinge had little Fabry-Perot action; 
i.e., it was acting like a grating instrument with two neutral density filters installed. 
The resulting instrument profile was almost as wide as that of the grating alone. 
Sagging of the meshes was suspected. 

It was decided that the sagging could be a real problem, and the only clear 
alternative to using metal mesh was to use mesh deposited on quartz or silicon. The 
transmission of crystal quartz in the FIR is roughly 75% for a 1mm thick sample. 
Fused quartz has a much lower transmission (Kimmitt 1970). 

The properties of a mesh are altered when the mesh is placed on a substrate. 
The most significant change is that the diffraction edge is shifted to longer wave- 
lengths by a factor that is equal to the index of refraction of the subtrate (Durschlag 
and DeTemple 1981). In other words, the spacing of the mesh is in effect increased 
by a factor of n, where n is the substrate’s index of refraction. 

The mesh that we are using was photodeposited on optically flat crystal 
quartz by Nellie Whetten at the National Research and Resource Facility for Sub- 
micron Structures here at Cornell. The grid is made of aluminum, with spacings of 
25 n m and line widths of 10 fim. The index of refraction for crystal quartz at 2 K 
is 2.120 at 150 ^m (Loewenstein, Smith, and Morgan 1973). These plates are then 
equivalent to #480 mesh. The quartz is slightly wedged to prevent unwanted FP 
effects within the crystal itself. 



The radiation that is reflected by the mesh is shifted in phase by an amount 
that is almost equal to it radians. It is, however, closer to 0.95 tt, owing to non-zero 
transmittance. Therefore, this must be taken into account in the placement of the 
meshes; and the mesh planes are adjusted accordingly. 

3.6 Error Analysis 

Associated with any instrument are the complementary errors that exist due 
to imperfections in either design or construction. The Hinge has several sources of 
errors which could almost totally frustrate the Fabry-Perot action if the errors were 
large enough. In this section, we shall consider these sources of error said see how 
significant they are. 

One source of error is that the actual mesh plane is displaced from where it 
should be by an amount 6. We will assume that the displacement is perpendicular 
to the plane. To understand how much the displacement will affect resolution, we 
need to consider each illuminated ruling of the grating as a part of a microetalon. 
There are N such microetalons. Let us first consider the situation when 6 = 0. For 
one edge of the beam, the distance between the plates is spanned by Ni wavelengths 
of radiation with wavelength Ao* At the other edge of the beam, there are N 2 such 
wavelengths. 

When the mesh plane is displaced by a distance 6, then each microetalon 
will transmit a slightly different wavelength. The range in wavelengths can easily 
be calculated. For the Ni^ microetalon, the transmitted wavelength is 

Ni Ai = Ni A 0 + 6 . (3.35) 

Likewise, for the microetalon, 

N 2 A 2 = N 2 A 0 + 6 . (3.36) 

We can therefore estimate the bandwidth by finding the difference between Ai and 
A 2 . This is 
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At 153 /im, Ni = 700 and N 2 = 892 for a beam 1.7 cm in diameter. If we want the 
resolution to be at least 10000, then 6 must be 2.0 mils or less (lmil equals 0.001 
inch). 

Another source of error is having the plane of the rulings not intersect the 
mesh planes along the rotation axis, but displaced by an amount 6. It can be 
shown that the displacement of the grating will produce a shift in the spacing of 
the microetalons, 6a. This shift is independent of position on the grating and is 
given by 

6a = 6 (cosa + cos/? - ^ , (3.38) 

\ cosoe / 

where a and 0 are the angles of incidence and diffraction, respectively. As above, 
we have the following range of wavelengths that are transmitted: 


Ni Ai = Ni A 0 - 6a 


N 2 A 2 = N 2 Aq — 6a , 


(3.39) 


from which we find 

AA = “’(£-£)• (3 - 40) 
Again, for Aq = 153 nm and a resolution of 10 4 , we find that 6a = 2.0 mils, or 
the grating may be displaced by as much as 4.2 mils. This is well within ordinary 
machining capabilities, as is the mesh displacement criterion. 

The grating offers two other sources of error. The first is if the rulings are 
displaced along the grating surface by an amount 6. That is to say, the rulings 
are not an integral number of spacings from the pivot axis. We measure from the 
center of each ruling to determine distance. The displacement will cause a range of 
wavelengths to be transmitted, since the difference in spacing between each succesive 
microetalon is no longer constant. Again, with similar analysis as above, we find 
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(3.41) 


where g is the ruling spacing (96 rulings per centimeter) . For a resolution of 10000 
at 153 fim, we have a machining tolerance of 1.3 mils. 


The other source of error from the grating is a side-to-side tilt. This situation 
will result in the beam walking off of the plates. This walkoff alone will limit the 
resolution. However, along with the walkoff, each succeeding bounce off the grating 
will result in a longer cavity spacing. The walkoff will lead to a lower transmission 
for obvious reasons. If we assume that there are now M traversals in the cavity 
before the beam walks off, then for each microetalon, we have from (3.2) 


I = Iq(1 — ii M ) 2 


1 + sin2 (^)' 

■ 1 + sin2 (f ) J ’ 


(3.42) 


where 4> is the phase delay as given in (3.1) and (3.30). The factor on the right 
will determine the resolution. For the canonical resolution of 10000 at 153 /*m, we 
find that the number of traversals within the cavity is relatively independent of the 
reflectivity. Six bounces, or twelve reflections off the grating, are necessary. Let ip 
be the tilt angle of the grating. With each reflection off the grating, the deviation 
of the radiation will double. After 12 reflections, we find that the deviation is 
«32V» meters. The maximum deviation that cam be tolerated is roughly 0.4 inch, 
allowing for the fact that the mesh is oversized. This implies that ip must be roughly 
3 x 10 -4 rad. Over the span of the grating width, 1.5 inches, this tolerance means 
that the grating axis must be positioned within 0.48 mil. The post optics and the 
exit aperture will dictate that this tolerance be even more stringent, and we should 
expect to be able to place the grating axis within 0.25 mil of where it should be. 

Each reflection off the tilted grating will produce a somewhat longer path 
length than the previous reflection. This longer path length will tend to spread the 
FWHM of the transmitted peak. To keep the spread under one part in 10000, this 
means that the grating must not be tilted by more than 0.001 rad, or a machining 
error of no more than 1.5 mils. This constraint is much more relaxed than that for 
the previous effect, which means that this is the less important of the two effects 
that a tilted grating will produce. 

The final set of errors that can arise from machining imperfections is due 
to the meshes being tilted. There are two directions of tilt, along the direction of 
dispersion and perpendicular to it. We shall first consider the latter. 
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If the mesh is tilted side-to-side, then we have the same problem that we had 
before with the grating. However, this time the radiation hits the mesh only once 
per bounce as opposed to hitting the grating twice per traversal through the cavity. 
This suggests that the tilt error can be no more than 1.6 mils in order to prevent 
the beam from walking off the mesh. However, again, the exit aperture will dictate 
a tolerance which is about a factor of 2 more stringent. 

For a tilt along the direction of dispersion, the tolerances are approximately 
a factor of 10 smaller. This is due to the fact that a small change in angle will result 
in a large movement along the grating owing to the steep diffraction angles. The 
angular error grows with each bounce as well. At 153 j*m, a tilt of 2 x 10~ 4 rad in 
one of the meshes will cause the beam to walk across the radius of the second mesh. 
Since the mesh diameter is 1.2 inches, this corresponds to a tolerance of better than 
0.24 mil. 

The machining constraints that are imposed by these errors are not very 
demanding in general. There are a few exceptions where tolerances of less than a 
mil are necessary, but these tolerances are still not beyond the skills and techniques 
of a very good machinist. Because of this, the machining costs of the Hinge were very 
reasonable. In fact, they were much less than the price for the ruling of the grating, 
though that price was dictated by the need to have the grating plane accurately 
pass through the hinge axis. 

3.7 Simulations 
3.7.1 Programming Considerations 

In order to determine the instrument profile, it is necessary to know what 
(3.34) looks like. There is no straightforward analytical solution to this equation, 
so numerical computation is required. 

There are several ways to calculate (3.34). In practice, it is desirable to have 
a minimum number of calculations, especially trigonometric functions. Written as 
is, (3.34) requires 2M trigonometric calculations, where M is the number of rulings 
illuminated by the beam. This number is ~ 200. We can reduce this by a factor of 



80 


at least two or more. 

In order to do that, we go back to (3.32) and (3.33). Instead of first sum- 
ming over the number of reflections, we should sum over the number of rulings 
illuminated. Equation (3.33) then becomes 

oo M— l 

E t ot = E 0 T £ £ fE» e -i(n+§)*m . (3 43) 

n=0 m=0 


We can rewrite <f> m as 


47 r 

<j>m = 9 { Ni + m) (sina + sin/?) 


— C + m * , 


(3.44) 


where 


47 r 

? = — g Ni (sina + sin/?) 
47T 

/c = — g (sina -I- sin/?) . 

A 


(3.45) 


This gives us the following expression for the total electric field: 


£** = Eq T±R* 3in| i" + i. ) jfe 
“o sin[(n+ 2 )f] 


(3.46) 


We calculate the intensity to be 
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In practice, one need only sum over a small and finite number of reflections instead 
of over am infinity of them. It turns out that the finesse is a good measure of how 
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many bounces over which one needs to sum. If we consider N bounces, then one 
sees that we need only calculate 4N trigonometric functions. For a finesse of 30, this 
method gives the computer almost a factor of two fewer trigonometric functions to 
calculate. 

From the dimensions of the Hinge and its optical configuration, we can cal- 
culate Ni and M. From a and /?, we can then figure out f and k. It is then a simple 
matter to perform the sums in (3.47) and to calculate the intensity as a function of 
wavelength for meshes with reflectivity R. 

3.7.2 Results 

The transmission profiles of the Hinge have been calculated for various reflec- 
tivities. Figure 3.8 shows two profiles, one for R = 0.7 and the second for R = 0.9. 
These were calculated for the grating tuned to 154.6179 ^xm radiation. As one can 
see, there are several interesting features of these profiles. 

The first thing to notice is that there is a narrow central spike which is due to 
the combined constructive interference from all of the microetalons. The sidelobes 
are suppressed by the action of the grating. As one can tell, higher mesh reflectivity 
yields better sidelobe suppression. This is because a higher reflectivity increases the 
number of bounces within the FP cavity and causes the effective bandpass of the 
grating to decrease. 

The sidelobes are spaced approximately 0.1 /xm apart. One can calculate the 
resolution from the simulation profiles and find an equivalent ordinary FP system. 
At this wavelength, we find a good fit for a plate spacing of 12.28 cm. The actual 
pathlength between plates for the central portion of the beam is 12.44 cm, but this 
pathlength ranges between 10.94 and 13.95 cm across the size of the beam. The 
equivalent FP system will result in a Qequiv* 8 0.10 /xm, which is what we observe in 
figure 3.8. 

Since there is a range of orders going through the etalon cavity, there is 
also a range of Q values. For the simulations at 154.6179/xm, Q will range from 
0.086 to 0.109 /xm. This means that the sidelobes should be much broader than the 
central peak. Since the Q values cover a relatively broad range, this can also explain 
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FIGURE 3.8 Theoretical tranmission profiles of the Hinge for mesh reflectivities 
of 70% (top) and 90% (bottom). 
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the apparent bifurcation of some of the sidelobes as different parts of the beam 
come into or go out of phase with other parts. 

The resolution of the system closely follows that for the equivalent FP system. 
This is shown in figure 3.9 as a function of mesh reflectivity. For reflectivities 
approaching 95%, the resolution nears 10 s . As the reflectivity goes down below 
20%, the resolution of the Hinge drops below that for the equivalent FP system. 
The grating profile dominates at low reflectivities. 

If we have high resolution, it is important that the sidelobes are suppressed 
in order to avoid confusion. In this system, higher reflectivities not only give higher 
resolution, but automatically give better sidelobe rejection due to the increasing 
number of times that the beam hits the grating. For 70% reflectivity, the first 
sidelobe is 67% of the central peak. At 90% reflectivity, this number is down to less 
than 20%. 

The grating is not 100% efficient, as was assumed in the above simulations. 
If the grating is 90% efficient in the one polarization, it acts as a neutral density 
filter and effectively lowers the reflectivity of the meshes by — 10%. 

3.7.3 Errors 

Errors due to machining imperfections have been incorporated into the trans- 
mission profile programs. This is done by simply following the beam of each mi- 
croetalon separately as it encounters tilted meshes or displaced meshes and gratings. 
We have already discussed the various sources of error in section 3.6. 

Figure 3.10 is a transmission profile with one mesh that is tilted along the 
direction of dispersion by 0.10 mil. It is tilted so that the beam will walk down the 
grating. What we see is that the total transmission is down by a factor of 4 or 
more, and the peaks have shifted to somewhat longer wavelengths. The resolution 
is also degraded, as well as the symmetry of the profile. 

It is evident that the transmission profile is very sensitive to tilt errors. These 
errors can arise in several ways, including a misalignment of the optics. Hence, it 
is important to carefully align the optics in order to minimize any degradation to 
the system. 
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FIGURE 3 . 9 Theoretical resolution as a function of mesh reflectivity for the Hinge. 
The solid curve represents the resolution for an equivalent ordinary FP system. At 
low reflectivities, the grating profile dominates, causing lower resolution for the 
Hinge as opposed to the equivalent system. 



85 



Wavelength (in >tm) 


FIGURE 3. 10 Transmission profile of the Hinge if one of the meshes is tilted by 
0.01 radian. The transmission is decreased due to walkoff problems, the transmitted 
radiation is shifted towards longer wavelengths (as directed by the tilt), and the 
profiles are widened and asymmetric. 
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Displacement shifts of the grating do not affect the profiles very much. Basi- 
cally, the profile will shift in wavelength depending on whether the cavity is short- 
ened or lengthened by the grating displacement. 

3.8 Laboratory Results 

3.8.1 Spectra Through the Interferometer 

In order to determine the actual instrumental profile, we made two sets of 
measurements in the laboratory. The first set was made using the interferometer. 
The second was using a narrow line source. 

Meshes deposited on crystal quartz were installed into the instrument. These 
meshes have reflectivities on the order of 85%. With the interferometer, we were 
able to determine the Hinge profile with a resolution of over 4000 at 153 /xm. Figure 
3.11 shows the result. Each bin of the spectrum is 0.038 jxm wide. To within a 
fraction of that, the central peak is centered in one bin and is less than one bin 
wide. We know from this that the resolution is better than 4000, but we do not 
know how much better. 

Unfortunately, we cannot get more resolution out of the interferometer. The 
spectrum is noisy out in the wings, primarily due to the fact that the system sen- 
sitivity was not very good at this stage. This was primarily caused by the detector 
losing stress and the quantum efficiency plummetting at these long wavelengths. 

During this same time, we concentrated our efforts at obtaining information 
at 142 ^m, where the signal was strongest. This spectrum is shown in figure 3.12. 
As one can tell, the sidelobes are about 30% of the signal. This is consistent with 
a reflectivity of 85%. The wings of the spectrum are also close to 0. 

3.8.2 Absorption Spectrum of CO 

The interferometer allows us to look at the gross characteristics of the Hinge. 
In the earlier tests that were run, it was noted that there was a short wavelength 
leak (about 20 /xm less than the central peak). In order to combat this, the CsBr salt 
filter was installed to purge any radiation below 138 /xm from entering. The cause 
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FIGURE 3 . 1 1 Laboratory transmission profile of the Hinge at 153 nm. These data 
were obtained using the interferometer to determine the instrument profile of the 
Hinge. The instrumental resolution of the interferometer was about 4000. Clearly, 
the Hinge’s bandpass is narrower than what the interferometer can measure. 
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FIGURE 3.12 Laboratory transmission profile of the Hinge at 142 /*m. These 
data were obtained in the same manner as the 153 nm data in figure 3.11. Note the 
sidelobe structure. 
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of the leak is most likely radiation that is bouncing off of one of the mesh/grating 
holder apertures and hitting the grating at a somewhat shallower angle. 

In order to determine the actual resolution, it is necessary to have a source 
of narrow line emission or absorption. This is accomplished by placing a cell of CO 
gas in the beam. There is a transition from the J = 16 — * 17 rotational levels of 
CO at 153.2667 /urn. This transition would be seen in absorption. By stepping the 
grating, it is possible to sweep the instrumental profile along this narrow absorption 
line and determine the transmission profile. 

The Hinge can be stepped in bursts of 3.9 x 10 4 pm. This was determined 
using the interferometer to calibrate the Hinge from 140 pm to past 153 /xm. By de- 
termining the halfwidth of the absorption feature, we cam determine the resolution. 

Figure 3.13 is the result of this endeavor. The data were fitted by a Gaus- 
sian (although a Lorentzian should be the true shape). The FWHM is found to 
be 0.00716 fim, or the resolution is 21400. This is also consistent with a mesh 
reflectivity of 85%, if the grating efficiency is 93%. This is most likely the case. 

Since the pressure in the cell was close to atmospheric at the time, it is quite 
possible that the actual resolution of the system is better. The cell we were using 
had a slow leetk, and in order to integrate for the length of time we needed for a 
good signal-to-noise ratio, it was easiest to simply keep the cell at ambient pressure. 
Due to this, the CO line may be pressure broadened. Owing to a lack of data on the 
collision rate of N 2 -CO for this transition, it is not possible to accurately calculate 
the pressure broadening. Estimates of the N 2 -CO collision rates from tabulated 
He-CO collision rates yield a pressure broadening width which is over a factor of 
10 larger them that observed. 


3.9 Summary 

We have demonstrated that the Hinge is a viable instrument capable of pro- 
ducing high resolution (R « 21000) spectral observations in the FIR. The instru- 
ment is compact, highly rugged, and reliable. We have cooled it down several times, 
realigning the optics each time, and have found similar results with every cooldown. 
The Hinge is also stable in the sense that the grating may be tuned to within 1/3 
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FIGURE 3 . 1 3 Laboratory tr ansmis sion profile of the Hinge at 153 /xm. These data 
were obtained by observing a narrow absorption feature due to CO and sweeping 
the grating of the Hinge through the feature. Assuming that the absorption line is 
not significantly pressure broadened, teh resolution of the Hinge is roughly 21500. 
The solid line is a profile fit to the data. 




to 1/2 of a resolution element. 

An instrument such as this is very attractive for use in space missions where 
a Fabry-Perot spectrometer can be used. This instrument is very simple in that 
it has only one moving part. That part is a grating and the motion is a simple 
rotation. There are no needless complexities that are associated with regular FPS 
systems where the movable plate must be kept parallel to the fixed plate. 

The sidelobe rejection is adequate in the Hinge, and rivals that found in other 
FIR FPS systems. With this instrument, FIR emission and absorption lines can 
be resolved, line profiles can be determined, and a greater understanding of FIR 
sources can be achieved. 



IV. Molecular Rotational Lines and Radiative Transfer 

4.1 Introduction 

Within many astronomical systems, molecules are present, sometimes in 
great abundance. They are remarkably durable and exist in environments which 
are not conducive to either their formation or survival (Dalgarno 1985a). Owing to 
their complex energy level structure, they serve as excellent probes of the physical 
conditions that exist within regions in which they are found. Molecular emission 
or absorption lines have been detected from the interstellar gas in external galaxies 
and from molecular clouds in our own Galaxy. These observations have helped in 
the understanding of these sources. Other areas of astronomical study which have 
benefited from molecular spectroscopy include masers, circumstellar shells, proto- 
stars, Hll regions, planetary nebulae, Herbig-Haro objects, sun spots, comets, and 
planetary atmospheres. We are beginning to gain a much fuller overview of star 
formation regions, stellar mass loss, galactic structure, and planetary atmosphere 
dynamics from these observations. 

The most abundant molecule almost certainly is H 2 . Most of this cannot be 
directly observed, and it is only through other molecular studies (mainly CO) that 
an estimate of the H 2 abundance is known. From this, one concludes that molecular 
hydrogen makes up 10-20% of the mass of the inner disk of the Galaxy (Shull and 
Beckwith 1982) and is probably responsible for 50% of the hydrogen gas content 
beyond the solar orbit (Li, Riley, and Wolfendale 1983). 

Other molecules which have been found in abundance include CO, OH, CN, 
CS, HCN, HCO, and NH 3 , to name a few. As of 1984, 64 molecules and 46 isotopic 
variants had been observed in a host of sources. These range from the simplest 
diatomics to cyanodecapentayne (HCnN) (Rydbeck and Hjamlarson 1985). The 
universe is rich in molecular chemistry. 

Different transitions, which occur in different parts of the electromagnetic 
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spectrum (ranging from the ultraviolet to radio), will eminate from different parts 
of the region in which the molecules are found. The molecules serve an ancillary 
role as probes. Often, they will control the structures, both thermal and ionization, 
of the systems in which they are found. Their presence can dramatically alter the 
dynamical instabilities, and hence the evolution, of these systems. 

In this thesis, we will be mainly concerned with the shocked region of the 
Orion Molecular Cloud (OMC-1) surrounding the Becklin-Neugebauer object and 
the Kleinmann-Low nebula (BN/KL). The shocked region is most likely a result of a 
massive stellar wind from very young stars plowing into the surrounding molecular 
cloud. The gas in the shock is hot (~ 1000 K) resulting in significant populations 
for the higher rotational levels of molecules such as CO and OH. Next to H 2 , these 
molecules are some of the most abundant in the region. 

The major portion of the radiation from these molecules will result from 
rotational transitions in the FIR. The millimetric lines are suppressed owing to 
the small values for the spontaneous transition coefficients, and the mid- and near- 
infrared lines are suppressed owing to the small level populations. For CO, the peak 
intensity occurs for J— 15, which is near 160 ^m (McKee et al. 1982). Theoretical 
studies of the shocked region in Orion suggest that CO emission should be significant 
for a wide range of J values, probing a temperature range of a few degrees Kelvin 
to a few thousand degrees Kelvin (Draine and Roberge 1984; Chernoff, Hollenbach, 
and McKee 1982; Draine and Roberge 1982). The FIR observations of CO and 
other molecules cam enhance our understanding of these shocked regions. 

By observing many rotational emission lines of CO and OH, it is possible to 
piece together a fairly constrained model of the emitting region. Since different lines 
probe different parts of the shocked region, our understanding of the area’s overall 
structure can be greatly enhanced. It is our goal here to undertake this task. A 
significant part of the groundwork for the understanding of the CO emission lines has 
been laid down by Draine and Roberge (1984); Chernoff, Hollenbach, and McKee 
(1982); Draine and Roberge (1982); McKee et al. (1982); and Watson (1982). Not 
until in the past several years have enough OH lines been measured to facilitate an 
analysis similar to that tried with the CO line intensities. 
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4.2 Molecular Physics 
4.2.1 Theory and Nomenclature 

Since this thesis deals primarily with CO and OH, we shall limit ourselves to 
the physics of diatomic molecules alone. Interested readers may refer to Herzberg 
(1945) and Townes and Schawlow (1975) for more information concerning poly- 
atomic molecules. We shall generally follow Herzberg (1950) for the discussion of 
diatomic molecules. 

For a diatomic molecule, we first consider the case of a rigid rotator consisting 
of two masses, m\ and m 2 , separated by a distance r. Let us assume for the moment 
that the molecule is spinning about an axis perpendicular to the internuclear (figure) 
sods. The moment of inertia is given by 


Ib = V ' 2 . 


(4.1) 


where n is the reduced mass of the system. The rotational energy is then written 
as 


E = 



(4.2) 


where u is the angular velocity of rotation. Eigenstates of the Schrodinger equation 
exist for 


h 2 J[J + 1) 

8tt2/ b 


(4.3) 


where J is a positive integer. The total angular momentum, excluding spin, can be 
written as 


|J| = ■v/J(7+I)^ . 


(4.4) 


Equation (4.3) can be rewritten as 


E = BJ{J + 1) cm -1 , 


(4.5) 


B = 


h 

8tt 2 c Iq 


where 


(4.6) 
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B is known as the rotational constant. 

In general, the total angular momentum of a rotator is not necessarily per- 
pendicular to the figure axis since the electrons will contribute to a nonzero moment 
of inertia about the line joining the nuclei. This moment of inertia is known as / A . 
In most cases, J A is much smaller than Iq. However, the corresponding angu- 
lar momenta are nearly the same because the electrons ;orbit the nuclei at higher 
frequency. 

The component of J that is parallel to the figure axis is called A. This is a 
measure of the electronic angular momentum, and as such, we have 

W = A £ • (4.7) 

Here, A is the angular momentum quantum number of the electrons about the figure 
axis. If we let L be the total electronic orbital angular momentum, then one can 
show that the component of L about the figure axis is a constant of the motion. L 
will precess about this axis with a constant component Ml(A/27t), where 

M l =L,L-1, (4.8) 

Hence, from (4.7) and the above discussion, we have 

A = \M l \ . (4.9) 

The values of A can range from 0, 1, 2, • • • , L. For A = 0, 1, 2, 3, • • *, the corresponding 
molecular states are designated as E, IT, A, $,•••, quite analogous to the atomic 
spectral nomenclature. 

Likewise, the component of J that is perpendicular to the figure axis is desig- 
nated N. It represents the rotation of the nuclei. Figure 4.1 shows a vector diagram 
of the diatomic molecule and the various momentum vectors. As can be seen, the 
direction of A cam be in one of two senses. For a given N, there are two modes 
of motion of the system corresponding to the two directions of A for each value of 
J. This double degeneracy is evident from (4.9), where for a non-zero A, Ml may 
have two values, A and —A. 
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FIGURE 4. 1 Vector diagram for the rigid rotor. The actual rotation is around J 
and is depicted by the curved arrow. N is the component of the rotation about an 
axis perpendicular to the intermolecular axis while A is that component about the 
intermolecular tods. After Herzberg (1950). 
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We can treat the electronic spin similarly. The component of the spin in the 
direction of the internuclear axis is given by E. If the resulting total spin is S, then 

E = 5,5 -v . (4.10) 

As opposed to A, E can be either positive or negative. 

The total angukir momentum along the internuclear axis is obtained by 
adding A and E. Since both of these vectors lie along the same axis, we only 
need to add algebraically. The sum is given by 

n = |A + E|. (4.11) 

The multiplicity of the state is given by 2 5 + 1. Hence, as in atomic nomenclature, 
the multiplicity is written as a superscript to the left of the molecular state desig- 
nation and (A + E) is a subscript to the right. For example, if A = 2 and 5 = 1, 
then the possible states are 3 A3, 3 A2, or 3 Aj. 

The symmetry parity of the rotational states depends on the value of A. For 
A = 0, the even J levels have positive parity, or the eigenfunctions do not change 
under an inversion. The odd J levels have negative parity. If A # 0, then each 
rotational level is doubly degenerate, one eigenfunction being positive, the other 
negative. The selection rule for transitions associated with dipole radiation is this: 
transitions between like parity states are forbidden. 

When we include spin in the coupling mechanisms, we find that J becomes 
a resultant of the molecular rotation and the total electronic angular momentum. 
Hence, we may write 

j = n,n + i,n + 2, . . . . (4.12) 

With a rotation of the molecule, the states that are doubly degenerate owing to a 
non-zero A will split in energy. This splitting is known as A-doubling. 

One can also observe hyperfine structure, produced by nuclear spin, in the 
rotational spectrum of diatomic molecules. If we let T be the total nuclear spin, 
then the resulting quantum numbers are 


F = J + T,J + T -T\ . 


(4.13) 
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The coupling strength between different angular momenta and rotations is 
dependent upon the molecule in question and how fast it is rotating. For this 
reason, the above analysis does not always hold. There are five major cases, known 
as Hund’s cases (a) through (e), whereby different couplings hold for each case. We 
have followed Hund’s case (a), which describes fairly well the OH molecule. 

4.2.2 Non-rigid Rotors and Rotation Constants 

So far, we have been treating diatomic molecules as rigid rotators. This is, 
however, only an approximation. We know that molecules also display vibrational 
spectra, which means that the distance between the two atoms can change. 

If one considers the diatomic molecule to consist of two masses connected by 
a stiff spring, one sees that for faster rotation, the two masses will be pulled further 
apart and the moment of inertia, /q, will increase, decreasing the rotational energy. 

To accurately fit the rotational spectra of diatomic molecules, we write the 
energy of a level as a sum of the vibrational contribution (which is a function of 
the vibrational quantum number, v) and the rotational contribution (a function of 
J). Then we can expand in the appropriate variables to find (Mantz et al. 1975; 
Herzberg 1950) 


E = u e (v + j) - u e x e (v + ^) 2 + u e y e (v + i) 3 

+ J(J + l)[i?e - cc e (v + -) + 7e(v + -) 2 - 6 e (v + -) 3 H ] 

+ J 2 (J + l) 2 [-^e — /? e (t> + ~) + *e(t> + ^ ) 2 ~ 

+ J\J + l) 3 [/f e - Vt{v + i) + •••] + J\J + l) 4 (L e -•••], (4.14) 

where the subscript ‘e’ refers to the equilibrium values. The different variables 
we defined in the standard references on molecular spectra and won’t be redefined 
here. Terms in the top line represent the vibrational contribution. The other terms 
represent the contributions from the rotation of the molecule. This equation may 
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also be rewritten as (Townes and Schawlow 1975) 

£ = E y U (” + j)' ^ + !)' • (4-15) 

1J 

where the summation is carried out from 1 to infinity over j and 1. The coefficients, 
Yjj, are known as Dunham coefficients. In practice, only the lower Dunham coef- 
ficients are necessary. Either set of coefficients may be experimentally determined 
and used to calculate energy levels and transition frequencies. 

For further information on the physics of molecular spectra, see Herzberg 
(1950) and Townes and Schawlow (1975). 

4.2.3 Selection Rules and Einstein A Coefficients 

A molecule cannot undergo an arbitrary radiative transition from one level 
to another. As for atoms, there are allowed and forbidden transitions for dipole 
radiation. The selection rules for the allowed transitions are given as follows: 

AJ = 0, ±1 
AF = 0, ±1 
Av = 0, ±1 
AA = 0,±1 

(4-16) 


The transition probability is represented by the Einstein A coefficient, which 
gives the lifetime of an atom or molecule in the excited state before it decays to 
a specific lower energy state. It is well known that the transition probability of 
spontaneous emission is 
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(4.17) 
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where t / nm is the energy difference between levels n and m in wavenumbers, and 
R 11 ™ is the dipole matrix element. If A J = — 1, the transition to a lower energy 
level occurs by changing J by 1; and Townes and Schawlow (1975) show that 

|RDra|2= " 2 ijTi' (4 - 18 > 

where \i is the permanent dipole moment of the molecule. This dipole moment 
depends, to some extent, on the molecular state (£, IT, etc.). For the most part, we 
shall only be concerned with rotational levels within the ground vibrational state 
of the ground electronic state. The energy differences between rotational states is 
on the order of a few meV; between vibrational states, ~1 eV; and between different 
electronic states, several to several tens of electron volts. For the temperatures with 
which we shall be dealing, the molecules are, in general, in their ground vibrational 
and electronic states. 


4.3 The Molecules 

4.3.1 CO 

The major thrust of this work is to understand the physical conditions in 
the shock region of OMC-1 that is the source of CO and OH emission from highly 
excited rotational states. This is done by collecting data on a number of different 
transitions and then using radiative transfer to determine temperatures and densi- 
ties. In order to do this, we shall anticipate the transport theory and present here 
the sources of the molecular data necessary to determine such quantities as level 
energies, spontaneous transition probabilities, collisional cross-sections, etc. 

The CO molecule has a straightforward rotational spectrum structure. The 
ground molecular state is a *£, with A = 0 and S = 0. This means that J is 
identical to the nuclear angular momentum, N. Figure 4.2 shows am energy level 
diagram for the X E (t; = 0) rotational lines. The values for the energies were derived 
from the rotational and Dunham coefficient data of Mantz et al. (1975) and Todd 
et al. (1976). 

To find the Einstein A coefficient, we may use the results of (4.5), (4.17), 
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FIGURE 4.2 Rotational energy level diagram for CO in its ground electronic and 
vibrational states. Transitions which have been observed in OMC-1 are indicated 
by arrows along with their respective wavelengths in microns. 
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and (4.18). Equation (4.5) is a good approximation for finding the energy levels. 
Combining all of these, we find that the spontaneous emission rate is 




512* 4 „3,,2 J* 
3h*c 3 * 2J+1 



(4.19) 


for B in units of ergs. If J is large, (4.19) may be approximated by 


Aj-+j - 1 = Aq J 3 , 


(4.20) 


where Aq is found to be 1.118 x 10 -7 s -1 . The dipole moment of CO is 0.112 debye 
(D) (Muenter 1975). 

Transitions among levels can take place not only radiatively, but through 
collisional processes as well. In the interstellar medium, most collisions are with 
H 2 , the most abundant constituent. Collisional excitation rates for CO by H 2 have 
been calculated by Green and Chapman (1978) and McKee et al. (1982). Since 
the potential energy surface (PES) was not known very well at that time, these 
rates were extrapolated from He-CO collision rates, being multiplied by a factor 
dependent on the velocity difference between H 2 and He. Green and Thaddeus 
(1976) argue that scaling by the velocity ratio works for H 2 , but only when the 
potential surface is roughly spherical. This describes para-H 2 (even j values), but 
it is possible that ortho-H 2 (odd j values) may be abundant as well, and then 
the scaling arguments lose ground. However, owing to the small dipole moment 
of CO, the difference between ortho- and para-H 2 may be small. Recently, Flower 
and Launay (1985) have calculated para-H 2 ~CO collision rates which tend to agree 
with those of Green and Thaddeus (1976). Schinke et al. (1985) have recently 
incorporated a more accurate H 2 -CO PES in calculating the collision rates. Their 
rates differ significantly from those of McKee’s group for the large AJ transitions. 
Experimental results at low temperatures would be very useful. It should also be 
noted that most of the calculations are for temperatures of 300 K or less. This is 
not always the case, especially in the shocked regions where H 2 vibrational spectra 
indicate T ~10 3 K. 
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4.3.2 OH 

The ground electronic state of OH is a 2 n state, where A = 1 and 5 = 1/2. 
H can take on two values, 1/2 or 3/2. This leads to a splitting of the levels into two 
ladders, one for fl = 1/2 and the other for H = 3/2. Figure 4.3 is a schematic of 
these two ladders, showing the A-doubling due to the parity. The hyperfine levels 
within each parity level are not shown. 

The Einstein A coefficients for the various A-doubling transitions within a 
rotational level have been calculated by Destombes et al. (1977). Brown et al. (1982) 
have published similar data for the various rotational transitions including hyperfine 
structure. The energy levels of OH found by Brown’s group are more reliable since 
they were measured by laser magnetic resonance. The transition probabilities are 
consistent between the two groups. 

Collisional excitation rates for H 2 -OH interactions have not, until recently, 
been known very well. Dewangan and Flower (1981, 1983) made an improved 
calculation over that of Bertojo, Cheung, and Townes (1976); however, recent ex- 
perimental results (Andresen, Hausler, and Lulf 1984) suggest that the theoretically 
expected cross-sections are significantly in error. Schinke and Andresen (1984) have 
recently performed another set of calculations for H 2 -OH collisional excitation, and 
their results match the trends in level populations observed experimentally. The 
actual values of the theoretical cross-sections, however, still do not agree entirely 
with experimental results. Again, these results assume that the hydrogen is in a 
para state. 


4.4 The Chemistry and Physics of Shocks 

In this section, we will discuss some of the physical and chemical properties 
of shocks. The treatment of shocks includes some very complex and often little- 
understood areas of study, and it is not within the framework of this thesis to 
discuss these in depth. We do wish to provide a general overview of shocks in the 
interstellar medium as they are understood today. This will follow the outline of 
Dalgarno (1985b). 
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FIGURE 4.3 Rotational energy level diagram of OH in its ground electronic and 
vibrational states. Transitions which have been observed in OMC-1 are indicated 
by arrows along with their respective wavelengths in microns. 
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The first question that must be addressed is the definition of a shock. A shock 
consists of a supersonic propagation of a disturbance in a fluid medium, usually a 
gas. Since the information concerning the disturbance cannot propagate upstream 
faster than the local speed of sound, the medium cannot anticipate and react to the 
disturbance until the shock actually arrives. Upon the arrival of the shock, the gas 
is compressed and heated. The boundary layer separating the upstream gas from 
the hot, compressed, shocked gas is known as the shock front. 

Shocked regions are associated with many different astronomical phenomena. 
They are most commonly linked to stars. In stellar formation processes, strong 
winds driven by mass loss and radiation pressure, will plow into the ambient parent 
gas cloud and cause a shock wave to form. Massive stars will continue to feed the 
shock wave throughout most of their lives via radiation pressure or stellar winds. 
Cataclysmic stellar deaths (supernovae) will produce enormous mass losses which 
also result in shock waves. Collisions of interstellar clouds are yet another source 
for shocks. 

Various factors affect the physics of the shock. The shocked material will 
be compressed and heated. This heating can result in dissociation and ionization, 
whereupon recombination will produce photons capable of ionizing and dissociating 
the upstream material. This will change the nature of the shock as it progresses 
through the cloud medium. Magnetic fields complicate further the shock physics. 
Magnetohydro dynamic (MHD) considerations suggest that magnetic precursors can 
preceed the shock and compress the material before the shock arrives. 

Much of the kinetic energy in the shock front is transformed into thermal 
energy. It is very simple to roughly determine the temperature, T a , of the shocked 
region for a shock velocity, v 8 , by equating the kinetic energy to the random thermal 
energy. For atomic hydrogen, this results in 

r 8 « 40vjK, (4.21) 

where v 8 is measured in km s _1 . If the shock is slow, the temperature will rise, 
but not enough to actually cause much dissociation. Instead, the gas will become 
warm and the chemical composition will evolve. For feist shocks, dissociation and 
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ionization will result, with new molecular populations being born out of the cooling 
downstream gas. 

The shock regions are sites of complex chemical and physical interactions. 
Dissociation is collisionally induced, but how fast molecular dissociation progresses 
depends upon density and radiative transfer. For low density clouds, a shock ve- 
locity of about 50 km s _1 will result in dissociating 90% of the H 2 . Since atomic 
hydrogen reacts very efficiently with CO to produce C and O, CO dissociates at 
v B at 40 km s -1 . Higher densities decrease the radiative stabilization of molecular 
hydrogen by increasing the probability of collisionally dissociating a molecule in a 
vibrationally excited state, resulting in nearly complete dissociation at shock veloc- 
ities of 25 km s~ . At this velocity, the CO will have a better chance of surviving 
owing to the lower shock temperature. When magnetic fields are present, the hy- 
drogen may not dissociate until velocities on the order of 45 km s -1 are reached 
(Draine, Roberge, and Dalgarno 1983). 

The chemistry of the shocked region is governed mainly by reactions with 
atomic and molecular hydrogen. Prasad and Huntress (1980) have listed over 1400 
reactions among 137 species. Improved reaction rates and overlooked reactions have 
been listed by Mitchell (1984) and references within for 113 chemical species. 

Once the shock has passed, the temperature of the region is raised and chem- 
ical reactions begin. These reactions, and the abundances of the resulting products, 
depend to a large extent on the physical parameters of the region, such as tempera- 
ture and density. The gas is primarily cooled by fine-structure transitions of O I and 
the rotational and vibrational lines of H 2 , CO, OH, and H 2 O. Hence, the thermal 
and chemical structures are interdependent. 

A schematic of the oxygen and carbon chemical reaction families is presented 
in figure 4.4. In brief, the oxygen will tend to combine with carbon to form CO 
or else with hydrogen to ultimately form H 2 O. Small amounts of O and OH will 
also be present, with OH enhancement present near the shock front (Draine and 
Roberge 1982). For more details on the chemistry and physics of shock regions, see 
Dalgarno (1985b); Draine (1985); Mitchell (1984); Draine, Roberge, and Dalgarno 
(1983); Chemoff, Hollenbach, and McKee (1982); Draine and Roberge (1982, 1984); 
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FIGURE 4.4 A schematic diagram of the carbon and oxygen chemistry that is 
believed to occur in star formation regions and molecular clouds. After Prasad and 
Huntress (1980). 
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Prasad and Huntress (1980); and references within. 


4.5 Radiative Transfer 

4.5.1 Definitions 


In just about any area of astrophysics, one is dealing with the transport 
properties of a medium. For molecular shocks, the interaction of the radiation field 
with the gas is very important in determining the chemical and thermal structure 
of the region. This interdependence requires a knowledge of radiative transfer. 

It is desirable that the terms for the description of radiation and its mea- 
surable quantities be properly defined. We will be using the terminology of Gray 
(1976) and Rybicki and Lightman (1979). 

We begin by considering a radiating surface of area dA. The specific intensity 
is a measure of the amount of radiant energy, dE v , between u and du, emitted at 
an angle, 6, from the normal to the surface, and contained within a solid angle <ffl 
in the time interval dt. This can be written as 


dE u 

cos 0 d A dCl dt du 


(4.22) 


Averaging over direction, we obtain the mean intensity, J u , which may be written 
as 

Jv = — [ I u dn . (4.23) 

47r J 

The integration is performed over the whole unit sphere centered on the area, dA. 

The flux is defined as the net energy flow across an area dA, in a time dt, 
and in a spectral range du, or 


jJEy 

v dA dt du ’ 


(4.24) 


where the integration is again performed over the unit sphere. Combining (4.22) 
and (4.24), we find that 

F u = f Iv cos 0 dCl . 


(4.25) 
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As radiation passes through a medium, emission and absorption processes 
may add or subtract energy from it. The absorption of energy is characterized by 
an absorption coefficient, a „, in units of cm - *. The change in the specific intensity 
is given by 

dly = —ctylyds , ( 4 . 26 ) 

where ds is the length of the medium through which the radiation has traveled. If 
n is the number density of absorbers, and a v is the cross-section of each absorber 
for radiation of frequency u , then it can be easily shown that 


Oly — TlOy ■ ( 4 . 27 ) 

For absorption processes, the change in the specific intensity is proportional to the 
intensity itself. 

The emission coefficient j v can be defined in much the same way as the 
absorption coefficient. In general, j u is independent of the specific intensity. We 
may write the change in intensity as 

dly = jy ds . ( 4 . 28 ) 

The units of the emission coefficient are energy per unit time per unit solid angle 
per unit volume per unit frequency interval. There axe two processes that can 
contribute to j v , either real emission or scattering. 

It is quite simple to incorporate the emission and absorption into the equation 
of radiative transfer, which is 


dly _ - . 

— — —Clyly+Jy • 


( 4 . 29 ) 


The absorption and emission of radiation is governed by quantum mechanical 
considerations. If we assume that we have a system which has two energy levels, 
1 and 2, then there are three processes that were identified by Einstein that can 
take place. The first is that of spontaneous emission. Here, the system begins in 
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the higher energy level, 2, makes a transition to level 1, and in so doing, emits 
a photon corresponding to the energy difference. This transition occurs without 
outside influence, and the transition probability rate is denoted by A 21 . The units 
are s _1 . 

If the system is in level 1, then there is the chance that it can absorb a 
photon of the appropriate energy to boost it into level 2. The probability rate 
for this process must be proportional to the photon density. We can write this 
probability as 

probability rate = B 12 / Ju du , (4.30) 

Jo 

where 4>{u) is the line profile and B\2 is the Einstein coefficient for absorption. 

The system can interact with an external radiation field in another way, 
and that is via the process of stimulated emission. In this process, the system 
is in the excited level and a photon of appropriate energy will interact with the 
system so as to induce it to emit a photon and jump to the lower energy level. The 
probability rate is defined much like that for absorption, and the Einstein coefficient 
for stimulated emission is written as f?2i- 

4.5.2 Relations and Special Cases 

In this section, we look in more detail at the equation of radiative transfer 
and several special cases. 

We start by assuming that we have pure emission and no absorption. In this 
case, a v — 0 and the solution to (4.29) is 

Iv{»)=Iv(*o)+ ['m*')**' , (4.31) 

J a 0 

where we integrate from some initial so to s. 

Next, we consider pure absorption without emission. In this case, j y = 0 and 
the solution for the equation of radiative transfer is 

Me) =Mso)e" Ma) , (4-32) 
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where 

M«) = f ap(s')ds' . (4.33) 

J » o 

The variable r„ is the optical depth, and it is a measure of the probability of a 
photon being absorbed within the medium. The greater the optical depth, the 
greater the chance for absorption. Typically, if the medium has an optical depth 
less than unity, it is considered to be transparent; whereas for optical depths of 
unity or greater, the medium is opaque. It should be apparent that 

7"|/ — a u ds . (4.34) 


If we divide (4.29) by a v , then we have 


dlv _ t i j* 
~7~~ — ~ i v t — • 
dT u CL%/ 


(4.35) 


The ratio of the emission coefficient to the absorption coefficient is known as the 
source function S v . This allows us to rewrite (4.35) as 


dl v 

dr v 


— Iu ■+■ 1 Si/ . 


(4.36) 


The general solution to this equation is 

Iu{r v ) = MO) e”* + r e-^'^ S u {t' u ) dr' v . (4.37) 

Jo 

We can interpret the first term on the right hand side as the initial intensity dimin- 
ished by the absorption. The second term is then the integrated source diminished 
by absorption. 

For a constant source function, the solution for (4.37) is 


Iv — Sp + e T>/ (MO) — Sy) • 


(4.38) 


We see that the intensity approaches the source function as the optical depth in- 


creases. 
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Rybicki and Lightman (1979) have offered a simple physical interpretation 
of the source function and the transfer equation. If we consider (4.36), then for 
I„ < S v , /„ will increase along the ray. If I v > S„, then /„ will decrease along the 
ray. In either case, I v will tend towards S v . The intensity will actually become 
equal to the source function if the optical depth is sufficiently large. 

From Kirchoff’s law for thermal emission, the source function is dependent 
only upon the temperature of the medium, and is given by 

S„ = Bv(T) s ^ \e~ hv l kT - lp 1 , (4.39) 

c 

where h is Planck’s constant, c is the speed of light in vacuum, k is Boltzmann’s 
constant, and T is the temperature of the medium. From (4.38), we see that for 
sufficiently large optical depths, the intensity will also approach B V {T). This is the 
condition for blackbody radiation. 

It is often more convenient to rewrite the equation of transfer in terms of the 
Einstein coefficients. Let us first consider a two level system that is in thermody- 
namic equilibrium at temperature T. If ni and ri 2 represent the number density of 
atoms in levels 1 and 2, respectively, then the steady state solution dictates that 

ni B\ 2 J = ri2 A 21 + n 2 i?21 J • (4.40) 

Since the line profile is usually very narrow, we may approximate the integral in 
(4.30) as being equal to J. Solving for J, we find 



In thermodynamic equilibrium, the level population ratio is given by 

!li = il e ^2i/ T , (4.42) 

*2 92 

where hv2i is the energy difference between the two levels and g x is the degeneracy 
of the i th level. However, in thermodynamic equilibrium, J u = B v , and hence we 
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see that 


giBu = giBn 

2h i/|i . 

A21 = — j — -® 2 l ( 4 - 43 ) 

c 


From the definition of the emission coefficient, we can see that in terms of 
the spontaneous emission coefficient, 

3v = ~~^ n 2 ^21 4>{ u ) • (4.44) 

4 7T 

Likewise for the absorption coefficient, we may write 


hi / 21 

a„ = — — 4 >{u)(niBi 2 - n 2 B 2 1) 

47T 


(4.45) 


It should be remembered that the stimulated emission is proportional to the inten- 
sity and is, in many respects, like negative absorption. For this reason, it is included 
in the absorption coefficient. Using (4.29), (4.44), and (4.45), we obtain 


dJ v 
da 


—~7~{ n \Bi2 - n 2 B 2 i) H”) *v + ~r~ n 2 ^21 H") 
47T 47 T 


(4.46) 


We can rewrite this using the optical depth and the source function. The absorption 
coefficient and the source function may then be rewritten as 

<*V = ni B12 (1 - 

4ir \ g 2 n\) 


Su — 


T 12 A 21 


ni5i2 - 112 B21 
= 2 hi/ 3 

'(as- 1 ) -1 - 


(4.47) 
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4.5.3 Collisional Excitation and Multi-Level Systems 

Most syst ems do not have just two levels, but usually many levels that must 
be considered. Collisional excitation and. de-excitation processes also add to deter- 
mine the final equilibrium level population of atomic and molecular systems. Let 
be the collision rate coefficient for producing a transition from level j to level k. 
The units are cm 3 s -1 . If n c is the number density of the colliders, then n c 7jk is the 
probability rate for a collisional transition. Let us consider the level. In a steady 
state, the rate of change of nj is zero. That is to say, the processes which contribute 
to increasing the j th level exactly match those which contribute to decrease it. This 
may be written as follows: 

Y "k(”c7kj + Skj-^) + Y n k^kj = n j EH + £jk«A*j) + Y 
k#j k>j *-k^j k<j 

The right hand side of (4.48) represents the processes which will increase the j** 1 
level population. They are collisional excitation and de-excitation into the level, 
stimulated emission and absorption, and spontaneous emission from higher levels. 
The left hand side represents the debit half of the ledger, and includes the same 
processes; only this time depopulating the level. 

The collisional rate coefficients 7^ and are related to each other through 
detailed balance. This is simply derived, and results in 


. (4.48) 


7jk=7kj^e-^-^ 

J <7j 


(4.49) 


where 2?j is the energy corresponding to the j th level. 

It is possible to solve (4.48) for a system that has reached a steady state 
if all of the collisional and Einstein coefficients are known. The level populations 
will then be known. Use of (4.46) and (4.47) will then enable one to solve for the 
intensities and fluxes that one would expect from the medium. In general, this is a 
non-trivial process if the number of levels is greater them 2. 
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4.6 Escape Probability Theory 
4.6.1 Castor’s Method 

For regions that contain velocity gradients, the radiative transfer is compli- 
cated still further owing to doppler shifts that vary throughout the medium. If the 
flow velocity is much greater than the thermal width, as in the case of OMC-1, then 
it seems clear that a photon’s chance of being absorbed diminishes quite rapidly 
the further it moves away from its source. Hence, the radiative transfer becomes 
a local problem and it is possible to treat the chance of a photon travelling out of 
the region as an escape probability. Much work has been done on this subject, first 
introduced by Sobolev (1957). We shall mainly follow the relevant papers of Castor 
(1970) and Goldreich and Kwan (1974). 

Since the escape probability formalism was first established in order to under- 
stand spectral line formation and profiles as a result of radial flow in the atmospheres 
of Wolf-Rayet stars, we shall derive the relevant formalism assuming spherical ge- 
ometry. The minor revisions for plane-parallel geometry will be discussed later. 
One should keep in mind that we will be dealing with radiative line transfer and 
not continuum transfer. 

We begin by defining the geometry. We assume that we have an optically 
thick internal source of radiation which constitutes the core of the region and has a 
radius r c . Owing to the spherical geometry, we need only consider radiation moving 
in one specific direction; viz., the z-axis. This radiation will be a function of the 
cylindrical coordinate, p. The geometry is shown in figure 4.5. It should be clear 
that p 2 + z 2 = r 2 . We can then write the equation of transfer, (4.29), as 

= a(u,p,z)(^I(u,p,z) - S(r)^ . (4.50) 

The sign convention is a little different in (4.50) from that in (4.29), due to the way 
in which the coordinate systems are set up. 

There exist two regimes for solutions. The first one is for p > r c or z > 0. 
If p > r c , then the core will not contribute directly to the observed intensity if 
scattering is negligible. The second regime must include the contributions from the 



FIGURE 4.5 A schematic of the geometry used for the escape probability formal- 
ism. P and z define the coordinates, and r c delineates the radius of the optically 
thick continuum core. After Castor (1970). 
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core. In the first regime, we find using (4.37) that 


I(v,p,z) = j S ^(p 2 + z ,2 )l^ exp[r(»/, p, z) — r(v,p, z')\ dr(v, p, z') . (4.51) 


If p < r c and z < 0, then the solution is 




/(^P,*) = J s([p 2 +z , 2 )l')exp\T(v,p,z)-T(l/,p,z , )]dT(l/ t p,z , ) 

+ / c exp{-[r(i/,p,-(r? -p 2 )*) -r(v,p,z)}} , (4 


(4.52) 


where J c is the intensity from the optically thick core. The optical depth for a ray 
traveling at constant p is defined as 


'(i'.P.s) = / a(v y p,z')dz' . 
J-OO 


(4.53) 


Since we are assuming spherical symmetry, the velocity of a region is dependent 
only upon the distance from the center. Hence, for a line profile 0(Ai/) in the rest 
frame, and a rest frame absorption coefficient at line center equal to ao(r), we have 


= a 0 (r) <t>(^+ ’ 


(4.54) 


where vq is the frequency of the radiation at line center. Since <f>(Av) is usually a 
very narrow function, most of the contributions to a(i/, p, z) will come from a rather 
small region where A V « 0. This occurs for z « zq, where zq satisfies 


r ((P 2 + *§)*) 
c(p 2 + zJ)J 


vq-v 


(4.55) 


We can then approximate the optical depth equation, (4.53), by 

r(v,p,z) ss ao((p 2 + *o)*) J ~ *&) dz> • (4.56) 
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Substituting x for the argument of <f>, we cam define the function y to be 

y(x) = [ 4>{x') dx' , (4.57) 

J—OO 

with the boundary conditions y(— oo) = 0 amd y(oo) = 1. Substituting into equation 
(4.56), we find 


r(v,p,z) = 


ao 


((p 2 + *o)*) 


(*UWL 




+ *m z 

c r 


- V 0 ) = r(t/, p, oo) y . (4.58) 


Evaluating the partial derivative in (4.58) amd dropping the subscript 0, we obtain 
for the optical depth at infinity 


r(t/,p,oo) = ac(r) 


^ ( ,) 1+ £(ita£d ^ 


(4.59) 


Knowing the optical depth, we cam now substitute equation (4.58) into (4.51) amd 
(4.52) to find the intensity. Due to the fact that the major contribution to the 
intensity will have its source in a very small range of z, we may approximate the 
integration to find 


/(i/,p,2) « S^(p 2 +Zq)^) 1 - ex P{ 1 ‘(^,P,oo) [y(^) ~ 1]} , (4.60) 


for p > r c or z > 0. For p < r c , z < 0, the solution is 


I{u,p,z) m s((p 2 +4) *) jl — exp jr(i/,p,oo) \y{z) - y(-(r 2 -p 2 )s) } j 

+ / c exp{r(i/,p,oo) jy(z) - y(-(r 2 - p 2 )*)] } . (4.61) 


In order to evaluate (4.48), the detailed balamce equation, we need to deter- 
mine J, the meam intensity. This is evaluated in cylindrical coordinates by 

J = ~ J dz j 4>(u + — ¥-p-z - uq^J l(v, (r 2 - z 2 )5,z) du . (4.62) 
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Equations (4.60) and (4.61) can then be substituted into (4.62). Since most of the 
contributions to the integrals come from values of u and z that satisfy z ss zq , we 
can expand the source function as well as (p 2 + Zq J 1 / 2 in a power series. The part 
of (4.62) that involves the source function then reduces to 


hJ'A s(r) C{ l-exp[r(i/,p,oo)(y-l)]}^(Ai/)di/| . (4.63) 


We need to digress here for just a moment and consider the form of equation 
(4.38). This equation can be rewritten as 


/„ = e~ Tv J„( 0) + (1 - e~ rv ) S„ . (4.64) 

Looking at the contribution from 7(0), we see that the coefficient e~ Tv represents 
the probability that the background photons incident on the medium will escape 
the medium. 

We may apply a similar understanding for the contribution from S (r) to J 
in the Wolf-Rayet star model. If we denote the escape probability by /3, we may 
write (4.63) as (Elitzur 1984; Rybicki 1984; Canfield, Puetter, and Ricchiazzi 1981; 
Castor 1970) 

J = (1 - 0)S{r) . (4.65) 


From this, we see that 
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(4.66) 


where \i is the cosine of the angle between the ray and the radius, or equivalently, 
z/r; and 
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where 

mM = ' <4 ' 68) 

Likewise, with the contribution from the core, we cam find the escape probability 
for that radiation and this is (Castor 1970) 

1 yl 1 _ e-'oK**) 

p c = — / dp 

2 ytl-rf/r 2 ) 1 /* WM) 

It is then very easy to find the total J, which is given by 



J = (l- 0)S(r)+0 c I c 


(4.70) 


For plane-parallel atmospheres, the results are similar to those derived above 
for spherical geometries. Scoville and Solomon (1974) have calculated the escape 
probability for such a geometry using the method of Lucy (1971), who used the 
Eddington approximation (Stacey 1985). The result is that for the two-stream 
approximation, the emergent flux requires n = l/\/3. This yields the following 
expression for the escape probability: 


0 = 


1 - e -3r 0 (z) 
3r 0 (s) 


(4.71) 


where 


t 0 (z) = a 0 («) 


c dz 


v 0 dV{z) 

Equation (4.70) simplifies somewhat to become 


(4.72) 


J = {l-0)S{z) + -0I o +0I BB , 


(4.73) 


where S(z) is given by (4.47), Iq is the background radiation field impinging on 
only one side of the cloud, and I BB is the 2.7 K microwave background radiation 
field (Elitzur 1984). 
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4.6.2 Detailed Balance and the Escape Probability Method 


The goal of the escape probability method is for us to be able to determine in 
a self-consistent fashion the level populations for the atomic or molecular systems 
under consideration as well as the line intensities or fluxes that we should observe. 
To this end, we shall construct once again the equation of detailed balance. In 
keeping with convention, we shall redefine nj to be the population density of the j th 
level divided by the degeneracy. 

We begin by looking at the level and determining the rate of change of the 
population due to the various processes that occur. For now, let us be concerned 
only with the levels that lie above j. Spontaneous emission from these upper levels 
will tend to populate the level. Since the radiative transfer of the line photons is 
now a local problem, we see that a spontaneous emission will result in an actual 
increase in the population of level j if the photon can escape. That is to say, 


dnj\ 

^ / spoilt 



(4.74) 


Let the background radiation field be the cosmic blackbody field. Stimulated 
emission from level k to j will also help to populate level j. Similar arguments 
concerning the true amount of population due to this process lead us to the result 
that 

(iff) ’ ( 4 - 75 ) 

\ ai / stim k> j x 9} ' 

where F(T BB^j) the photon occupation number, equal to [e^/^BB _ l] -1 , 

for a blackbody temperature Tbb = 2.7 K. For absorption from level j to k, we can 
write in the same spirit 

(^)abs = ~ ^ kj Akj nj 

k>j 


Adding (4.74), (4.75), and (4.76) together, we find the radiative processes involving 
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the higher levels will contribute to a true population rate for level j of 

^ = XJ(— ) ^kj ^kj [ n k ( n k — n j) •^ , (^BB»* / kj)] ■ (4-77) 


Collisional processes will also help to either populate or depopulate the levels. 
For the higher levels, we have either population of level j through de-excitation or 
depopulation of level j through excitation. Using (4.49), we can write the collisional 
term as 



= E & C kj 
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n k ~ 


rij e At ^j/* r kin t 


(4.78) 


where C k j is the collisional transition rate from level k to level j (= n c ) divided 
by the degeneracy of level j, and T ^ is the kinetic temperature of the region. We 
can determine the corresponding rates for the lower lying levels, and the final total 
rate is given by 



^ [(“•) ^kj ^kj n k + ( n k ~ n j) ^(^BB.^j)] 
k>j L 9 > 


+ 0kQci 


C kj (n k -nje Ki/Wkin) 


^jk A jk [ n j + ( n j ~ n k) ^(^BB.^jk) 

k<j L 


+ 9k Cfc («j - ^k e h ^ /kT ^ 


(4.79) 


Equation (4.79) will give us m equations with m unknowns for an m-level 
system. We solve for nj, and from that, we can calculate the emergent intensity. 
Using (4.38) and (4.47), we find that 
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123 

Filially, it is also possible to include the role of internal sources of radiation, 
such as dust grains imbedded in the cloud. See Elitzur, Goldreich, and Scoville 
(1976) for details. 



V. Shocked OH and CO Emission Lines 
From The Orion Molecular Cloud 

5.1 Introduction 

Regions of recent star formation evince massive bipolar flows, which are ap- 
parently a common phenomenon (Bally and Lada 1983) . These flows show velocity 
widths in the low J CO emission line wings of up to 130 km s —1 , and should be 
hypersonic. One then expects shock waves as a result of the flows plowing into the 
ambient molecular cloud surrounding the star formation region or else adjacent to 
it, as in the blister models of Balick, Gammon, and Doherty (1974) and Zuckerman 
(1973) . The presence of vibrationally excited H2 in almost all of these regions 
suggests that there is shock heating occuring as well (Kwan 1977). 

Rotational and vibrational emission lines of H2, CO, and OH serve as excel- 
lent tracers of different parts of these shocked regions (Chemoff, Hollenbach, and 
McKee 1982), as do other molecular constituents, such as NH3, CS, and SO (Keene, 
Blake, and Phillips 1983; Townes et al. 1983; Genzel et al. 1982; Plambeck et al. 
1982; Goldsmith et al. 1980; Ho and Barrett 1978). 

For typical preshock densities and flow velocities, the shocks should induce 
. CO emission from a broad range of rotational levels from J = 1 to J « 35 (Draine 
and Roberge 1984; Chernoff, Hollenbach, and McKee 1982; Draine and Roberge 
1982; McKee et al. 1982; and references therein). These levels range from several 
degrees Kelvin to several thousand degrees Kelvin (see figure 4.2), and therefore 
sample different regions within the molecular cloud. In fact, the low J transition 
emission features of CO usually display a very narrow spike from the ambient qui- 
escent gas and a very broad plateau from the shocked gas. 

One of the most extensively studied star formation regions is the Orion Molec- 
ular Cloud (OMC-1). During the last 20 years, OMC-1 has been the subject of 
observing programs at a variety of wavelengths from the near infrared to the radio. 
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Of interest to us is the region of the cloud that contains the Becklin-Neugebauer 
object as well as the Kleinmann-Low nebula. This region appears to be a site of 
recent star formation which has triggered a massive outflow of material. This has, 
in turn, contributed to the production of H 2 O, OH, and SiO masers, which can be 
used to outline the shock front (Plambeck et al. 1982; Genzel et al. 1981). 

Since the OMC-1 is relatively nearby, being at a distance of 500-530 parsecs 
(Mezger and Hoglund 1967; Becker and Fenkart 1963), we are given a unique chance 
to study shock chemistry and physics through the observation of molecular rota- 
tional transitions in the FIR. By comparing the intensities of different rotational 
lines involving the same molecule, it is possible to derive a set of constraints on the 
actual temperature and density of the emitting region. 

The CO lines that have been measured thus far suggest a two-component 
model for the shocked and post-shock regions. The warm component has a temper- 
ature on the order of 750 K and a molecular hydrogen density of about 3 x 10® cm -3 . 
For the hot component, these values are 2000 K and 1 x 10® cm -3 , respectively (Wat- 
son 1982). 

This thesis is concerned mainly with the initial detections of several FIR OH 
rotational lines in emission from the BN-KL shock region in OMC-1. The hydroxyl 
radical is much more readily dissociated than CO, and is therefore much more 
sensitive to the physical conditions in this region. Because of this, we do not expect 
the OH to be as extensive as the CO, but to be found in significant abundance in 
a relatively thin sheet within the shock front (Draine and Roberge 1982). 

Along these lines, we have also made the first astronomical detection of the 
J = 15 — ► 14 rotational transition of CO at 173 ^m. The significance of this line 
shall also be discussed below. 

5.2 A Brief Overview of the Orion Nebula 

The Orion Nebula is barely visible to the naked eye under the best conditions 
as an unresolved patch of light. It is not mentioned in any records, medieval or 
ancient, and it is interesting to note that Galileo also apparently never pointed his 
telescope in that direction (Gingerich 1982; Burnham 1977). Nicholas Peiresc is 
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credited with discovering the nebula in 1611. It wasn’t until 1656 that a drawing 
of the nebula was published by Christiaan Huygens. It is claimed that Orion so 
stirred Sir William Herschel when he first looked toward the heavens with his new 
telescope in 1774, that he dedicated the rest of his life to such studies. It was in 
1771 that Charles Messier published his first listing of clusters and nebulae, and 
the designation M42 is probably the most well known catalog name of this object. 
In 1880, Henry Draper took the first successful photograph of a nebula, the subject 
being Orion. 

Jumping ahead 100 years, we find a consistent picture of the Nebula begin- 
ning to emerge from the data collected thus far. Visually, the nebula is illuminated 
by the Trapezium stars and 6 2 Orionis. These stars are relatively young, and there 
are stars still in the process of formation within the region. The nebula is composed 
of the surrounding Hn region, ionized by Orionis, which forms a blister on a 
much larger giant molecular cloud (Werner 1982; Pankonin, Walmsley, and Harwit 
1979; Balick, Gammon, and Doherty 1974; Zuckerman 1973). This molecular cloud 
has an estimated mass of 2.5 x 10 4 -1.0 x 10 5 M© (Tucker, Kutner, and Thaddeus 
1973). 

To the north and west of 6 1 C Orionis lies the infrared nebula known as 
the Kleinmann-Low (KL) nebula (Kleinmann and Low 1967). This near-infrared 
source of radiation is most likely the recent birthplace for several hot O and B 
stars. This nebula was discovered after a compact 2 nm source in the vicinity 
was found by Becklin and Neugebauer (1967). These two objects are about 11 7/ 
apart on the sky, and are collectively known as BN-KL. The coordinates of BN are 
a 1950 = 5 h 32 m 46.7 B , 6 1950 = -5°24'l7". 

In the region of BN-KL, there is an infrared source known as IRc2. This 
source is now believed to be the central engine that is supplying energy to the 
massive outflow that has been observed. The dynamic processes were first found 
to be phenomena associated with outflows centered on BN or IRc2 through milli- 
metric CO observations by Solomon, Huguenin, and Scoville (1981). This outflow 
was determined to be molecular and bipolar by Erickson et al. (1982). Directly 
surrounding IRc2 is an expanding central disk that is about 70 astronomical units 
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in size (see figure 5.1). This disk is inferred from interferometric observations of the 
SiO maser (Wright and Plambeck 1983) as well as from the placement and velocity 
centroid of the so called H 2 O “shell” masers (Genzel et al. 1981). 

Surrounding this inner disk is a torus or doughnut of gas that is expanding 
at low velocity (~18 km s -1 ) and is most apparent in thermal emission from SO 
(Plambeck et al. 1982) and in the low velocity H 2 O and OH masers (Downes et al. 
1981). This flow is in the NE-SW direction, very much along a ridge of quiescent 
gas. The torus is somewhat tilted with respect to our line and on the order of 0.02 
pc in size. 

Perpendicular to the ridge, and nearly centered on IRc2, hot high velocity 
gas has been found to exist on a scale of roughly 0.1 pc. It was initially discovered 
by Thaddeus et al. (1972) in an H 2 S line. Observations of SO, SO 2 , HCN, and 
SiO also showed evidence for high velocity wings (Dickinson et al. 1976; Snyder et 
al. 1975; Wannier et al. 1974; Gottlieb and Ball 1973). Earlier CO observations, 
showed the high velocity gas, but it was overlooked due to the intense, low-velocity 
central peak. This gas is believed to be part of a high velocity outflow in the NW- 
SE direction centered on IRc2. Looking at the H 2 2/im emission map of Beckwith 
et al. (1978), we see that the emission forms a double-lobed structure centered 
on IRc2 (see figure 5.2). This structure is perpendicular to the quiescent ridge 
and is on the order of 0.1 pc in extent. The line profiles for H 2 , CO, OH, and 
other molecules show high velocity plateau features with widths ranging from 50 
to 130 km s _1 . The H 2 data suggest temperatures of 2000 K and densities on the 
order of 10 5 cm -3 . The velocities, temperatures, and densities make sense if we 
assume am hypersonic outflow into the ambient molecular cloud, forming a shock 
wave which heats and excites the various molecules. 

The BN object was at first considered to be a protostar, but it is now consid- 
ered to be a young hot star that is losing mass (Downes et al. 1981). This mass loss 
has enshrouded the star in a cloud of gas and dust. Recent radio observations at 
the VLA indicates a size of no more than 0 , (07, or 5 x 10 14 cm (Moran et al. 1983). 

Continuum maps in the FIR show a single maximum centered on the BN-KL 
region (Keene, Hildebrand, and Whitcomb 1982; Werner et al. 1976). The source 
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FIGURE 5.2 Locations of H 2 O masers in BN-KL. The LSR radial velocities are 
given next to the maser locations. Solid contours are from the 20 nm. continuum 
data. Dashed contours are those of the t; = 1 -+ 0) S(l) transition of H 2 from 
Beckwith et al. (1978). From Downes et al. (1981). 
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of this continuum is most likely dust grains that absorb the starlight from IRc2, 
BN, and KL, and then reradiate it, mostly from the infrared to the submillimeter 
part of the spectrum. Since the source size (HWFP) is roughly 1* for the continuum 
source and for the shocked gas, and this size corresponds to the diffraction limited 
beam size in our observations, we look at the BN-KL object to observe the shocked 
gas. 

For other histories and overviews on the Orion Nebula, see Stacey (1985); 
Glassgold, Huggins, and Schucking (1982), and Watson (1982). 

5.3 OH Observations Toward BN— KL 

5.3.1 Inspiration 

The BN-KL region has been the subject of many recent millimeter and sub- 
millimeter observations in the rotational transitions of OH and CO, as well as in 
a plethora of transitions of diatomic and more complex molecules (Rydbeck and 
Hjalmarson 1985; Schultz et al. 1985; Johansson et al. 1984; Storey et al. 1981; and 
references therein, to name a few). The CO and OH measurements suggest that 
the emitting region resides in the shocked and post-shock gas. Measurements of the 
various ro-vibrational transitions of H 2 have shown that njj 2 ~ 2 x 10 s cm -3 within 
the shock region and that the temperatures are on the order of 2000 K (Nadeau 
and Geballe 1979; Beckwith et al. 1978). The Berkeley group found that the CO 
transitions are consistent with a temperature ~750K and molecular densities ~3 x 
10 6 cm -3 for the warm component of the emitting region (Watson 1982). This 
suggests that most of the observed CO lies within the post-shock gas closest to the 
shock front. 

Since OH is chemically more active than CO, we do not expect the OH 
to coexist extensively with the CO. More likely, the OH will be abundant in a 
relatively confined sheet in the shocked or post-shock gas, with smaller abundances 
throughout the rest of the region. The OH may also tend to be more clumpy than 
the CO. 

Recent observations of the 2 n 3 / 2 J = j — ' ► | transition at 119 /*m (Watson 
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1982; Storey, Watson, and Townes 1981), along with several A-doubling radio lines, 
led us to observe four more OH transitions in the FIR. From this set of observations, 
it should be possible to apply radiative transfer theory to place constraints upon 
the temperature and density of the OH-emitting region. 

The FIR observations of OH are important not only in understanding the 
physical parameters of the OH-emitting region, but also in attempting to understand 
the pumping mechanism of the OH masers that are associated with BN-KL and 
other H II regions. Much has been written about the pump mechanism, and it is still 
not known if radiative pumping or collisional pumping dominates (see, for example, 
Genzel 1985; Elitzur 1982; Reid and Moran 1981; and references therein) . By 
observing the two parity transitions for each FIR transition, we may be able to 
shed some light on this issue. 

5.3.2 The 163 pm Observations 

In 1983 Much, we made sub-millimeter OH observations toward the BN- 
KL region. The 91cm bent-Cassegrain telescope on board the KAO was used 
as our observing platform. The observing altitude was 13-14 kilometers and the 
diffraction-limited beam size was l(0 square. The instrument used was our interfer- 
ometer/grating-spectrometer, capable of an instrumental resolution A/AA ~ 1200, 
and described in detail in chapter II. The Orion data from this flight series consti- 
tuted the first detection of the 2 II 1/ /2 *7 = § — * j OH rotational transitions and was 
reported in Viscuso et al. (1985a). The laboratory wavelengths for these transitions 
are 163.12 and 163.40 ^m (Brown et al. 1982). In 1984 January and 1985 Jan- 
uary, we reconfirmed the 163 nm flux measurements with an improved instrumental 
resolution of «4000. 

As we saw in chapter II, the spectrum one derives from an interferogram is a 
convolution of the instrumental profile of the spectrometer with the source profile. 
The instrumental profile, as best as we can tell from laboratory measurements, is 
approximately Gaussian. 

The OH results from the 1983 March flight series are displayed in figure 
5.3. The instrumental resolution at the time was roughly 1200 at 163 /mi. The 
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FIGURE 5.3 Spectrum of BN-KL at 163 /im. The smooth dashed curve represents the expected continuum profile. 
The CO line is due to the J = 16 — ► 15 rotational transition. This line is much stronger than the two OH lines from the 
2 n i/2 J = j — » j transition. Regions of telluric water vapor contamination are marked. The noise may be estimated 
from the departure of the continuum data from the expected profile. Data were obtained in 1983 March. 
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interferometer was not upgraded until the next year. At this resolution, it is possible 
to separate the two sets of lines arising between the different parity states, but the 
hyperfine transitions cannot be separated. The hyperfine transitions within one 
parity transition occur at wavelengths that differ by one part in 10 5 , on average. The 
two parity transitions (see figure 4.3) are separated by 0.275 nm , and are therefore 
resolvable at this instrumental resolution. 

Figure 5.3 displays the 1983 March spectrum which is drawn on top of the 
instrumental profile (dotted line). Since the continuum level is roughly constant 
over this small wavelength range, the dotted line represents the instrumental profile 
times a scaling factor that simply reflects the constant continuum level. 

Werner et al. (1976) have found the continuum intensity from the BN-KL 
region to follow a blackbody distribution for a temperature of 70 K and with an 
emissivity that is inversely proportional to wavelength and is unity at 20 pm. At 
163 nm, for a l' beam, we find that the continuum flux level we measure is 5.77 x 
10“ 17 W cm -2 in each bin at a resolution of 1200. 

Looking at figure 5.3, the first obvious feature is the J = 16 — ► 15 rotational 
transition of CO at 162.8116 /xm. This line was first observed by Stacey et al. (1983), 
and the flux in the line is 7.8 x 10 -17 W cm -2 . Stacey et al. (1983) give a slightly 
smaller estimate, but the two agree within the errors, most of which are due to 
determining where the continuum level actually lies on the spectrum. 

The next two bumps to the right of the CO line are the two OH lines. 
Although the lines are not necessarily significant (i.e., the signal-to-noise is not 
greater than 5), they are suggestive since they do happen to lie in the correct 
positions and are both raised the same amount. If we also look at the fit to the 
continuum, we see that the fit is very good except for those two lines and the two 
places where telluric water vapor features lie. There is very little signal in the wings, 
which indicates that the overall system noise is low. Another indication that the OH 
lines were real came from the removal of the CO line from the original data. Since 
the wavelength of the CO line is well known, its contribution to the interferogram 
can be calculated and removed. This removal preserved the OH features, indicating 
that they were not caused by the dominant CO line. 
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The interferometer was upgraded in 1983 December, and was successfully 
flown the next month on the KAO. The 163 Atm OH line measurements were con- 
firmed at an instrumental resolution of » 3700, albeit not with much improvement 
in the signal-to-noise ratio. In 1985 January, we again looked at BN-KL and dedi- 
cated another 30 minutes of observing time to obtain five 163 Atm spectra at mod- 
erately high resolution (R«2800). The result is shown in figure 5.4. The data 
superimposed upon a Gaussian instrument profile which should define to reason- 
able accuracy the continuum level from Orion. The flux levels are 1.75 x 10 -17 and 
2.1 x 10~ 17 W cm -2 for the 163.12 and 163.40 /im lines, respectively. These values 
are higher than the 1.2 x 10~ 17 W cm -2 derived from the 1983 March data (Viscuso 
et al. 1985a). 

These data show a poor fit to a Gaussian at the lower and upper ends of the 
wavelength scale. Here, the signal seems to be much higher, especially on the short 
wavelength side. The reason for this is the fact that the spectrometer is not aliasing 
properly for such a small free spectral range at this wavelength. This can be seen in 
figure 5.5, which shows the data superimposed on a laboratory calibration spectrum 
run near this wavelength. The calibration clearly shows significant asymmetry, with 
the lower end of the spectrum being effectively raised. The free spectral ranges for 
the calibration spectrum and the data are slightly different, and that most likely 
explains why the match is not perfect. However, the calibration spectrum does 
show the effect of imperfect aliasing, and explains the raised ends of the spectrum. 

Subsequent confirmation attesting to the existence of the 163 Atm lines has 
been offered by Melnick (1985), who used the Berkeley group’s Fabry-Perot spec- 
trometer to detect the 163 /xm OH lines. Their instrumental resolution at the 
time was 9000 at this wavelength. Flux measurements were only secured for the 
163.12 ^m line, giving 1.3 x 10 -17 W cm -2 (Crawford et al. 1986). 

5.3.3 The 84 /xm Observations 

In a cooperative venture with Dr. Edwin Erickson and his group at NASA/ 
Ames Research Center, measurements of the ^n 3 / 2 J = j — 1 ► § transition of OH at 
84 /xm were acquired on 1983 November 22-23 using the NAS A/ ARC Cooled Grat- 
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FIGURE 5 . 4 Spectrum of BN-KL at 163 pm. The smooth solid line represents the 
expected continuum profile. The interferometer/spectrometer was used to obtain 
these data in 1985 January with an improvement in the system resolution by a factor 
greater than 3 over that used to obtain the data in figure 5.3. The integration time 
was about 50 seconds per point. The signal scale is arbitrary. The rise in the profile 
on the left is due to an incomplete aliasing of the order by the instrument. The 
CO line does not appear as strong in these data due to its placement far out on the 
instrumental wing. 
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FIGURE 5.5 Spectrum of BN-KL at 163 /tm. These data are the same as those 
in figure 5.4, except that the instrumental profile was replaced with a laboratory 
blackbody spectrum. Note the better profile matchup on the left. The sampling 
rates used for both spectra are different, and that accounts for the imperfect match 
between the two profiles. 
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ing Spectrometer (CGS) facility instrument (Erickson et al. 1984). This instrument 
basically consists of a huge grating that is operated in high order (m ~ 75), and has 
a resolution of from 4000 to 10000, depending on wavelength. The whole assembly 
fits in a cryogenic barrel that looks like a somewhat scaled down version of a 55- 
gallon drum. At the time of operation, there was an array of 6 Ge:Ga detectors 
which were separated by 0.0184 nm. The diffraction overlap caused the resolution 
to be 26% lower than that found by simply using the detector separation. The 
resolution was 3400 and the circular beam size was l' in diameter. 

Figure 5.6 shows the spectra for the two parity transitions at 84.42 and 
84.60 pm (Brown et al. 1982). Since the detectors needed to be calibrated, the raw 
spectra were divided by a continuum spectrum at a nearby wavelength of 87.08 /im, 
a region free from any expected telluric or astronomical spectral features. Owing to 
different instrumental response at the different wavelengths, the baselines are sloped. 
Detector #5 was inoperable during the flight, which accounts for the missing point 
in the data. Line flux determination is again made by comparing the line strength 
against the continuum level. For the 84.42 pm line, the flux is 1.0 ± 0.3 x 10 -17 W 
cm -2 , while the 84.60 ^m line flux is 1.4 ±0.4 x 10“ 17 W cm -2 . These results were 
first reported in Viscuso et al. (1985b). Verification of these lines was published 
later in Watson et al. (1985). It should be noted, however, that the latter group 
had actually made their observations earlier. Their beam size is 30" compared with 
1 ; for ours, and the resulting flux level is down by roughly a factor of 2. 

The 84.4202 pm data are controversial at this point because of the nearby CO 
line at 84.4107 ^m (/ = 31 — ► 30). The two lines are separated by only 33 km s -1 , 
which is not enough for them to be resolved with the Berkeley FPS or the NASA/ 
ARC CGS. The measured flux from the line is consistent with that expected from 
the CO, given the data from all of the other CO lines that have been unambiguously 
measured. However, the flux is also consistent with that expected from OH given 
the flux from the other OH line at 84.60 /xm. More will be said about this later. 
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FIGURE 5.6 Spectrum of BN-KL at 84 Data were obtained in 1983 Novem- 
ber with the NASA/ARC facility cooled grating spectrometer. Dashed lines rep- 
resent the expected instrumental profiles for the emission lines. The solid lines 
represent the continuum baselines. The residual baseline slopes are due to instru- 
mental effects. 
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5.3.4 Other OH Observations 

Other various rotational transitions in OH have been observed toward BN- 
KL. Table 5.1 gives a list of observations that have been carried out and reported, 
including intensities and beam sizes where known. 

The A-doubling lines of the 2 n 3 / 2 J = § level have been well known for 
many years. Maser action is seen for the 1665 MHz line towards BN-KL (Reid 
and Moran 1981; and references therein). These masers are associated with the 
expanding doughnut of gas, and their formation is still not completely understood. 

The 2 n 3/2 J = \ F = 4 + — ► 4~ A-doubling transition (13441MHz) was 
detected in emission from BN-KL by Matthews et al. (1984). This feature is 
about 25km s -1 wide (23 ± 4), with a line antenna temperature of 20 ± 4mK. 
The F = 3 + — ► 3~ hyperfine transition at 13435 MHz was marginally detected with 
a line antenna temperature of 10 ± 5 mK. The 140 foot telescope of the NR AO was 
used for the observations, with a sensitivity of 4.5 Jy K -1 . 

The 2 n 3/2 J = | F = 4~ — ► 4 + and F = 5“ -» 5 + A-doubling transitions 
at 23817 and 23826MHz have not been detected toward Orion (Walmsley 1985). 
Guilloteau et al. (1984) report a possible detection of the 6035 MHz line due to the 
2 n 3/2 J = \ F = 3 — 3 + A-doubling transition. Knowles, Caswell, and Goss 
(1976) have published a detection of this line toward Orion, showing a relatively 
narrow emission with a flux density of 4 Jy. This may possibly be a weak maser 
line. 

The 2 fI 3/2 J — § — > > | rotational transitions of OH at 119.23 and 119.44 ^m 
have been detected by the Berkeley group (Watson et al. 1985; Watson 1982; Storey, 
Watson, and Townes 1981). The more recent estimates for the line fluxes are on 
the order of 7 x 10 -18 W cm -2 . 
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Table 5.1 OH Observations of OMC-1 

Transition 

Wavelength 

Flux 

(10 -17 W cm -2 ) 

Reference 

2 n x /2 J = | -*■ i 

163.40 Atm 

1.2±0.6 

1 


163.40 Atm 

2.1±0.7 

2 


163.12 Atm 

1.2±0.6 

1 


163.12 Atm 

1.7±0.7 

2 

2 n 3 / 2 j = I -» § 

119.44 Atm 

0.7±0.3 

3 


119.23 Atm 

0.83±0.3 

3 

2 ir 3 /2 j = J -*■ | 

84.60 Atm 

1.4±0.4 

4 


84.42 Atm 

1.0±0.3 

4 

2 n 3/2 j = ]f = 4 + 

— *• 4“ 2.24 cm 

9 X 10 -9 

5 


A list of the published OH observations to date toward the Orion Molecular Cloud. 
The errors in the 163 and 84 Atm lines are primarily due to uncertainty in the 
calibration as opposed to variances in the signal strengths. The reference numbers 
refer to the following articles: 1) Viscuso et al. (1985a); 2) this work; 3) Watson et 
al. (1985); 4) Viscuso et al. (1985b); 5) Matthews et al. (1984), 
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5.4 OH in Orion: A Simple Approach 
to a Complex Problem 

5.4.1 The Modus Operandi 

We now have a limited, but significant, set of observations involving the lower 
level transitions of OH found in the shock region of OMC-1. It is the purpose of 
this chapter to see if these observations are consistent with each other. 

For the most part, we will assume that there is only one overall area of the 
cloud that contains the OH emitting regions, an assumption that is not necessarily 
true. Another simplifying assumption that is invoked is that of uniformity among 
the emitting regions. Even though we allow for clumps when we introduce the 
filling factor, we assume that the clumps may be explained with one temperature 
and density. Therefore, we will assume that there is an average temperature and 
density which will represent the OH emitting regions. 

We will begin by analyzing the 163 /xm data, deriving from them a molecular 
hydrogen density. Then, we can determine the expected 119 /xm flux, using the new 
collisional cross sections of Andresen, Hausler, and Lulf (1984). The discrepancy 
between the expected and measured fluxes is resolved when we consider the effects 
of collisional de-excitation as a means of destroying photons. The 84 fim data will 
be analyzed along similar lines, and we will attempt to predict the 13.4 GHz flux as 
well. 

The analysis that will be carried out in this chapter does not carry with it 
any evolution due to shocks and shock chemistry. The latter considerations present 
a complex exercise that will hopefully be carried out at a later time. The present 
treatment, however, is still very important, since it can be used as a springboard in 
the understanding of the OH line fluxes, and also be used to check and provide limits 
to the MHD shock models that have been run recently by Chernoff and Draine. 

We did not have sufficient resolution to determine the line widths for the 84 
and 163 /xm line data. Since the 119 jxm lines show a width of 30 km s -1 , we shall 
assume that the other OH lines are also that wide. Recent results from Melnick 
(1985) suggest that the 163.12 /xm line is about this wide, or a bit wider. The 
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163.4 pm line in the data is cut off by a telluric water vapor feature, and so a width 
determination is impossible at this time for that line. 

5.4.2 The H 2 Density 


The 163 i*m data allow us to obtain an estimate of the H 2 density within the 
emitting region. The total 163 inn line flux is in the range 2.4-3. 8 x 10“ 17 W cm -2 . 
Using 520 parsecs as the distance to the Orion Nebula, we find that the total 163 /zm 
luminosity is 6.4-10.1 x 10 47 photons s -1 . Since we expect the OH to reside within 
the shock and post-shock region (from earlier arguments) , we can choose as a start- 
ing point for the temperature and density of the region those parameters that were 
found from the CO measurements, namely T ~ 750 K and 3 x 10® cm -3 . The 
dominant excitation mechanism at these temperatures and densities is collisional 
excitation by H 2 . Recent experimental values for the collisional cross-sections from 
the 2 n 3/2 J = j level to higher levels have been obtained by Andresen, Hausler, and 
Lulf (1984). These cross-sections are significantly different from those obtained the- 
oretically by Dewangan and Flower (1981,1983) and Schinke and Andresen (1984). 
The experimental cross-sections are listed in table 5.2. 

To find an expression for the molecular hydrogen density, we note that the 
collisional cross-sections from the ground state to excited states within the 2 n 1/ ^ 2 
ladder is roughly 3.6A 2 . Since the spontaneous emission rates are much higher than 
the collisional de-excitation rates, we may assume that most of the excited levels in 
the 2 n i y 2 ladder will result in a 163 /zm photon being released. We can then equate 
the collisional excitation rate directly to the 163 /zm luminosity. 

The oxygen to hydrogen ratio is given as 4.25 x 10~ 4 by number (Draine, 
Roberge, and Dalgarno 1983). According to the shock models of Draine and 
Roberge (1982,1984), the OH abundance will be most enhanced in a region that 
has a thickness of roughly 5 x 10 14 cm. The OH to O ratio in this region is close 
to \%. We will therefore assume that most of the OH emission will come from this 
thin region. If we let e = [OH]/[H 2 ], then c « 8.5 x 10 -6 . 
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Table 5.2 

H 2 -OH Collisional Cross Sections (in A 2 ) 

Final J 


2jI l/2 

1 + 
7 

— 

1.1 

1“ 

7 

— 

1.2 

3 + 
7 

— 

0.8 

3~ 

7 

— 

1.1 

5 + 
7 

2.5 

0.5 

5- 

7 

2.5 

0.75 

7 + 
7 

0.75 

0.15 

7 - 

7 

1.05 

0.3 

9+ 

7 

0.15 

0.04 

9“ 

7 

0.35 

0.08 


A list of the H 2 -OH collisional cross sections as obtained from experimental mea- 
surements for collision energies of 83meV (Schinke and Andresen 1984; Andresen, 
Hausler, and Lulf 1984). All of these measurements are for transitions from the 
2 n 3/2 J = § ground state. 



144 


Equating the excitation rate to the 163 luminosity, we may write 

(ovrtHj) {Nh 2 € f 0 ) = Li63 , (5.1) 

where a is the collisional excitation cross-section, equal to 3.6 A (see table 5.2); v 
is the average collision velocity; ngj 13 the molecular hydrogen density, N h 3 is the 
number of H2 molecules present in the emitting region; /o is the fraction of the 
OH molecules in the ground state; and L1Q3 is the photon luminosity in the 163 /xm 
lines. 

The observations were carried out using a l' square beam. If we let D be 
the size of the region represented by a side of the beam, and let d represent the 
thickness of the emitting region (we will assume this to be 5 x 10** cm), then we 
may write the volume of the emitting region to be 

V = ^ = D 2 dg { , (5.2) 

n H 2 

where gf is the filling factor for the beam. We assume that the OH-bearing region 
is uniform in the direction of d, as is suggested by the MHD models. The above 
two equations may be combined to yield 

«H 2 = (£l63)5 (D 2 dg { ovef 0 )~? . (5.3) 

For a kinetic temperature of 750 K, we find 

,n H2 « 2.6-4. 2 x 10 6 g f -1 ^ 2 / q 1/<2 cm -3 . (5.4) 

Using the OH data of Watson (1982), we calculate for the OH emitting region 
an expected molecular hydrogen density on the order of 6 x 10® cm -3 . We find a 
minimum expected value for the filling factor to thus be ~ 0.2. 

How good is the correspondence between collisional excitation into the 2 n 1 y 2 
ladder and 163 /im emission? If the 163 /xm line is optically thick, then you may 
not actually be measuring the collisional excitations into that ladder. To find the 
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optical depth, we need to find the cross-section for absorption. If nig3 is the number 
density of 163 jim photons, then the absorption rate must be a nie3 e. But this is 
equivalent to By J163, where By is the Einstein coefficient for absorption, and Jie 3 
is the mean intensity at 163 fim. The photon number density is related to J by 

= (5 - s) 


The cross-section for absorption can then be written as 


_ 9\ A 3 

• b * " (i 8irAv Aii ’ 


(5.6) 


where the various relations among the Einstein coefficients were used. The value 
of Ajj for the 163 /xm lines is 0.065 s -1 per parity transition (Brown et al. 1982). 
Assuming that Av is about 30 km s -1 (making use of the various line widths found 
for other molecules in the shocked region), we find the cross section to be 75 A 3 . 
The column density of OH in the ground state in the 2 n 1/ / 2 ladder cannot be greater 
than ngj id, the total OH column density, which is approximately 2.6 x lO 10 cm -2 . 
Hence, the optical depth for resonant scattering (absorption) at 163 /tm is roughly 
190. 

One must be careful, however, in interpreting this number. Although it is 
true that the 163 nm photon will be absorbed very quickly, it will also be reemitted 
via spontaneous and stimulated emission. The only way to destroy the photon is 
either via collisional de-excitation or emission of a 53 fim photon to the ground state 
of OH ( 2 n 1/2 J = | — ♦ 2 n 3/2 J = |). The ratio of the collisional de-excitation rate 
to the overall stimulated emission rate for the upper 163 nm level is on the order of 
10 -3 . This implies that only 1 in 1000 absorptions will result in a destruction of 
the photon. 

The number of absorptions or scatterings that a photon undergoes is pro- 
portional to the optical depth. Slater, Salpeter, and Wasserman (1982) have shown 
that the number of scatterings is nearly equal to twice the optical depth, or 


Wes « 2r abs . 


(5.7) 
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The reason, why the number of scatterings does not follow a random walk pattern 
(N^s <x r 2 ^) is that not only is there a random walk in space, but also in fre- 
quency. An optical depth of 190 will then result in roughly 380 scatterings and the 
proportion of these that lead to a collisional de-excitation is 0.001. Hence, most 
of the excited OH molecules will end up undergoing a photon-producing transition 
which eventually leads to the photon escaping. From this, we can assume that most 
excitations will lead to the release of a 163 nm photon. 

It should also be noted that if there is a strong velocity gradient within the 
medium, then we must be careful in determining the optical depth. Most of the 
absorption will be localized, and the effective column density of OH used will be 
much smaller, leading to a smaller optical depth and to a higher percentage of the 
photons escaping. 

One other component of the molecular cloud that we need to consider is 
the effect due to the dust component. The photometric observations of Werner et 
al. (1976) suggests an emissivity inversely proportional to wavelength and equal to 
unity at 20 /im. Hence, the optical depth of the dust at 163 /tm is estimated to be 
« 0.12, suggesting that the dust is optically thin and most of the 163 /tm photons 
will escape. 


5.4.3 The 119jim Flux: Problems and Solutions 

Is the observed 119^xm flux consistent with the 163 /zm data? We can obtain 
an estimate of the expected 2 II 3 / 2 J = § — 1 ► § flux in the following way. We stated 
above that collisional excitation is the dominant mechanism by which the excited 
states are populated. For the 2 n 1 / 2 ladder, most of the excited levels ( J > j) will 
radiatively decay down the ladder; and will ultimately result in the emission of a 
163 nm photon. The same argument can be used for the 2 n 3 / 2 ladder. The total 
collisional excitation cross-section for the 2 n 1/ / 2 J > j ladder is « 3.6 A 2 . For the 
2 n 3/2 ladder, a » 7.3 A 2 . Using (5.1), we note that the ratio of the 163 to 119 /im 
luminosities should be equal to the ratio of the collisional cross sections, since all 
of the other variables in (5.1) are independent of the transition we are considering. 



147 


Hence, if most of the photons escape, the 119 /xm luminosity should be 

Lh 9 « 1.3-2.0 x 10 48 photons s -1 . (5.8) 

From this figure, we find the expected flux at earth to be 6.7-10.3 x 10" 17 W cm -2 . 
We measure a total flux that is a factor of 5 to 7 smaller, implying that only about 
16% or so of the 119 pm photons actually escape the nebula. 

In order to better understand the different processes that can occur, we need 
to consider what is happening to the 119 jxm photons and the OH molecules. The 
spontaneous emission rate, Ay, for the 2 n 3 / 2 J = | — ► j transition is 0.138s" 1 
per parity state (Brown et al. 1982). This is one path that will lead to a transition 
from the first excited state to the ground state. A photon is released in the process. 

Another photon-producing process is stimulated emission. The emission rate 
is given by By Jug. Near the surface of the nebula, the mean intensity is simply 
related to the observed flux. The flux that we observe should be equal to the photon 
number density times the ratio of the object’s solid angle to An steradians, and then 
multiplied by the average component of the speed of light within the solid angle 
and along the path from the object to the observer, all of this then multiplied by 
the filling factor. This reduces to 


An XF\ 

nx ~ hc 2 6Clg, ’ 


(5.9) 


where n\ is the photon number density, F\ is the measured flux, and BO is the solid 
angle subtended by the beam. The mean intensity is related to the photon number 
density via 


J\ = 


h c 2 
An Av 


n A 


(5.10) 


Combining equations, we find 


J a *a 
A SVAvgf ' 


(5.11) 


There are two sources for J\, the OH 119 /xm line emission and the emission 
due to the dust grains. If is the filling factor due to the grains, then we find 
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v 

that the stimulated emission rate per parity transition is 

R* tun = 5.9 X lO" 4 [3.80 f - 1 + 8.5J7- 1 ) s _1 • (5.12) 

Similarly, we may find the absorption rate , which is 

flabs = 8.3 x 10" 4 [3.8(7 f _1 + B.5g~ l ] s" 1 . (5.13) 

We can obtain a very crude estimate for /q by assuming that only the first 
two levels are significantly populated. This assumption leads to 87% of the OH 
molecules residing in the ground state. We can expect a smaller fraction, and so 
will assume that roughly 75% of the molecules are in the ground state. The cross- 
section for absorption is given by (5.6), and is 43.6 A for the 119 nm photons. The 
optical depth for absorption is found to be 

^aba 55 > (5.14) 


where the range is ±23%. 

Once a 119 nm photon is absorbed, it may be reemitted via spontaneous 
emission. The photon will travel through the medium in a random fashion, being 
absorbed and then reemitted many times. As was stated above, the random walk 
occurs in space and in wavelength. According to Slater, Salpeter, and Wasserman 
(1982), the actual number of resonant scatters is very nearly twice the optical depth. 
For this transition, this remains valid until r^ a approaches values of 10 9 or greater. 
Hence, the number of scatterings is 

JV es (119) « r ia « 109g" 1/2 . (5.15) 

Collisional de-excitations can also lead to a transition to the ground state 
from the first excited state. The collisional excitation cross-section is roughly 2.5 A 2 
per parity transition (Andresen, Hausler, and Lulf 1984), which means that the de- 
excitation cross-section is roughly 2.1 A 2 . We have made use of equation (4.49) to 
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obtain this result. Therefore, the collisional de-excitation rate from the 2 n 3 j 2 J = | 
level to the ground state is given by 


*cd = ff cd V n H j 

« 2.2 x lO _4 0f 1/2 s* 1 • (5.16) 

The optical depth for true absorption is proportional to the number of reso- 
nant scatterings, N es« It is, in fact the number of resonant scatterings multiplied 
by the ratio of the collisional de-excitation (i.e., photon destroying de-excitation) 
rate to the photon conserving de-excitation rate, or 


’’true 88 ^ES 


■^cd 

■^stim "b ^ 


For the 119 nm radiation, this reduces to 


(5.17) 


10.7 

r ‘"“ * 1 + 2.23e + 62 5( ’ 


(5.18) 


where e is the ratio of the filling factor for OH to the filling factor for the continuum 
radiation from grains. 

We saw that from the 163 /*m data, we expect a factor of ~ 4 times as much 
119 Mm radiation as we actually observe. This implies that rtrue is on the order of 
1.8. We may solve for the OH filling factor, which we find to range from 0.05 to 0.17, 
depending on the actual 163 fxm flux. From this, we find the molecular hydrogen 
density to be approximately 1.2 x 10 7 cm -3 , nearly independent of the value that 
we choose for the 163 nm. flux. This implies a column density for hydrogen to be 
roughly 6 x 10 21 cm“ 2 and for OH to be 5 x 10 16 cm“ 2 . These numbers compare 
favorably with previous results of Townes et al. (1983) and Watson (1982). 


5.4.4 The 84 urn Flux 

We must now determine whether the 84 Mm data are consistent with the 
163 and 119 Mm lines and with the parameters which we have determined from the 
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latter meats urements. The Einstein coefficient for spontaneous emission of the 84 nm 
transition is Ay = 0.521s -1 (Brown et al. 1982). From (5.6), we can determine the 
cross-section for resonant absorption of 84 pm photons to be 

<7 abs (84) = 53.4 A . (5.19) 

We can then estimate the optical depth for absorption. To do this, we need to know 
the column density of OH molecules that are in the first excited state ( 2 Tig J = §)• 
The fraction of the population in the ground state was estimated to probably be 
75%. will assume that the fraction in the first excited state, / 1 , is roughly 20% 
or les' then follows that the optical depth for absorptiorTis 

r abs(84) « 54 . (5.20) 

m 

The collisional excitation cross section from the ground state to the J — j 
state is given as 0.9 A 2 (Andresen, Hausler, and Lulf 1984). We may estimate the 
collisional de-excitation rate to be roughly twice that number, or 1.8 A 2 . Then, the 
collisional de-excitation rate is roughly 7 x 10 -4 s -1 . The stimulated emission rate 
can be calculated as in the previous section, and we find that for 84 /im radiation, 
the stimulated emission rate is 

A,tun(84) * 3.1 x 10‘V 1 + 3.8 g~ 1 } s~ l . (5.21) 

The optical depth for true absorption is then found to be 


niue * 0-14 . (5.22) 

This implies that roughly 86% of the photons will escape the region, while the 
remainder will be destroyed. The collisional excitation rate from the ground state to 
the J = j state is approximately 3.8 x 10 -4 s -1 , whereas the spontaneous emission 
rate is over a factor of 1000 larger. We will assume that collisional excitation is 
the predominant process which populates this excited upper level. Then, the flux 
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that we expect to see from each of the 84 pm parity lines is found to be between 1.2 
and 3.3 x 10 -17 W cm -2 . We are nearer to the lower limit in the actual detections, 
which indicates that the 163 nm flux is probably closer to 1.2 x 10“ 17 W cm -2 . 

5.4.5 The 13.4 GHz Data 

One final set of OH data which we must consider is that of the 2 n 3 / 2 J = 
2 F = 4 + — ► 4“ A-doubling transition at 13.4 GHz (Matthews et al. 1984). The 
13.4 GHz line was found to have a width of 25 km s -1 , and the field of view was 
2(l. We may assume that most of the flux was coming from a more confined region 
comparable to our beam size of ~1 ; . The flux estimate for this particular transition 
is 0.09 Jy. This corresponds to a photon flux of 1.0 x 10 -2 photon cm -2 s -1 . From 
the 3K background, there is a photon density of 7 x 10“ 10 cm -3 . The spontaneous 
emission rate is 9.3 x 10 -9 s -1 , while the stimulated emission rate is 3.4 x 10 -8 s -1 
(Destombes et al. 1977). One can determine the expected flux from this transition 
by the following method. Assume that all of the 2 II 3 / 2 7 = j F = 4 + -+ 4" flux 
is coming from a l' region centered on BN-KL. We know that for the J= | level, 
the population fraction is about 20%. If we assume that the J = j level population 
is down by smother factor of ~5, then we expect the level to contain about 4% of 
the OH. Knowing the emission rate, we then can deduce that the expected flux is 
~1.2 x 10 -2 cm -2 s -1 . This takes into account the fact that only one polarization 
was observed. The measured flux is about 20% less than this value, suggesting that 
at least to first order we have a consistent picture of the OH emitting region. 

5.5 Shocked vs. Post-Shock Gas: A Quandary 

5.5.1 Statement of the Problem 

The FIR rotational lines of CO have been used to trace out the hot shocked 
gas ( T ~ 1000 K). These lines have excitation temperatures on the same order. 
According to Draine and Roberge (1982), an enhancement of OH in the shock front 
is to be expected, primarily due to accelerated production. Since the line widths of 
OH (as determined by Watson 1982) axe similar to those of the hot, shocked CO, 
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we may assume that these two molecular species are coexistent. 

This is not surprising, but an equally, if not more, important component was 
not included in the previous analysis. In the post-shock region which extends for 
roughly 10* 7 cm in radius (Genzel 1985), the gas is warm ( T ~100K, as determined 
by numerous observations of molecules such as HCN, SiO, NH3, etc.) and the 
OH abundance is roughly 1/40 that in the shock front. However, since the lower 
lying transitions that have been measured in OH have excitation temperatures of 
only 250 K or less, these molecules can contribute a great deal to the emission 
mechanisms. Since the total column density in this region is greater than that from 
the shock front, this component may actually dominate. 

This section will deal with this component of the region, proceeding through 
the same analysis as before, but for a warm gas in a region of size 1 x 10 17 cm. 

5.5.2 The Post-Shock Treatment 

Once again, we begin with the 163 fim flux. Using the same arguments as 
before, we find that 


n H2 » 1.87-2.35 x 10® g f 1/2 f Q 1/2 cm -3 . (5.23) 

Much of the difference in the numbers lies with the collisional excitation rate, which 
is proportional to t7, the average relative collision velocity. This is decreased by a 
factor of about 3, due to the decrease in temperature. 

When we redo the analysis for the 119 nm line emission, we find the following 
to occur. Since the thickness of the region, d, is much larger, the optical depth for 
resonant scattering is much larger. This will in turn drive the filling factor down, 
since the optical depth for true absorption remains unchanged. This is, indeed, 
what we see. The filling factor now varies over a range between 0.041 and 0.074. 
The molecular hydrogen density, however, is nearly constant over this range, and is 
given by 

mh 2 » 1.0 x 10 7 cm -3 . 


(5.24) 
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This leads to upper limits for the H 2 and OH column densities of 

Nr 2 = 1.0 x 10 24 cm -2 

N OE = 2.1 x 10 17 cm -2 (5.25) 

Again, the 84 /*m line is consistent with this picture, and implies that the 163 fim 
flux is probably closer to its lower value of 1.2 x 10 -17 W cm -2 . The 13.4 GHz flux 
is also predicted by this model. 

Although the hydrogen column density is high, this is due to the assumption 
that the region is uniform in the direction of d, which is now 10 17 cm long. Most 
likely, this is not the case. If the region is uniformly clumpy, then we can reduce 
the column density by a factor of ~ 3. The mass of hydrogen associated with this 
OH-emitting region is roughly 13 M©. This agrees rather well with estimates from 
CO lines and MHD models. More about this will be discussed in the next chapter. 

We see that we can explain the OH observations with two very different 
models. It is not clear which one is actually correct, or if both regions can be 
contributing to the observed fluxes. More data are needed in order to understand 
the OH emitting regions, especially data from the higher level transitions. MHD 
models of OMC-1 that have been run indicate that it is the post-shock gas that 
contributes most of the 119 nm photons. This is most likely also true of the 84 /zm 
and 163 ^m photons. These transitions may help trace out the OH emitting regions 
in the warm gas. 

5.6 A Measurement of Hot, Shocked CO in Orion: 

The J = 15 -*• 14 Observation 

5.6.1 Inspiration 

Carbon is one of the most abundant elements with Z>2 in the interstellar 
medium. Although it is not known how much carbon is tied up in grains, CO is one 
of the end products of carbon chemistry. CO is ubiquitous in the Galaxy, a fact 
that has enabled several surveys of the galactic plane to be conducted, looking at 
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the 2.6 mm J = 1 — ► 0 transition. 

The bond between the carbon and oxygen atoms is strong, requiring about 
11 eV to dissociate. Because of this, CO is the dominant carbon-bearing molecule 
in the ISM. CH, on the other hand, has a bond energy of about 4.2 eV. 

There has been a host of CO observations toward OMC-1, covering transi- 
tions with excitation temperatures of several degrees Kelvin to over 3000 K. These 
data tell us much about the molecular cloud and its structure. The high J transi- 
tions must arise in the shocked gas and from any post-shock gas that has not cooled 
considerably. 

The millimeter observations of the low J transitions of CO have shown a 
narrow spike and a broader plateau region in the spectra. The spike has a width 
of «5km s -1 , whereas the plateau region has a FWHM from 30 to 150 km s -1 
(Schultz et al. 1985; Phillips, White, and Watt 1982; Zuckermam, Kuiper, and 
Rodriguez Kuiper 1976). The spike is ubiquitous in a region covering many axe 
minutes in angular extent (Zuckermam and Palmer 1975). In contrast, the plateau 
region which is centered on the BN object is about 36" in size when measured in 
the J = 1 — ► 0 transition (Solomon, Huguenin, and Scoville 1981). Estimates of 
the size of the plateau region by interferometric means using CO brings this value 
down to roughly 17" (Masson et al. 1984). For higher J transitions, the source 
size seems to be appreciably larger, — 1 ; .5 (Storey et al. 1981). The spike feature 
is not pronounced toward the BN region for J > 6. This is indicative of the spike 
originating in a cool, quiescent region. The fact that high J transitions ( J = 34) 
have been seen (Watson et al. 1985) implies that the plateau feature originates in a 
hot part of the cloud. Since H 2 measurements show regions of peak intensity near 
the edge of the plateau emitting region (Solomon, Huguenin, and Scoville 1981; 
Beckwith et al. 1978), this strongly suggests that this hot component is due to 
shocked or post-shock gas at temperatures up to 2000 K. 

McKee et al. (1982) have shown that the higher J transitions for CO should 
be optically thin. One of the conclusions of their paper is that the CO transitions 
for J ~ 15 are rather insensitive to temperature and density, once the critical values 
of both parameters are reached. Such is the case in Orion. Hence, it is possible to 
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obtain from the measurement which we report in this paper the mass of hot shocked 
CO in the nebula. 

The CO observations that were available in 1982 led Watson (1982), Chernoff, 
Hollenbach, and McKee (1982), and Draine and Roberge (1982) to attempt to fit 
the data with a model. Watson’s model consisted of a two component system 
(warm and hot, corresponding to post-shock and shock) while the other two models 
employed MHD shocks with shock chemistry and evolution. The results showed 
that a fit to the data was possible. Since the low J transitions arise from the whole 
molecular cloud — hot, warm, and cold gas — and since these lines are optically thick, 
the fits to the data from these regions are sometimes not very good. 

We were inspired to make the 173 fim observation for three reasons. The first 
is that according to the fit to the data, the shocked region has a temperature of 
roughly 1000 K. In this temperature region, the intensity of this CO transition is 
nearly independent of density over the range that is expected for Orion (McKee et 
al. 1982). We can determine the amount of molecular hydrogen that is associated, 
with the hot CO by finding the amount of CO present, and then extrapolating to 
hydrogen. This was first done by Stacey et al. (1983), using the J = 17 — ► 16 
transition at 153.3 jim. 

The second reason for observing the 173 /xm line was to extend the FIR 
observations that have taken place to a longer wavelength. Before this observation, 
no detections of spectral lines longward of 163 nm had been made in the FIR. 
This detection constitutes the first beyond 163 /xm and, as a result, is the first 
astronomical detection of the 173 j*m CO feature. 

The third reason that this line was observed was for calibration purposes for 
our instruments. Since this line was expected to be bright, it was an obvious choice 
to use in order to help us calibrate the grating spectrometer and interferometer. 
During the same flight series, we also observed several other objects, some at wave- 
lengths out to 175 nm. Since these lines were expected to be weak, the calibration 
of the grating spectrometer was essential to have optimum throughput. 
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5.6.2 The Observations 

On the night of 1985 January 23-24, we observed the BN-KL region using 
the KAO, flying at an altitude of 12.5 km. The beam-size was l' square. During 
the integration time, some mild turbulence was experienced, causing the beam to 
wander about the source. This meandering was only occasional, and was no larger 
than roughly 15" at worst. Since the CO source is believed to be slightly larger 
than l' , this should not have affected the results very much. However, turbulence 
in the atmosphere means changing boresight water vapor, and that can introduce 
fluctuations in the signal strength. 

Figure 5.7 shows the resulting spectrum at this wavelength. Telluric water 
vapor features occur near 173.34 and 173.72 /xm. In addition, a broad feature lies 
between 174.24 and 175.09 /xm. The calibration of the spectrum was obtained by 
calibrating on the continuum of Orion (Werner et al. 1976). This calibration was 
confirmed by a calibration of Venus at this wavelength. The total flux in the line is 
4.85 x 10“ 17 W cm -2 . 


5.6.3 Interpretation 

The excitation temperature of the J = 15 rotational level is roughly 10 3 K. 
We therefore expect the 173 /xm emission to arise from the hot, shocked CO in and 
near the shock front. The size of the CO emitting region is roughly l(5 (Storey et 
al. 1981). This implies that the total 173 /xm flux that we should measure is about 
80% higher than our observed value of 4.85 x 10 -17 W cm -2 . From this, we deduce 
that the total photon luminosity of OMC-1 at 173 /xm is —7.5 x 10 48 photons s -1 . 

For temperatures and densities that are likely to exist in the shock region, 
McKee et al. (1982) find that the specific intensity of the CO molecules is «1.78 x 
10 -20 erg s -1 str -1 molecule -1 . From this we can deduce that there must be 
—1.2 x 10 53 CO molecules in the emitting region. If we assume cosmic abundance 
of carbon and that approximately 25% of the carbon is in the form of CO (the rest 
being in grains and other carbon molecules), then the total mass of H 2 is found 
to be ssl.OM©. This compares favorably with the 1.5 Mq which was found for the 
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Wavelength (/*m) 


FIGURE 5.7 Spectrum of BN-KL at 173 nm. The J = 15 — ► 14 CO emission line 
at 173.631 /xm arises in the hot shocked gas surrounding IRc2. Data were obtained 
in 1985 January with an integration time of 20 seconds per point. All data were 
obtained using a V beam. 
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J = 17 — * 16 transition at 153.3 /im (Stacey et al. 1983), and with the 1.0 M© given 
by Watson et al. (1985). 

5.7 Discussion of Results 
5.7.1 The Present Findings 

The CO observation fits in well with the rest of the CO data that has been 
gathered to date on Orion. Table 5.3 lists these data and they are displayed in 
figure 5.8. We see that the 173 /im line lies on the curve that is a model fit to the 
data (Watson et al. 1985). Since the region around J ~ 16 is nearly independent 
of density, this line tells us little by itself. It is only when we use it with the rest of 
the data that we can determine temperatures and densities, and specific intensities 
from McKee’s models. We have been able to determine the mass of hot, shocked 
gas that is in the Orion nebula. This mass turns out to be roughly 1M©. We may 
determine the filling factor of CO by using the results of Chernoff, Hollenbach, and 
McKee (1982). They have found that the shock region contains about 8M© of gas. 
This implies that the filling factor is roughly 0.13. 

If we assume that the shock front forms a roughly spherical shell around the 
outflow region, with a size of ~ 10 17 cm, then we cam determine the thickness of this 
region by using the njj 2 densities obtained for the shock region by Watson et al. 
(1985) and comparing the volume of the shell to that containing 1 M© of hydrogen. 
This thickness turns out to be ~3 x lO 13 0f -1 cm, where gf is the filling factor for 
the CO. Plugging in, we get a shell thickness of 2.4 x 10 14 cm. That is to say, the 
hot, shocked CO most likely lies within a shell of minimum thickness 2.4 x 10 14 cm. 
The reason why this is a lower limit is due to the fact that there may be clumping, 
as suggested by Wynn- Williams et al. (1984) and Geballe et al. (1986). If g z \ is 
the clumping factor (the ratio of the total clump volume to the total volume) , then 
the thickness of the shell will be inversely proportional to g c \. Draine and Roberge 
(1982) expect a shock front to have a thickness on the order of 2 x 10 15 cm. If this 
is the case, then we can estimate a clumping factor on the order of 0.1. 

The OH data do not lend themselves so easily to interpretation as do the 



Table 5.3 CO Observations of OMC-1 

Transition 

Wavelength 
(in nm) 

Beam Size 

Intensity 
(W cm -2 str -1 ) 

Reference 

1-0 

2600.76 

65 " 

5.0 x 10 -14 

1 



45 " 

3.3 x 10~ 14 

2 

2-1 

1300.40 

1 ' 

6.4 x 10“ 13 

3 



1(7 

1.8 x 10 -12 

4 

3-2 

866.963 

55 " 

3.6 x 10“ 12 

5 

6-5 

433.556 

35 " 

2.0 x 10“ 10 

6 



35 " 

5.8 X 10“ 10 

7 

7-6 

371.650 

34 " 

1.2 x 10 _1 ° 

8 

15-14 

173.631 

l ' xl ' 

5.7 x 10" 10 

9 

16-15 

162.812 

l ' xl ' 

7.6 x 10“ 10 

10 



1 ' 

7.5 x 10“ 10 

11 

17-16 

153.267 

l ' xl ' 

8.3 x 10“ 10 

12 

21-20 

124.193 

1 ' 

1.3 x 10 -9 

13 



60 " 

1.0 x 10“ 9 

14 

22-21 

118.581 

1 ' 

2.1 x 10" 9 

13 



44 " 

1.1 x 10" 9 

15 

26-25 

100.461 

44 " 

3.8 x 10 _1 ° 

15 

27-26 

96.7725 

60 " 

6.5 x 10“ 10 

14 

30-29 

87.190 

60 " 

2.4 x 10“ 10 

14 

31-30 

84.4107 

30 " 

1 . 8 - 2. 6 x 10 _1 ° 

15 

34-33 

77.0587 

44 " 

8.6 x 10" 11 

15 


A list of the published CO observations to date toward the BN-KL region of the 
Orion Molecular Cloud. The reference numbers refer to the following articles: 1) 
Zuckerman, Kuiper, and Rodriguez Kuiper (1976); 2) Scoville (1981); 3) Wannier 
and Phillips (1977); 4) Phillips, White, and Watt (1982); 5) Richardson et al. 
(1985); 6) Goldsmith et al. (1981); 7) Koepf et al. (1982); 8) Schultz et al. (1985); 
9) Viscuso (1985); 10) Stacey et al. (1983); 11) Crawford et al. (1986); 12) Stacey 
et al. (1982); 13) Watson et al. (1980); 14) Storey et al. (1981); 15) Watson et al. 
(1985). 
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ORION BN-KL CO J-J-1 Emission 



FIGURE 5.8 Plot of the CO intensity towards BN-KL versus the transition ( J — ► 
J — 1). The solid curve is a theoretical fit to the data from the MHD model of 
Draine and Roberge (1982). The fit is worse for low J transitions, owing partly 
to the increasing optical depth with decreasing J. Observational data are listed in 
table 5.3. 
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CO data. As we saw above, it is possible to construct two different models of the 
region that fit the data. Whereas the interpretation that the OH emission arises 
from the shocked gas was discussed in Viscuso et al. (1985a), Chernoff (1986) favors 
the post-shock gas as the source of the OH lines. It is unclear which regime makes 
the more significant contribution, but most likely both are important. 

Clearly, the OH lines resulting from highly excited levels will trace out the 
hot component of the gas. These lines, however, have not been detected with any 
degree of assurance. 

The CO data have helped enormously in directing the way for MHD shock 
calculations. The model that Watson used in his thesis is based upon a simple 
escape probability method, which neglects the shock structure and evolution. The 
fits were not strongly dependent upon temperature, leading to the approximation 
of an isothermal slab. 

Velocity gradients do exist within the flow region, as determined by H 2 O 
masers. Since much of the OH emission may be coming from the post-shock gas, 
we expect velocity gradients. In fact, it is most likely due to these gradients that 
the observed line widths (30 km s _1 ) are a factor of 80 bigger than those expected 
from random thermal motions alone. The treatment of the region is not necessarily 
amenable to such a simple radiative transfer model as that above. Following the 
CO work, it is desirable to perform a numerical radiative transfer study on the 
OH lines that have already been observed. Such a study can take into account the 
velocity gradient, as well as incorporate a multi-level molecular structure. The next 
chapter will deal with just such an endeavor. The goal is to be able to determine 
whether the OH emission is from the thin shock front or the much more extended 
post-shock region. 

We have demonstrated that it is possible to determine some of the parameters 
in the OH emitting region. The density of the molecular hydrogen is roughly 1.0 x 
10 7 cm -3 , almost independent of temperature. Collisional de-excitation plays a 
very important role in the destruction of the 119 pm photons, and a minor role in 
the frustration of the 84 and 163 nm flux. 
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5.7.2 Predictions for the Future 

With better signal-to-noise ratios for the 163 and 119 /xm lines, we could 
estimate the relative line fluxes between, the two different parity states. From the 
collisional cross-sections, we should expect the 163.4 ^tm line to be »70% as strong 
as the 163.1 txm transition. For the ^3/2 J - § § transitions, the 119.2 line 

flux will be «87% as strong as the 119.4/xm flux. Since the 2 n 1 /2 ladder decays 
slowly to the ground state, we also expect that the lowest levels in this ladder will 
be populated more than anticipated, and therefore, the 53 fxm and the 96 pm fluxes 
should be higher than the ratios of the Einstein coefficients alone would lead us to 
expect. 


5.8 Other Observations in Orion 

During the 1983, 1984, and 1985 winter flight series aboard the KAO, we 
also undertook a program whereby we searched for FIR emission from several other 
molecular species, namely CH and CH2. These two species are important links in the 
carbon chemistry. The chemistry of shock regions is not very well understood, but 
the CH and CH2 fluxes are not expected to be very large owing to small abundances. 
The carbon chemistry suggests that the CH would readily combine with any atomic 
oxygen to form CO. This would imply that very little CH would exist. Rydbeck et 
al. (1976) have found that CH is ~4 times less abundant than OH in extended Hll 
regions; however, with the shocked region in OMC-1, the CH/OH abundance ratio 
could be wildly different from 1/4. This suggests that the FIR emission lines from 
these species may be near the limit of detectability. 

The observations resulted in no positive detections of any of the lines. We 
were able to place an upper limit of 1.0 x 10“ 17 W cm“ 2 on the individual fluxes 
from the 2 n 1 / 2 J — § — 1 ► j CH transitions at 149.09 and 149.39 /xm (Brown and 
Evenson 1983) and the 3 Bi A^Kc = $12 ^03 CH2 transitions at 148.846, 148.933, 

and 149.188 /xm (Sears et al. 1984). 



VI. Escape Probability and OH in OMC-1 


6.1 Inspiration 

In chapter V, we saw how it was possible to apply radiative transfer to 
perform back-of-the-envelope calculations on the OH data presently available for 
OMC-1. The molecular hydrogen density could be determined quite well, but the 
temperature of the OH-emitting region remains uncertain to within factors of 10 or 
more. The filling factor is also not very well determined, albeit its fractional range 
is smaller than that for the temperature. 

This state of affairs is not surprising, especially in light of the simplified 
treatment used in chapter V. In order to more accurately determine the physical 
parameters of the OH-bearing region, it is necessary to incorporate interactions 
among all of the populated energy levels in the analysis. This is best done by 
numerical methods. We expect that the results of the numerical treatment should 
be close to those of the back-of-the-envelope treatment if the assumptions that were 
made in the previous chapter are valid. 

It is our goal to accurately fit the OH data from OMC-1 by using a one- or 
two-component model. This model will give us average molecular hydrogen densities 
and average temperatures for the OH-emitting region. The most straight-forward 
method used in solving the radiative transfer equations for multi-level, non-LTE 
(Local Thermodynamic Equilibrium) systems is the escape probability formalism 
outlined in chapter IV. 

6.2 The Radiative Transfer Code 
6.2.1 Modus Operandi 

The computer code that was written for this task began as an effort by 
Sean Colgan here at Cornell for analyzing OH maser and H 2 CO (formaldahyde) 
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data. This program was used as a template for our code. We have completely 
rewritten the code, implementing some streamlining changes and improving the 
overall rate of convergence to a solution. The program outline follows below. 

The equations to be solved are those of (4.79) and (4.80). We need to deter- 
mine the level populations for a given kinetic temperature and molecular hydrogen 
density. Other parameters that must also be specified are the total OH-column 
density, the OH line widths, and any dust or background sources of radiation and 
their associated parameters, and the charged particle density. The latter sources 
are known, and we include the dust parameters of Tielens and Hollenbach (1985) 
and Pipher, Duthie, and Savedoff (1978) in our calculations. 

"\e solution for which we search is the steady-st te solution, in which the 
populations of each energy level approaches a constant *ue. If we assume that we 
have k levels, then the rate of change in the population of each level must be zero 
in the steady-state solution, or 



i = 1, . . . ,k 


( 6 . 1 ) 


If we replace n, by /i(ni, . . . , n k ), then the set of equations in (6.1) can be rewritten 
as 


/l(ni,... ,n k ) = 0 
/ 2 (ni,... ,n k ) = 0 

/ k (ni,...,n k ) = 0 . (6.2) 

To solve this non-linear coupled set of equations for the level populations, we expand 
the functions in a Taylor series. We shall retain only terms of first order to facilitate 
the computation of the populations. If we begin with an initial point in solution 
space, {oi}, then we have 
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- . . v , dfl e , d fl c , , 5/l c 

0 = Mau...,a k ) + ^6i + — (, + ■■■ + — St 


0 = /2(fll “ k) + + + ” ' + jty^ 


0 = /k(l1 ak) + l^ {, + l^ <2 + "' + ^ <k 


(6.3) 


where the 6[ represent the deviations from aj. We can now compute the new ap- 
proximation to our solution from 


n i = a i + ^i i = l,...,k. (6.4) 

Using this new solution, we can then go back and repeat the process until the 
appropriate degree of convergence is achieved. One should note that this process 
does not guarantee convergence to a solution. 

The usual starting point is to assume that we have LTE. From here, we then 
calculate the partial derivatives and solve the resulting linear set of equations for the 
This solution is accomplished by using a decomposition and backsubstitution 
scheme. This scheme was borrowed from Press et al. (1986). We then find the 
new {r^}, and repeat the process. For more information concerning the solution of 
coupled-equations, see Acton (1970) and Press et al. (1986). 

Once we have found the solution to the level populations, we use equation 
(4.80) to derive the expected line intensities for all of the allowed transitions. 

6.2.2 Molecular Parameters 

The code as it now stands is capable of obtaining solutions for one of three 
molecular species, OH, CO, or H 2 CO. In this thesis, we are concerned only with the 
first two molecules. Some of the molecular parameters for OH and CO have already 
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been discussed in chapters IV and V. We briefly mention here other processes that 
are also necessary to consider. 

These involve collisional excitation and de-excitation rates, primarily for col- 
lisions with charged particles. For OH, charged-particle collisional excitation rates 
for transitions in the A-doublets have been calculated by Bouloy and Omont (1977, 
1979). Charged-particle collisional rates for the transitions between arbitrary levels 
have also been calculated by Goss and Field (1968). These are used in determining 
the effects of charged particles in the molecular cloud. 

For CO, new and improved collisional excitation rates by molecular hydrogen 
have recently been published by Schinke et al. (1985). They include rates from 
J = 0 to L. From these rates, we have been able to use the infinite-order sudden 
approximation (IOS) of Goldflam, Green, and Kouri (1977) to calculate the inter- 
level rates from J to J ' . The IOS approximation is valid for changes in rotational 
energy much less than the collisional kinetic energy. This means that for large A J, 
the method breaks down. However, DePristo et al. (1979) showed that it is possible 
to use the simple IOS formalism with a scaling factor to take the energy difference 
between levels into account and arrive at collision rates that match those found for 
more exact calculations. We plan to use these new rates in the near future to redo 
the fits of McKee et al. (1982) and Watson (1982) to the Orion CO data. 

6.3 Preliminary Analysis 

It is useful in the interpretation of the results from the code to first determine 
what we may expect to find under certain conditions. For now, we shall be concerned 
with the 163, 119, and 84 fim data, in particular, the 163:119 /im and 84:119 /x m flux 
ratios. In this chapter, we shall continue to use the convention that rij refers to the 
i th level population divided by the degeneracy. 

It is useful to first determine the flux ratios in the case of LTE. The OH 
will reach LTE if it is strongly coupled to the gas through collisions, i.e., collisional 
excitation and de-excitation dominate over radiative processes in the final solution 
for the level populations. Clearly, for a given temperature, if the density is low, 
then collisions will be few and spontaneous emission will depopulate levels. If the 
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density is very high, the collisions will occur frequently enough that the population 
of the levels is little affected by radiative processes. There is a critical density where 
the collisional excitation rate of a level equals the spontaneous emission rate. For 
level i, we may write this condition as 

) ] ffj n j n crit 7ji = ^ ] 9i n i -^ij > (®-5) 

j<» j<i 

where n cr j t is the critical density of the colliders, which in our situation are hydrogen 
molecules. The H 2 -OH collisional cross-sections are on the order of lA 2 (Schinke 
and Andresen 1984), meaning that the collisional excitation rates are on the order of 
10 ~ 12 <pl/2 cm 3 g-l^ w h e re T is the kinetic temperature. The Einstein A coefficients 
for the infrared transitions are on the order of 0.1-1 s -1 (Brown et al. 1982), leaving 
us with an estimate of the critical density of 

n cr it ~ 10 12 T -1 / 2 cm -3 . (6.6) 

For temperatures between 100 and 1000 K, we expect to approach LTE for densities 
greater than 3 x 10 10 cm -3 . From our previous analysis of the OH data and the 
CO results discussed in chapter V, we don’t expect densities much greater than 
10 7 cm -3 . Hence, the OH level populations most likely will not approach LTE. 

We would also like to know when the collisional rates become unimportant 
altogether. We cam estimate this lower limit to the critical density by looking at the 
transitions betweeen the significamtly populated levels that have small spontaneous 
emission rates. These are, for the most pairt, the A-doubling radio tramsitions. For 
the lower levels, these A coefficients are on the order of 10 -9 or 10“ 10 . Hence, the 
density when collisions begin to contribute to the level populations is roughly 

n min ~ 10 3 T -1 / 2 cm -3 . (6.7) 

For the expected range of temperatures in OMC-1, n min is on the order of 50 to 
100 cm -3 . This is at least three orders of magnitude below the inferred density in 
the region. 
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It is instructive to find the LTE ratios among the 163, 119, and 84 ^tm lines. 
This is accomplished in the following manner. Figure 6.1 shows the energy level 
diagram that we need to consider. This is a simpler version of figure 4.3, for we are 
not considering hyperfine splitting or parity splitting of the rotational levels. Using 
equation (4.80) and neglecting the contribution from the microwave background, 
we find that 

he 3 = / 119\ 3 / 1 — e~ fl63 \ / n 4 W n 1 -n 2 \ 
hl9 \1 63/ \1 - c~ f n9 J \ri3 - m J \ ri2 / 

where ri 19 is the optical depth in the 119 /im line, and so on. Since we are assuming 
LTE, we may simplify these equations by noting that in LTE, 


nj oc e Ei / kT , (6.9) 

where T is the kinetic temperature and E\ is the energy of the i 1 * 1 level. Using (6.9), 
(6.8) becomes 


hez _ /119\ 3 /1 -e _r i63\ / e (E2-Ei)/kT _ x \ 
hl9 ~ \ 163 ) {l-e-rntj \ e (E A -E 3 )/kT _ 1 ) ’ 

Ig 4 _ /119\ 3 / 1 -e~ f «< \ / e ^~ E ^/ kT - 1 \ 
/ll9 \ 84 / \1 — e _r n9 ) y e (Es~E 2 )/kT — ‘ 


( 6 . 10 ) 


In the case where all three lines are optically thick, the second term in the 
right-hand side of (6.l0) becomes unity as the source function approaches the black- 
body function (KirchofF’s law), and can be written as 


lm Ai9\ V M /r-i 

hi9 \163/ e 88 / r — 1 

/ 84 __ /119\ 3 e 12 °/ r - 1 
/l!9 “ V 84 ) e nO/T _ ! ’ 


( 6 . 11 ) 
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OH 2 n Vibrational Ground Sfote 



FIGURE 6. 1 Simplified OH level diagram used in the determination of LTE level 
populations. The A-doubling of the rotational levels has been ignored. The numbers 
next to the levels on the outside indicate the level number. (Ground state is level 
1, first excited state is level 2, etc.) 
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where the numerical values in the exponentials represent the energy difference be- 
tween the appropriate levels in degrees Kelvin. Figure 6.2 shows a diagram of the 
LTE flux ratios as a function of temperature in the optically thick regime. (The 
flux ratio differs from the intensity ratio by a factor of (119/163) or (119/84), owing 
to the inclusion of the line width in determining the flux. We assume a line width 
that is constant in velocity units.) Above 500 K, the remaining exponentials may 
be approximated with a MacClaurin series in first order, yielding 


Ae 3 
Al9 


/119\2 

*84 ^ 

/U9x 

V163/ 

A 19 ~ 

V 84 / 


( 6 . 12 ) 


for a hot, optically thick column of OH in LTE. This asymptotic behavior is shown 
in figure 6.2. 

For the flux ratios to be on the order of unity, we note that the temperature 
must be on the order of 30 K. If the kinetic temperature is warmer, then the 163 /tm 
flux is suppressed relative to the 119 /im emission and the 84 ^m flux is enhanced. 
Conversely, for lower temperatures, the opposite effect is observed. 

This analysis suggests that the OH radiation may be coming from the cool 
post-shock gas at a temperature of about 30 K. I must emphasize that this is merely 
a suggestion and not a conclusion. The above analysis assumes that the OH is in 
LTE and that all of the FIR lines are optically thick. The former premise is most 
likely false, because we must have densities at least 3 orders of magnitude higher 
than have been deduced from measurements to gain LTE. Most likely, the 119 and 
84 ^m lines are optically thick, but the 163 pm lines may have optical depths of 
order unity or less. 

At this point in the analysis, we must invoke the aid of numerical methods. 
There are many different contributors which determine the final populations of the 
OH energy levels, including the microwave background, the interstellar radiation 
field (ISRF) impinging upon the molecular cloud from IRc2, BN, KL, all of the 
other star formation pockets within the region, the Trapezium, and the overall ISRF 
from the rest of the Galaxy, the radiation field from the warm dust grains within 
OMC-1, and charged particle collisional excitation, in addition to the spontaneous 
emission and H 2 collisional processes. 



Log Flux Ratio 
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LTE FIR Line Flux Ratios for Optically Thick OH 



Log Temperature 


FIGURE 6.2 Flux ratios for optically thick OH in LTE. The ratios are for the 163 
to 119 #xxn and 84 to 119 emission lines. We assume line widths of 50 km s . 
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6.4 Numerical Results 
6.4.1 Limits to the Physical Parameters 

We need to determine some physical limits to the parameters that we can 
adjust. These parameters are fivefold: the kinetic temperature, the molecular hy- 
drogen density, the electron density, the radiation field from the imbedded stars, 
and the OH column density divided by the OH linewidth. We asstime that the 3 
FIR lines originate from the same OH, and so their linewidths in km s -1 should be 
roughly comparable to each other. 

The upper limit to the kinetic temperature comes from the H 2 and CO 
observations of Beckwith et al. (1978) and Watson (1982). This limit is between 
1000 and 2000 K, the temperature of the gas in and near the shock front. The OH 
will exist only in the molecular cloud, so we can discount the high temperatures 
(> 10 4 K) associated with Hll regions. We do not expect the gas to be coupled 
directly to the dust, which is heated by the radiation from the nascent stars, and 
so a lower limit to the gas temperature can be below 70 K. A lower limit that we 
will use for the cool post-shock gas is 15 K. 

Previous estimates of the hydrogen density are roughly 10 5 cm -3 in the pre- 
shock gas (Chernoff, Hollenbach, and McKee 1982; Draine and Roberge 1982) and 
~10 6 cm -3 in the hot, post-shock gas (Watson 1982). If we assume pressure equi- 
librium for the cool post-shock gas, the density could be as high as 10® cm -3 . In 
order to observe the effects of the collisional processes, we shall include densities 
from 10 3 to lO^cm -3 . At the lower end, the collisions should be negligible in the 
determination of the population levels. At the upper end, collisions are beginning 
to totally dominate the radiative processes, and we should begin to approach LTE. 

The electron density in the molecular gas is relatively small. However, it can 
still contribute significantly to the final level populations if the gas is cold. Tielens 
and Hollenbach (1985) estimate the electron population to be between 10 -4 and 
10 -6 of the hydrogen density within the molecular cloud. We will assume that this 
range is correct. Most of the electrons deep within a molecular cloud are generated 
through the ionization of He and H 2 by cosmic rays, as well as photoionization of 
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heavy metals such as Mg and Fe. Charge exchange reactions result in the He and 
H2 reg aining their neutrality while the metals become charged. 

The radiation field has several components. The first component is that due 
to the dust grains. This continuum component has been measured photometrically 
by Werner et al. (1976); Pipher, Duthie, and Savedoff (1978); and Erickson et al. 
(1981). The dust is heated by the radiation from nearby stars and converts that 
radiation to infrared radiation. It is this converted light that gets through to the 
inner regions of the molecular cloud. This inner region only ‘sees’ the young O 
and B stars by the grain radiation. To estimate the contribution from the general 
interstellar radiation field, we use the results of Mathis, Mezger, and Panagia (1983). 

Another parameter which needs upper and lower limits is the OH column 
density divided by the line width in km s” 1 . These two variables cannot be sep- 
arated. However, there are measured values of the line width, yielding somewhere 
between 30 and 50km s _1 (Crawford et al. 1986; Watson et al. 1985). If the OH 
emission originates in the post-shock gas, then we do not expect the OH-emitting 
region to be more than 0.2 parsec long (Genzel 1985; Werner et al. 1976). The 
post-shock OH abundance is on the order of 2.1 x 10~ 7 that of molecular hydrogen 
(see Chapter V). If the hydrogen density is on the order of 10 7 , then the resulting 
OH column density should be less than or about 10 18 cm -2 . We shall use values in 
the range of 5 x 10 14 to 5 x 10 18 cm -2 . The lower values will help to demonstrate 
the importance of optical depth in the final solution for the level populations and 
flux ratios. 

Another limit that we can place on the OH column density and abundance 
ratio is derived from the mass of H2 that is present in the shocked and post-shock 
region. If we assume that the OH emission arises from an area roughly l' in size, 
then we can show that the mass of molecular hydrogen present is given by 

= 3 x lO _22 0f M© , (6.13) 

where gf is the beam filling factor and e is the [OH]/[H2] ratio. 

We can estimate the mass of hydrogen from CO measurements. In the hot, 
post-shock gas, we know from the analysis in chapter V of this thesis that there is 
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roughly 1M© of H 2 . This is derived from the high J rotational transitions. The 
low J transitions will measure not only the hot gas, but also the cooler gas as well. 
Looking at the optically thin J = 6 — ► 5 transition at 433.556 nm, we can relate the 
measured line intensity to the column density of CO by 


T Nco ffi) Mshl/Q5 
7 65 = 71 


(6.14) 


where /(6) is the fraction of the CO in the J = 6 rotational level. From (4.19), we 
know that Aes is about 2.2 x 10 -5 s -1 . For temperatures ranging from 35 to 300 K, 
the fraction of the population in the sixth rotational level remains very close to 10% 
(assuming LTE). Koepf et al. (1982) found the corrected antenna temperature of 
this line to be about 180 K. The line width was 38 km s -1 , yielding an intensity of 
2.3 x 10 -10 W cm -2 str -1 . The deduced CO column density is 2.9 x 10 18 cm“ 2 . If 
the total size of the region is roughly 0.2 parsec, then the total molecular hydrogen 
associated with the CO is on the order of 8Mq. This assumes that the [CO]/[H] 
ratio is on the order of 10 -4 . Since this ratio may be much smaller if more carbon 
is tied up in grains than we expect, this represents a lower limit to the mass of 
hydrogen. For the moment, we shall assume that 10 M© of molecular hydrogen is 
present in the warm, CO-bearing post-shock gas. This mass compares favorably 
with that deduced from J = 1 — ♦ 0 CO observations (Kwan and Scoville 1976). 
The mass of hot, shocked and post-shock gas is seen to be negligible. 

Assuming that the observed OH is confined to a region roughly l' in size as 
the data suggest (Watson et al. 1985), we should expect to find the mass of hydrogen 
associated with the OH to be comparable to that for the CO. We find a lower limit 
for the hydrogen mass when we assume that most of the OH FIR radiation originates 
in the thin sheet within the shock front in which the OH abundance is enhanced. If 
this sheet is 5 x 10 14 cm thick, and extends for l' in two dimensions, then by using 
the density and filling factor that we found in chapter V, we obtain a minimum 
hydrogen mass of roughly 0.17 M©. Remembering that the OH is enhanced in the 
shock front as opposed to the post-shock gas, we find that this enhancement nearly 
cancels out the hydrogen mass difference when one is searching for constraints on 



the column density of OH. We find that 


JV 0 H ~ 1(j16 1 cm 2 . 


(6.15) 


6.4.2 Analysis of Results 

Various runs of the code were made covering the range of parameters listed 
in the preceding section. These runs were made for a 14 level system, which is quite 
adequate for temperatures up to about 1000 K. If the kinetic temperature is higher 
than 10 3 K, then there will be significant population of the higher levels. For these 
high temperatures, we used a 20 level system. 

The Newton method used to find the solution does not guarantee conver- 
gence. Although variations on the theme were tried in order to maximize the per- 
centage of converging solutions, this maximum did not exceed 90% for more than 
1000 runs. A limit of 50 iterations was placed in the code in order to minimize com- 
putational time. This makes it possible that some of these non-converging solutions 
would have converged if given enough time. Various self-checks in the code guar- 
anteed a non-diverging solution, but these could not guarantee proper convergence. 
There are several instances where the solution was found for one set of parameters, 
and a small change of only one input parameter led to a non-converging solution. 
In these cases, it is sometimes obvious what the final result should be. 

Two of the parameters were found to have an insignificant influence in the 
determination of the level populations. These are the electron abundance and the 
ISRF. Electron abundances of 10“ 4 or greater will affect the population significantly, 
but such high abundances are not expected in the region. The general ISRF (not to 
be confused with the radiation field due to the heated dust grains) also contributes 
little to the solution. We may therefore ignore these two mechanisms in the ensuing 
discussion. 

This leaves us with 3 parameters to vary, the kinetic temperature, the molec- 
ular hydrogen density, and the OH column density divided by the line width in units 
of velocity. The most straightforward method of approach to the problem is to solve 
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for the ratio of the line strengths. The ratios of 163 /im to 119/xm and 84 /xm to 
119 /j,m line radiation are what we shall consider. 

Figures 6.3 to 6.7 show the results of the runs. Displayed is the line flux ratio 
as a function of H 2 density at several selected temperatures. Each graph is for a 
different column density of OH. In figures 6.8 to 6.10, the plots are of the ratios as 
a function of density at several selected column densities. Each of these graphs is 
for a different temperature. 

One should note several trends in these data. For low column densities, 
one would expect the optical depths to be small. The resulting fluxes are then 
proportional to the ratio of the optical depths (see equation 6.8). For the 163:119 fim 
ratio, we observe the following behavior. For low hydrogen densities, we will have 
few collisions, and the radiative processes shall dominate. The flux ratio will then 
depend solely on the column density, since the kinetic temperature does not affect 
radiative processes. We will assume throughout the rest of this chapter that the 
OH line width is 50km s -1 (Crawford et al. 1986). It should be remembered that 
the column density parameter for the radiative transfer code is the actual column 
density divided by the line width. 

For low hydrogen densities, as the column density increases, we note that 
the 163:119 jxm flux ratio increases as well, until it reaches a stable value of «0.48. 
While it is increasing, the optical depth of the 163 nm line is also increasing. The 
119/xm line becomes optically thick rather quickly, unlike the 163 ^m line. 

When the number density is very high, then we expect the collisions to 
dominate in determining the steady-state level populations. The final 163:119 (j,m 
flux ratio should then be determined by the temperature and the column density 
of OH. The column density comes into play in determining the line optical depth, 
and hence the line flux as in (6.8). As all of the lines become optically thick, then 
the ratio becomes fixed for a given temperature. 

As the temperature increases, we see that the trends suggested by our opti- 
cally thick, LTE model in figure 6.2 are actually obtained. The 163:119 ratio 
decreases as the temperature increases while the 84:119 /xm line ratio increases. The 
actual ratios are quite close to those in our optically thick, LTE model, which is 
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expected. 

In between the two density extremes, we have a region of complicated be- 
havior. For example, if we examine figure 6.5, where the OH column density is 
5x 10 16 cm -2 , we see a trend with the increase in temperature. The peak 163:119 /xm 
flux ratio moves down as the gas becomes hotter. This peak also begins to occur at 
lower densities, although the density at the peak is weakly coupled to the temper- 
ature. The reason for the behavior entails looking at the different populating and 
de-populating processes in detail. At low temperatures, the collisional processes still 
dominate in the 119 /im transition. However, radiative processes play the major role 
for the 163 /xm line. As the gas warms up, collisional excitation and de-excitation 
become important to the 163 /xm line, and the collisional processes bring this line 
flux down. The collisional processes tend to suppress the 119 /xm line. Since for low 
temperatures, only radiative processes contribute to the 163 /xm flux while collisions 
dominate in the 119 /xm flux, we get this enhancement in the 163:119 /xm ratio. 

We do not find this behavior for the low OH column density results (figures 
6.3 and 6.4). That is to say, for a particular temperature, the ratio will either 
monotonically increase or decrease, but it will not have the characteristic bumps 
shown in figure 6.5. This monotonic behavior is the result of the 119 and 163 /xm 
lines being dominated by the collisional processes at the same densities. 

All of the above arguments also hold for the corresponding trends in the 
84:1 19 /xm ratio. 


6.5 Discussion 

We are now in- a position to fit the data with a model based upon the results 
presented above. From the OH data, we know that the 163:119 /xm flux ratio is on 
the order of 1.6 to 2.5, depending on the strength of the 163 /xm lines. Most likely, 
the ratio is nearer to the lower limit. The 84:119/xm flux ratio is roughly equal to 
1.6. 

We see immediately from the results of our numerical runs that we need to 
suppress the collisional contributions to the 163 /xm transition in order to have such 
a high flux ratio. In fact, it is only for temperatures less than 50 K that we can 
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FIGURE 6.3a Flux ratio of 163 to 119^01 OH line radiation as a function of 
density for several temperatures. The OH column density is 5 x 10 14 cm -2 . We 
assume the [e~]/[H] abundance ratio to be 5 x 10 -6 , and use the dust continuum 
radiation field as found by Pipher, Duthie, and SavedofF (1978). 
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FIGURE 6.3b Same as figure 6.3a, except the line ratio is for the 84 to 119 
line radiation. 
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FIGURE 6.4a Same as figure 6.3a, with an OH column density of 5 x 10 15 cm 2 . 
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FIGURE 6.4b Same as figure 6.3b, with an OH column density of 5 x 10 15 cm -2 . 
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FIGURE 6.5& Same as figure 6.3a, with an OH column density of 5 x lO^cm 
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FIGURE 6 . 5b Same as figure 6.3b, with an OH column density of 5 x 10 18 cm 
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FIGURE 6.6a Same as figure 6.3a, with an OH column density of 5 x 10 17 cm 2 . 
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FIGURE 6.6b Same as figure 6.3b, with an OH column density of 5 x 10 17 cm 2 . 
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with an OH column density of 5 x 10 18 cm -2 . 
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FIGURE 6.7b Same as figure 6.3b, with an OH column density of 5 x 10 18 cm -2 . 
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FIGURE 6.8a Flux ratios for 163 to 119 /im line radiation as a function of density 
for various OH column densities. These data assume a gas temperature of 30 K. 
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FIGURE 6.8b Same as figure 6.8a, except the line ratio is for the 84 to 119 /xm 
line radiation. 
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FIGURE 6.9a Same as figure 6.8a, except the gas temperature is 100 K. 
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FIGURE 6.9b Same as figure 6.8b, except the gas temperature is 100 K. 
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FIGURE 6. 10a Same as figure 6.8a, except the gas temperature is 1000 K. 



193 



FIGURE 6. 10b Same as figure 6.8b, except the gas temperature is 1000 K. 
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have such a high 163 /xm flux. 

By using the two ratios, we can determine the temperature of the 0H- 
emitting region to be between 35 and 40 K. No other temperatures give a reasonable 
fit. At these temperatures, the best fit is obtained for an OH column density of 
1.0 x 10 17 cm -2 and a molecular hydrogen density of 4 x 10 7 cm" 3 . These values are 
highly indicative of cool, post-shock gas. If the gas were pre-shock, then the density 
would have to be lower by a factor of 100. That the OH FIR emission originates in 
this cool gas confirms the suspicion of Melnick (1985). 

It is possible for shocked gas at a temperature of 2000 K to cool rapidly to 
temperatures of 40 K or less. This cooling is done primarily by the CO emission 
lines. From McKee et al. (1982), we can estimate the cooling rate of the gas as a 
function of temperature. A rough analytic fit shows that the cooling is proportional 
to T 1 * 5 . If a region is heated to 2000 K, then in a matter of 100 years or less, the 
temperature of the gas (hydrogen plus CO) will be cooled to 40 K by the CO alone. 
The shock has been moving across the BN-KL complex for several thousand years, 
leaving behind post-shock gas, most of which has cooled down to 40 K or less. 
Therefore, the temperature that we have deduced for the OH-emitting region is 
quite consistent with the picture of the post-shock region. 

The expected flux at 119 /zm from this range of values is about 3.2 times that 
actually observed, implying that the beam filling factor is about 31%. If we assume 
an [OH]/[H 2 ] ratio of 2.1 x 10 -7 (Draine and Roberge 1982), then from the density 
determination, we find that the OH region is roughly 10 16 cm deep. 

The CO data suggest that the pre-shock gas has a density of ~10 5 cm -3 , 
which rules out the pre-shock gas being the source of the OH radiation. The ex- 
pected temperature of the OH suggests that we are in a cool part of the post-shock 
g as. The infrared CO transitions imply that the post-shock gas just downstream 
from the shock front is roughly 750 K with nu 2 «2x 10 6 cm~ 3 . Simple pressure 
equilibrium dictates that a 40 K gas will have a density of roughly 4 x 10 7 cm -3 . 
It is unlikely that the density in the OH region is much greater them this estimate. 
No observations to date suggest higher densities in or near BN-KL, except on very 
small scales (Genzel et al. 1982). 
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The above arguments, coupled with the model fit to the data, suggest the 
following values for the physical parameters of the OH emitting region. 

T gM « 37 ± 3K 

nn 2 » 4 ± 2 x 10 7 cm -3 

Nq jj w 1.0 i.4 x 10 17 cm 2 

g( « 0.36 ± 0.04 

£>OH « 1.1 ±0.6 x 10 1 ® cm , 

where gf is the filling factor and £?qh I s the length of the OH-emitting region. This 
assumes that the OH lines are 50 km s -1 wide. 

From this fit to the data, we can estimate the expected line fluxes for all of 
the infrared transitions among the 14 lowest energy levels of OH. These are listed 
in table 6.1. Most of the infrared fluxes are close to 10 17 W cm -2 , well within the 
detectability of present day instruments. The 13.4 GHz line flux that is estimated 
from this set of parameters is within a factor of 2 of that actually observed. ' 

We also note that in the Orion Nebula, OH masers are known to be present. 
These level inversions are found to occur for temperatures substantially greater them 
the 40 K region that we apparently are observing with the FIR OH emission lines. 
According to conventional thinking, Hll /OH masers are thought to arise in very 
compact dense areas of the molecular cloud surrounding the associated H II region 
(Genzel 1985; and references within). Since these areas are small, their filling factor 
should be small, and we don’t expect them to significantly add to the FIR OH 
emission. OMC-1 is thought to be clumpy (Wynn-Williams et al. 1984; Geballe et 
al. 1986), and most likely the OH masers are in these clumps. 

We can calculate the amount of hydrogen that is associated with the OH, . 
and we find it to be about 50 M©. The OH abundance ratio could very well be in 
error by a factor of 2 or 3, since the abundance is dependent on the initial [0]/[H] 
ratio in the pre-shock gas and on the temperature and density history of the gas 
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Table 6.1 Predicted Infrared OH Line Fluxes for OMC-1 

Transition 


Wavelength 

Flux 

(10" 17 W cm' 2 ) 

2lI 3/2 ^ = I “* ! 


119.44 Atm 

0.81 



119.23 (im 

0.83 

2 n 3/3 .' = i-f 


84.60 ixm 

1.2 



84.42 Atm 

1.2 

2 iW = !-i 


65.28 Atm 

1.1 



65.13 

1.1 

2 n 1/2 J=|^J 


163.40 fi m 

0.87 



163.12 Atm 

0.90 

2n i/ 2 J = | -*■ § 


98.74 Atm 

1.9 



98.72 Atm 

1.9 

2 ni/2 J = i -* 2 n 3/2 j = 

3 

I 

79.18 /tm 

2.2 



79.12 Atm 

2.2 

2 Hi/2 J = | -* 2 n 3/2 J = 

3 

2 

53.35 Atm 

1.8 



53.26 Atm 

1.8 

2n i/2 j = 1 -> 2 n 3/2 j = 

3 

2 

34.63 Atm 

0.39 



34.60 /xm 

0.39 

2 n J/2 J = i - 2n s/2 J = 

5 

2 

96.37 Atm 

0.90 



96.31 Atm 

0.89 

2 n,/2J = |- 2 n 3/ 2J = 

5 

2 

48.82 Atm 

1.2 



48.70 Atm 

1.1 

2 n, /2 J = |- 2 n 3/ 2J = 

7 

2 

115.39 Atm 

0.030 



115.15 Atm 

0.031 
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in relation to the chemistry. Because of this, it is quite likely that the abundance 
may be greater than what we assume it to be (2.2 x 10 -7 ) , and this will help drive 
down the amount of hydrogen associated with the OH. 

We should note that the mass of hydrogen in the shock region is roughly 
20 M© if the average density is 10 5 cm~ 3 . The average density is most likely more, 
since the CO data suggest densities an order of magnitude greater than that. This 
leads us to believe that a mass of 50 M© is consistent with the overall picture of the 
BN-KL region. 

We propose a one-component model for the FIR OH emission lines from 
OMC-1. This model suggests that the bulk of the FIR line emission comes from the 
cool post-shock gas that has a kinetic temperature of 37 K. The molecular hydrogen 
density is roughly 4 x 10 7 cm -3 and the OH column density is on the order of 
1 x 10 17 cm -2 . This implies that this cool post-shock gas is roughly 10 16 cm deep, if 
the [OH]/[H 2 ] ratio is 2 x 10~ 7 . The beam filling factor for a one arc-minute beam 
is about 0.36. 

These parameters fit within the current framework of thinking concerning 
the region of OMC-1 containing BN-KL and IRc2. However, one should view the 
above values of the physical parameters as merely suggested by the data. Better 
absolute values of the fluxes are desired, as well as the acquisition of other infrared 
transition lines of OH. More accurate collisional cross-sections are also needed, 
as well as measured values of cross-sections between different excited levels. The 
experimental values for the H 2 -OH cross-sections as reported by Andresen, Hausler, 
and Lulf (1984) differ significantly from those reported by Schinke and Andresen 
(1984), although both are based upon the same data set. 

The radiative transfer code that was used has its own limitations, as well. 
With it, we may determine the average effective temperature and densities in a 
uniform cloud. As of yet, no structure of the region is incorporated into the program. 
We do not take into account any temperature gradients that may exist, or density 
gradients or chemical or shock evolution. We hope to incorporate these at a later 
time. 

One final note on the OH data. The measured values of the 84.4 and 84.6 /zm 



198 


lines are comparable. However, the former line is nearly coincident with a CO ro- 
tational line (J = 31 — ► 30), from which we expect to measure a flux roughly 
comparable with what has been measured in that spectral region. Current instru- 
mental resolution cannot differentiate between the two lines. The crucial question 
concerns which molecule we actually saw at 84.4 /zm: CO or OH? If the OH is opti- 
cally thick, then there is a good chance that it will absorb the CO-emitted photons, 
and suppress the contribution from the CO. For the present set of parameters, the 
optical depth at line center of the OH is w 12. In the OH wing where the CO will 
be, the OH optical depth will still be much greater than unity. Most likely, the 
CO emission is simply being absorbed by the OH, and thus the measured 84.4 (im 
line flux is a measure of the OH line flux, as originally proposed by Viscuso et al. 
(1985b), and argued against by Watson et al. (1985). 

The results of the numerical method need to be reconciled with those found in 
chapter V. This reconciliation is based upon the failure of the back-of-the-envelope 
calculation to adequately treat the OH data from OMC-1. This failure has several 
sources. The first source comes from the assumption that the collisional excitation 
rates dominate. This may be true for the 119 and 84 /zm lines, but the 163 /xm line 
is radiatively pumped by the radiation field from the dust. This was not properly 
taken into account in the analysis in chapter V. Not all of the radiative de-excitation 
processes were taken into account, as well. Because of these competing mechanisms, 
we were unable to accurately calculate the level populations in chapter V. One of 
the most important reasons for the failure of the back-of-the-envelope analysis is 
that the radiative processes occur on a local scale. That is to say that due to the 
velocity gradient, a photon cam be absorbed and/or reemitted only in a localized 
region, amd not throughout the length of the cloud. For these reasons, it is only by 
using numerical radiative tramsfer techniques that we can begin to correctly amalyze 
the Orion data. 

That is not to say that the analysis presented in the previous chapter is not 
useful or cannot be accurate. That analysis is certainly powerful in its simplicity 
and ease of physical understanding. It is possible to gain insight into the workings 
of the different and competing excitation and de-excitation mechanisms working 
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on the molecule of interest. If the region to be studied has low temperatures (to 
keep the number of populated levels small) and a simple structure (no imbedded 
sources of radiation, for example), then the back-of-the-envelope calculations should 
be fairly accurate. 


6.6 Summary 

We have used the escape probability formalism to develop a radiative transfer 
code with which we have been able to fit the FIR OH emission line data from OMC- 
1. This fit confirms the suspicions of Melnick (1985) that the OH-emitting region is 
not associated with the hot, post-shock gas as previously believed (Viscuso 1985a), 
but is instead part of the cool ( T « 40 K), post-shock gas further downstream from 
the shock front. This is suggested by the fact that the 163 pm flux is enhanced 
relative to the 119 /zm line emission. Radiative transfer models fail to predict the 
anomolous 163 jxm flux for gas temperatures greater than 45 K. 

The mass of cool, post-shocked gas is on the order of 50 Mq» This much cool, 
processed gas can exist because the CO is a very efficient coolant in the ISM. The 
gas should be able to experience a rapid cooling from 2000 K to 40 K in 100 years 
or less. 

The molecular hydrogen density is found to be about 4 x 10 7 cm -3 , and' the 
OH column density is roughly l x 10 17 cm -2 . The beam filling factor is 36%. These 
findings are not in complete agreement with the analysis presented in chapter V. 
The densities are comparable, but the similarities stop there. 

Much more work will have to be done in order to better understand the OH 
emission from OMC-1. Better absolute calibration of the line fluxes is necessary and 
we also need to expand the present data set to include higher spectral and spatial 
resolution. More observations of other rotational transitions of OH in the infrared 
will also help us to further constrain the physical parameters of the OH-bearing 
region. 

Analysis similar to the above, using numerical methods, proves to be ex- 
tremely helpful in determining the section of parameter space where the model is 
consistent with the data. Analysis such as presented in chapter V, is important in 
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obtaining a basic understanding of the data, and for determining, as a first guess, 
the physical parameters of the region. It is only with mathematically more powerful 
approaches such as numerical radiative transfer codes that we can begin to bring 
the whole picture together. Methods such as these should prove invaluable in the 
analysis of spectral data from ground-based, airborne, and space-based observato- 
ries. 



VII. Future Efforts 


7.1 OH in OMC-1 

The OH data available for OMC-1 thus far suffer several serious weaknesses. 
The uncertainty in the absolute fluxes makes determining the 163:119 jxm and 
84:119 /xm flux ratios difficult. There are problems with beam sizes as well. Since 
the diffraction-limited beam size is on the order of 1 ; , which corresponds to a linear 
size of 5 x 10 17 cm at the distance of Orion, we are unable to determine the true 
size of the region. The data indicate that the beam filling factor is on the order of 
0.35, suggesting that either the OH-emitting region is clumpy or, also likely, that it 
extends as far as the shock front, which has a size of roughly 2 x 10 17 cm (Genzel 
1985). 

Another problem with the OH data is the lack of high spectral resolution. 
The data collected by our group were taken with instrumental resolutions of 4000 
or so. This is not enough to provide us with adequate spectral line profiles, since 
higher resolution results on the 163.12 ^m, 119.23 /mi, and 119.44 /xm profiles by 
the Berkeley group suggest linewidths of «50km s -1 wide (Crawford et al. 1986; 
Watson et al. 1986). 

The last problem with the OH data set is the lack of observations in other 
transitions. Those that cross the ladders at 53 /xm ( 2 IIjy 2 J = § — ► 2 n 3 / 2 J = §), 
79 /xm ( 2 n 1/2 J = j — ♦ 2 n 3 / 2 J = |), and 96 /xm ( 2 IT 1/ / 2 ./ = §—> 2 n 3 / 2 J = \) are 
essential in order to understand the population structure of these lower levels (see 
figure 4.3). The lack of reliable H 2 -OH collision cross-sections also cause problems 
with the understanding of the OH data. 

Future plans should include higher spectral and spatial resolution in the 
observations. The high spectral resolution can be obtained with Fabry-Perot in- 
struments such as the Hinge discussed in chapter III and the Berkeley group’s 
FPS (Watson 1982). High spatial resolution awaits the completion of the 3-meter 
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class Stratospheric Observatory for Infrared Astronomy (SOFIA), which is likely to 
become operational in the early 1990’s if NASA approves funding. The European 
Space Agency’s Infrared Space Observatory (ISO) and NASA’s Space Infrared Tele- 
scope Facility (SIRTF) will also be launched in the next decade. These facilities 
will be able to provide high spectral and spatial resolution. 

It is important to be able to map the spatial extent of the OH. We suggested 
in chapter VI that the OH masers are found in compact, high-density regions with 
temperatures of 100 K or more. These regions axe not expected to contribute much 
to the lower level rotational transitions which are dominated by the warm (40 K) 
OH, but they will contribute most of the flux for transitions between highly excited 
levels. It should be possible to separate these warm and hot components by looking 
at their spatial extent (the hot component should be confined to or near the shock 
front while the warm component should be more extended) as well as by the flux 
ratios between the high J rotational transitions and the lower J transitions. 

We hope to be able to obtain high-resolution data on the cross-ladder tran- 
sitions mentioned above within the next year or two. The higher spatial resolution 
must await the next generation of infrared observatories, namely SOFIA or a 3- 
meter balloon-borne system. 


7.2 CO in OMC-1 

The CO data are in somewhat better shape than the OH data for OMC-1. 
The CO data, as can be seen in figure 5.8, cover transitions from levels ranging in 
excitation temperature from 5 K to 3000 K. The recent fits to the CO data have, 
until now, been based upon the H 2 -CO collisional cross-sections of McKee et al. 
(1982). These cross-sections were calculated by using the He-CO potential energy 
surface (PES), and then multiplying the resulting de-excitation rates by a factor 
reflecting the difference in relative velocity reflecting the reduced mass difference 
between He and H 2 . 

This procedure neglected the fact that for large A J transitions, the difference 
in the reduced mass results in an overestimation of the cross-sections by factors of 10 
or more (Green 1986, private communication; Schinke et al. 1985). Schinke’s group 
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has used a more accurate H 2 -CO PES to calculate the collisional excitation rates. 
We hope to employ these new rates to analyze the CO data for OMC-1 in order to 
determine the temperature and density of the CO-emitting region. Since the earlier 
rates were overestimated for the high J levels, we predict that the temperature 
and/or density of this region will increase. 

7.3 Future Objectives for the Hinge 
7.3.1 Airplane-Based Observations 

It is now possible to obtain the high spectral resolution needed for the OH 
observations as mentioned above. Since the Hinge has a resolution of 15 km s -1 , 
this will be more than adequate to determine useful line profiles, which will help 
greatly in the interpretation of the dynamics of the region. This is also true for the 
CO lines in the FIR. 

There is also interest in observing other molecular species which are predicted 
to be somewhat abundant by chemical models. These include CN and CS, which are 
products of the complicated carbon, nitrogen, and sulfur chemistry cycles (Prasad 
and Huntress 1980; Mitchell 1984). 

The needed molecular data for these two species are found in the following 
literature. For CN, the energy levels can be calculated by using the Dunham coeffi- 
cients obtained experimentally by Skatrud et al. (1983) and Cerny et al. (1978). The 
dipole moment of CN is reported by Thomson and Dalby (1966) to be 1.47 D, and 
is supported by the theoretical results of Lavendy, Gandara, and Robbe (1984). For 
CS, the Dunham coefficients and rotational constants are found to be in excellent 
agreement among four different groups (Bogey, Demuynck, and Destombes, 1982; 
Bergeman and Cossart, 1981; Bustreel et al. 1979; and Todd and Olson, 1979). The 
dipole moment for CS is reported by Winnewisser and Cook (1968) to be 1.958 D. 

To find the expected flux from OMC-1, we need to calculate the specific 
intensities for CN vs. those for CO. One finds that the CN intensities are a factor 
of almost 165 times higher than those for CO. This is, in large part, due to the 
higher dipole moment of the former (1.47 vs. 0.112). The population levels for 
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the two molecules are almost identical, so the intensities basically differ by their 
A coefficients. Assuming that the collisional cross-sections for CN and for CO are 
similar, then we should also expect that the CN transitions for J ~ 16 are density 
independent, as they are for CO. It is therefore only necessary to find the column 
density of CN and calculate the expected fluxes in the lines. 

For CS, the excitation temperatures for J ~ 38 (the FIR transitions) are ~ 
1700 K. We therefore don’t expect a significant population in these levels. However, 
due to the large degeneracy and the large value of /x, the specific intensities are 
still factors of 40 higher than CO transitions at these wavelengths. In fact, if 
CS is in LTE at 750 K, the maximum intensities occur at these far-infrared (FIR) 
wavelengths. This is also true for CO and for CN, making the FIR lines the best 
lines for detecting these molecules in dense regions where the level populations can 
approach LTE. 

The CO observations that have been made toward the Orion Molecular Cloud 
may fair outnumber the stars within this vast complex. Towards BN-KL, the col- 
umn density for CO is 4.1xl0 17 cm -2 (Viscuso 1985) which agrees with the values 
determined by Koepf et al. (1982) and Watson et al. (1985). These values for Nqo 
yield fluxes for the FIR transitions of CO on the order of several xlO -17 W cm -2 , 
which are easily detected. 

For CN, the observations suggest that the column density toward OMC-1 
is ~ 1 x 10 15 cm -2 (Wooten et al. 1982; Churchwell 1980; Turner and Gammon 
1975; and Jefferts, Penzias, and Wilson 1970). This yields fluxes on the order of 
3.4 x 10~ 17 W cm -2 for the FIR lines ( N ~ 18). It should be noted that the CN 
abundance in Orion and other regions is about 6 orders of magnitude higher than 
chemical equilibrium models predict ( Wooten et al. 1982). The CN abundance may 
also vary among sources due to time-dependent chemical processes as discussed by 
Allen and Knapp (1978). If this is the case, then CN observations are important in 
determining chemical processes that are occuring in molecular clouds. This could 
also help in explaining the apparently non-equilibrium values of CN. 

The CS observations indicate a column density on the order of 1.7 x 10 14 cm -2 
(Goldsmith et al. 1980; Liszt et al. 1974; Penzias et al. 1971). From these densities, 
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the CS FIR flux will be ~ 3 x 10 -18 W cm -2 . This will be detectable for NEP’s 
better than 2xl0 -14 W Hz -1 / 2 with several seconds of integration. Padman et al. 
(1985) find a CS column density of 7xl0 13 cm -2 , but as they point out, the CS 
emitting regions may be clumpy and this would tend to enhance the density. 

We should note that it is not clear that these molecules will approach LTE in 
their level populations. At the temperature and density of the shocked gas in OMC- 
1, the CO is not in LTE, although it is close to it for the low J values. It should 
be noted that the charged particle collision excitation mechanism is proportional to 
/i 1 - 5 (Goss and Field 1968), so that charged particles will begin to play a major role. 
Collisions with C II and electrons will help bring the CN and CS closer to LTE since 
these molecules have a much larger dipole moment than CO. The critical electron 
density, the ratio of the spontaneous emission rate out of a particular level to the 
collisional excitation rate into that level, is ~ 10 3 cm -3 for CO, while for CN, this 
value is down by a factor of 100, to about ~ 10 cm -3 . For CS, the A coefficient 
is larger due to the higher J number (CN and CO are nearly identical in energy 
strucutre) , and so there isn’t as much of an advantage in the higher dipole moment 
reducing the electron critical density. 

The CN and CS FIR emission lines will be able to tell us much about the 
chemistry of the shock region in Orion. Both species exhibit high velocity wings 
in their profiles of 30 km s -1 or more. They should be relatively bright and their 
detection would constitute the first FIR astronomical detection of these lines. 

7.3.2 Space Missions for the Hinge 

As was mentioned in chapter III, the Hinge is well-suited for space missions. 
Not only can one obtain high spectral resolution, but the package is very small 
and lightweight and has only one moving part. This makes the Hinge ideal for 
the SPARTAN and EUREKA payloads. The ruggedness of the system, relying as 
it does on a rotation of a grating about one axis without the plates moving, also 
assures better survival through the rigors of launch. 

The Hinge should be amenable to redesign for operations at other wave- 
lengths, either shorter or longer. We are looking into the possibility of constructing 
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an etalon for use at shorter wavelengths around 63 and 88 /j, m, where neutral and 
ionized oxygen fine-structure lines contribute significantly to the cooling of the in- 
terstellar medium (Stacey 1985). 

With the Hinge or a similar instrument in space, it would be possible to do 
spectral line surveys of selected areas of the sky, or possibly even of the entire sky, 
depending on how large a cryogen supply can be provided. The recent success of 
IRAS points to the possibility of making significant advances with high-resolution 
infrared spectroscopy from space. 
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